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The Two Higgs Doublet model and its phenomenological implications are dis- 
cussed. A brief survey on the present status of this model is given. In particular. 
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we concentrate on the Two Higgs Doublet Model with Flavor Changing Neutral 
Currents. First, we develop some new parametrizations of the 2HDM type III, in 
such a way that the relation of it with models type I and type II become apparent, 
based on two of these parametrizations we get some bounds on the mixing vertex 
involving the second and third lepton generations, as well as some lower bounds 
on the mass of the pseudoscalar Higgs boson; such bounds are obtained from the 
g — 2 muon factor. 

Further, by using a parametrization in which one of the vacuum expectation 
values vanishes, we constrain some lepton flavor violating vertices, assuming that 
the lightest scalar Higgs mass m/jO is about 115 GeV and that the pseudoscalar 
Higgs is heavier than h^. Specifically, based on the g—2 muon factor and the decay 
width of /i — > 67, the following quite general bounds are obtained: 7.62 X 10^^ < 
< AAA X 10"^ < 2.77 x 10"^''. Additionally, based on the processes 
T ^ /i7, and r HHH, bounds on ^T-^ and are also gotten, such constraints 
on these parameters still give enough room for either a strong suppression or 
strong enhancement on the coupling of any Higgs boson to a pair of tau leptons or 
a pair of muons. Moreover, upper limits on the decay widths of the leptonic decays 
T —>■ 67, and r — > eee are calculated, finding them to stay far from the reach of 
near future experiments. 

Finally, the Flavor Changing Charged Current decay — > Ve^^n is considered 
in the framework of the 2HDM type III as well. Since FCNC generates in turn fla- 
vor changing charged currents in the lepton sector, this process appears at tree level 
mediated by a charged Higgs boson exchange. From the experimental upper limit 
for this decay, we obtain the bound |^^e/^J?±| ^ 3.8 X 10~^ where mH± denotes 
the mass of the charged Higgs boson. This bound is independent on the other free 
parameters of the model. 

On the other hand, as an addendum to this work, we study the production 
and decays of top squarks (stops) at the Tevatron collider in models of low-energy 
supersymmetry breaking. We consider the case where the lightest Standard Model 
(SM) superpartner is a light neutralino that predominantly decays into a photon 
and a light gravitino. Considering the lighter stop to be the next-to-lightest Stan- 
dard Model superpartner, we analyze stop signatures associated with jets, photons 
and missing energy, which lead to signals naturally larger than the associated SM 
backgrounds. We consider both 2-body and 3-body decays of the top squarks and 
show that the reach of the Tevatron can be significantly larger than that expected 
within either the standard supergravity models or models of low-energy supersym- 
metry breaking in which the stop is the lightest SM superpartner. For a modest 
projection of the final Tevatron luminosity, L ~ 4 fb~^, stop masses of order 
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300 GeV are accessible at the Tevatron collider in both 2-body and 3-body de- 
cay modes. Wc also consider the production and decay of ten degenerate squarks 
that are the supersymmetric partners of the five light quarks. In this case we find 
that common squark masses up to 360 GeV are easily accessible at the Tevatron 
collider, and that the reach increases further if the gluino is light. 
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Chapter 1 
Introduction 



The Standard Model (SM) of particle physics has been very succesful in 
describing most of the smallest scale phenomenology known so far. How- 
ever, it possesses some problems whose solution could imply physics beyond 
its scope. In order to motivate the study of the two Higgs doublet model 
(2HDM) which is one of the simplest extensions of the SM, it is necessary 
to discuss two important issues. Perhaps the two most fundamental ideas 
from which the success of the SM comes are, (1) the extension of the gauge 
invariance principle as a local concept (inspired in classical electrodynamics) 
and (2) the implementation of the Spontaneous Symmetry Breaking (SSB) 
phenomenon. 

The introduction of local gauge invariance generates the so called gauge 
bosons as well as the interactions of these gauge bosons with matter (fermions) , 
and also among the gauge bosons themselves (the latter only when the gauge 
group is non abelian). On the other hand, the combination of local gauge 
invariance with SSB leads naturally to the Higgs mechanism which in turn 
generates the masses of weak vector bosons and fermions. I shall give a brief 
survey of both ideas emphasizing in the Higgs mechanism, since the Two 
Higgs Doublet Model (2HDM) is an extension on the symmetry breaking 
sector. 



1.1 Local gauge invariance 

It is well known from classical electrodynamics that Maxwell's equations 
are invariant under a "local" gauge transformation of the form.- — > An+ 
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d^X{x) where is the four-vector potentiaL On the other hand, taking the 
free Dirac Lagrangian 

£free = ^i^r^^~m)^ (1.1) 

We can see that such Lagrangian is invariant under the global phase shift 
^ — > e*^^'. Nevertheless, inspired in the local gauge symmetry in electrody- 
namics explained above, we could ask the following question, is it possible to 
extend the global symmetry and demand it to be local? if yes, what are the 
physical consequences?. It is straightforward to check that such locality could 
be accomplished by replacing the "normal" derivative by the "covariant 
derivative" = d/^+iqAi^, where Afj_ is a four vector field that transforms as 
Afj, — > A^ + d^X (x) when the local transformation of ^ — > exp {—iqX (x)) ^ is 
realized. The Lagrangian (1.1) becomes 



£ = ^{i-f^D^ - m)^ = *(i7^a^ - m)* - g^^.^T^'* = £free - J^A^ 

it is easy to see that this new Lagrangian is invariant under the combined 
transformations ^ — > e"**'^^^)^', — + d^X [x) . Now, if we interpret 
An as the four vector electromagnetic potential, then J'* is the four-vector 
electromagnetic current. To complete the Lagrangian of Quantum Electro- 
dynamics (QED) we just add the kinetic term that describe the propagation 
of free photons 



£qED — £free — J^A^ — -F'^'^ F^,, 

Fni/ = d^Ai, — diiA^ 

And the kinetic term for the free photons (which leads to the Maxwell's equa- 
tions) is also local gauge invariant. Therefore, the matter-radiation coupling 
has been generated from the imposition of "locality" to the gauge princi- 
ple. Additionally, to preserve locality we have introduced into the covariant 
derivative a four- vector field (the four- vector potential A^), which is called a 
gauge field, and also a parameter q which acts as the generator for the local 
group transformations U {x) — exp (— i^A (a;)). In this case, to analize the 
symmetries we have used the unidimensional rotation group in the complex 
space, this group is called U (1) , the group of unitary complex matrices 1x1. 
In the electroweak SM we use besides, the group of unitary matrices 2 x 2 of 
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determinant one (known as SU (2)), the latter is a non abelian group whose 
generators obey the Lie algebra of the rotation group in three dimensions 
SO (3) (i.e. they are isomorphic). After applying local gauge invariance to 
the whole electroweak group SU (2) x U (1) , four gauge fields appear and 
they in turn generate after some additional transformations the three weak 
vector bosons and the photon. 



1.2 Spontaneous symmetry breaking and Higgs 
mechanism 



Using local gauge invariance as a dynamical principle is not enough to pre- 
dict particle physics phenomenology since it leads to massless gauge bosons 
that do not correspond to physical reality. SM predicts that these vector 
bosons acquire their masses from a Spontaneous Symmetry Breaking (SSB) 
phenomenon explained below. 

Sometimes when the vacuum (minimum of the potential) of a system is 
degenerate, after choosing a particular one^, this minimum is not invariant 
under the symmetry of the Lagrangian, when the vacuum does not have 
the symmetry of the Lagrangian we say that the symmetry has been sponta- 
neously broken. When this phenomenon occurs some other massless particles 
called Goldstone bosons arise in the spectrum. However, if the Lagrangian 
posseses a local gauge symmetry an interrelation among gauge and Goldstone 
bosons endows the former with a physical mass, while the latter dissapear 
from the spectrum, this phenomenon is called the Higgs mechanism [1]. To 
explain the mechanism we shall use a toy model describing a couple of self 
interacting complex scalar fields (0 and 0*) whose Lagrangian is local gauge 
invariant 



^Quantum field theory demands the vacuum to be unique such that perturbation ex- 
pansions are calculated around that point 



£ \D'^<P\' - V (0» - ^F.^F'^^ ; V (0) 



lfj'\cpf+l^i<tfcpr (1.2) 



where 
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is a complex field and 

This Lagrangian has already the local gauge invariance described by the 
simultaneous transformations 

(P{x) ^ e-''^^^^U{x), A^{x) ^ A^{x) + d^X{x) (1.3) 

Observe that the imposition of locality generates the interaction of the com- 
plex scalar fields with a four vector field. If /i^ < 0, The potential V {</)) 
posseses a unique minimum at = which preserves the symmetry of the 
Lagrangian. However, if > 0, the Lagrangian has a continuum degenerate 
set of vacua (minima) lying on a circle of radius n/X 

any of them might be chosen as the fundamental state, but no one of them is 
invariant under a local phase rotation^. According to the definition made 
above, the symmetry of the Lagrangian has been spontaneously broken. 
Choosing a particular minimum: 

= ^ ^ ^ ; (02) = 

we say that the field 0i has acquired a Vacuum Expectation Value (VEV) 
(0i). It is convenient to introduce new fields 

r} = (t)i-v ; C = 02 

and expanding the Lagrangian in terms of these new fields we obtain: 



£ = 



1 a'^v'^ 



-2iqv {d,0 A^ + [q [ry {8,^ - i {d,r^)] A^ + vq^ {r)A^A^) 



^The set of all ground states is invariant under the symmetry but the obligation to 
choose (in order to set up a perturbation formalism) only one of the vacua, leads us to the 
breaking of the symmetry. 
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The particle spectrum consists of 



1. A field 1] with mass ^/2fi. 

2. A vector boson that has acquired a mass qv > by means of the 
VEV. 

3. A massless field ^ called a Goldstone boson. 

However, the Lagrangian above looks dissapointing because of a term of 
the form (5^^) which does not have a clear interpretation in the Feynman 
formalism. Fortunately, we are able to remove the unwanted would be Gold- 
stone field out, by exploiting the local gauge invariance of the Lagrangian. 
Writing Eq. (1.3) in terms of 4>i and 02 

^ = e^^(^V = [01 cos e (x) - 02 sin 6 (x)] + i [0i sin 6 (x) + 02 cos 6 (x)] 

where 9 (x) = —qX (x), and using 

9(x) = - arctan ( (1.4) 

we get 0' to be real^. The gauge field transforms as A'^{x) — Af^{x) + 
d^X{x). However, this gauge transformation does not affect the physical con- 
tent of Afj_{x) so we drop the prime notation out from it. Using the local 
transformations defined by (1.3) and (1.4), (i.e. this particular gauge), the 
Lagrangian reads 



i(a,77)(9^)-/.V 



1 o^v^ 
— F^^F"" + - — A. A" 



+ 



+ ^q'v (riA.A^) + ^v'A.A^ - Xptr,' - -r?^| + 

so we have got rid of the massless field C, and all its interactions, especially 
the "disgusting" term (d^C,) A^. On the other hand, we are left with a massive 
scalar field r] (a Higgs particle) and a massive four vector field A^ (a massive 
"photon"). 

^Observe that the invariance under a phase rotation of the complex field ^ e^^cf) 
is equivalent to the invariance under an SO (2) rotation of the real and imaginary parts 
(j)! 4'! cos — ^2 sin ^ ; (/>2 — > (/>i sin 6 + (j)2 cos 6 
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By making a counting of degrees of freedom we realize that one degree 
of freedom has dissapeared (a massless would be Goldstone boson) while 
another one has arisen (a longitudinal polarization for the four vector boson 
i.e. its mass). Therefore, it is generically said that the photon has "eaten" 
the would be Goldstone boson ^ in order to acquire mass. This result is 
known as the Higgs mechanism. Notwithstanding, it worths to emphasize 
that as well as the massive vector boson, the Higgs mechanism has provided 
us with an additional physical degree of freedom that corresponds to a scalar 
field describing the so-called "Higgs particle" . 

We can note that the Higgs mechanism is possible because of the conju- 
gation of both the local gauge invariance principle and the SSB. For instance, 
if we implement a SSB with a global symmetry, what we obtain is a certain 
number of (physical) massless Goldstone bosons, it is because with a global 
symmetry we do not generate vector bosons that "eat" such extra degrees of 
freedom. Technically, the number of Goldstone bosons generated after the 
symmetry breaking is equal to the number of broken generators (Goldstone 
theorem [2]). 

In SM, the Higgs mechanism creates three massive vector bosons {W^, Z) 
and one massless vector boson (the photon), as well as a Higgs particle which 
has not been discovered hitherto. 



1.3 The Higgs mechanism in the Standard 



The SM of particle physics [3], picks up the ideas of local gauge invari- 
ance and SSB to implement a Higgs mechanism. The local gauge symmetry 
is SU (2)^ X U {1)y and the SSB obeys the scheme SU (2)^ x U (1)^ ^ 
U (1)q where the subscript L means that SU (2) only acts on left-handed 
doublets (in the case of fermions), Y is the generator of the original U (1) group, 
and Q correspond to an unbroken generator (the electromagnetic charge). 
Specifically, the symmetry breaking is implemented by introducing a scalar 
doublet 



It transforms as an SU (2)^ doublet, thus its weak hypercharge is F = 1. 
In order to induce the SSB the doublet should acquire a VEV different from 
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zero 




(1.5) 



What is new respect to the toy model above, is that the original local sym- 
metry SU (2)^ X U is non abelian^. Its generators are Tj, Y corresponding 
to SU (2)^ and U (l)y respectively, the generators Tj are defined as 

where (jj are the Pauh matrices. Therefore, such four generators obey the 
following lie algebra 



When the symmetry is spontaneously broken in the potential (see below) 
the doublet acquire a VEV, we can see easily that all generators of the 
SU (2)^ X U (l)y are broken generators 



However, we can define an unbroken generator by the Gellman-Nijishima 
relation 



In such a way that the scheme of SSB is given by SU (2)^ x U (l)y — > 
U (1)q . According to the Goldstone theorem, the number of would be Gold- 
stone bosons generated after the symmetry breaking is equal to the number 
of broken generators (which in turn is equal to the number of massive gauge 
bosons in the case of local symmetries). Therefore, instead of working with 
four broken generators we shall work with three broken generators and one 
unbroken generator Q. This scheme ensures for the photon to remain mass- 
less, while the other three gauge bosons acquire masses. 



Let us examine the contributions that the doublet $ introduces in the 
SM. 



very important consequence of the non-abelianity of the gauge group is the genera- 
tion of self interactions among the associated gauge bosons, they appear when the kinetic 
term for the gauge bosons (Yang-Mihs Lagrangian) is introduced. 
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1.3.1 The Higgs Potential 



The Higgs potential generates the SSB as well as the self interaction terms of 
the scalar boson, the most general renormalizable potential invariant under 
SU (2)^ X U is given by 

= + (1.6) 

where /i^ and A are free parameters of the theory. Since A should be positive 
for the potential to be bounded from below, the minimization of the potential 
(1.6) leads to a SSB when /i^ < with the following scheme SU (2)^ x 
U (l)y — > U {\)q where Q is the electromagnetic charge. After the SSB the 
Higgs doublet acquires a VEV as in Eq. (1.5) from which the Higgs doublet 
gives mass to the Higgs particle and is able to endow the vector bosons and 
fermions with masses (see below). 



1.3.2 The kinetic term 

The kinetic term describes the interactions between scalar particles and vec- 
tor bosons, and provides the masses for the latter when the Higgs doublet 
acquires a VEV. The kinetic Lagrangian reads 



£Hn^{D,^){D>^^)^ ■ D, = d,-'j-YWi-ignW; (1.7) 

where WJ^ with i — 1,2,3 are the four- vector fields (gauge eigenstates), 
associated to the three generators i.e. the SU (2)^ symmetry. On the 
other hand, is the four-vector field associated to the Y generator i.e. 
the U (l)y symmetry, g and g' are coupling strengths associated to and 
respectively. After diagonalizing the mass matrix of the gauge bosons 
we obtain the following mass eigenstates 

= ^^^^^^ ; M^. = lg\^ (1.8) 

Ml - \v\g- + 9')- (1-9) 
4 cos"^ Ow 



J y sin 9w cos 9w J \ 



r3 



\ _ f cosdw - sin 9w \ f ^ij. 



;i.io) 



and the gauge boson A^j, (the photon) remains massless, it owes to the fact 
that this gauge boson is associated to the unbroken generator Q (the elec- 
tromagnetic charge) i.e. to the remnant symmetry U {1)q- 
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1.3.3 The Yukawa Lagrangian 

Finally, we build up a Lagrangian that describes the interaction among the 
Higgs bosons and fermions, the resultant SU (2)^ x U invariant La- 
grangian is called the Yukawa Lagrangian 

-£y - r/g^L$f/R + vfj^L^Dn + h.c. (Lll) 

where \E'i are left-handed fermion doublets, Ur, are the right-handed 
singlets of the up and down sectors of quarks^, rf^-^ are free parameters that 
define the vertices and consequently, the Feynman rules of the Lagrangian 
where i,j are family indices. The Yukawa Lagrangian yields masses to the 
fermions when the Higgs doublet acquire VEV. 

So in brief, Eqs. (1.6,1.7,1.11) describe the contribution that the Higgs 
sector gives to the SM. 



^For leptons we have the same structure except that wc do not have right handed 
singlets of neutrinos (corresponding to the up sector of leptons). 
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Chapter 2 

The Two Higgs Doublet Model 



2.1 Motivation 

Despite the SM has been very sucessful in describing most of the Elementary 
Particles phenomenology, the Higgs sector of the theory remains unknown so 
far, and there is not any fundamental reason to assume that the Higgs sector 
must be minimal (i.e. only one Higgs doublet). Therefore, we could wonder 
to know whether the Higgs sector is not minimal. Of course, evocating 
arguments of simplicity we may consider the next to minimal extension as the 
best candidate. The simplest extension compatible with the gauge invariance 
is the so called Two Higgs Doublet Model (2HDM), which consists of adding 
a second Higgs doublet with the same quantum numbers as the first one. 

Another motivation to introduce the second doublet comes from the 
hierarchy of Yukawa couplings in the third generation of quarks, the ra- 
tio between the masses of the top and bottom quarks is of the order of 
mt/rrib fa 174/5 ^ 35. In SM, the masses of both quarks come from the same 
Higgs doublet, consequently, it implies a non natural hierarchy between their 
corresponding Yukawa couplings. However, if the bottom received its mass 
from one doublet (say $i) and the top from another doublet (say $2)5 then 
the hierarchy of their Yukawa couplings could be more natural if the free 
parameters of the theory acquired the appropiate values. 

On the other hand, the 2HDM could induce CP violation either explic- 
itly or spontaneously in the Higgs potential (sec section 2.2.1). However, we 
shall restrict our discussion on a CP conserving framework. A recent com- 
prehensive overview of the conditions for the 2HDM to be CP invariant (or 
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violated) can be found in Ref. [46]. 

An extra motivation lies on the study of some rare processes called Flavor 
Changing Neutral Currents (FCNC). It is well known that these kind of 
processes are severely supressed by experimental data, despite they seem not 
to violate any fundamental law of nature. On the other hand, Standard 
Model (SM) issues are compatible with experimental constraints on FCNC 
so far [4, 5], with the remarkable exception of neutrino oscillations [8]. In 
the case of the lepton sector this fact is "explained" by the implementation 
of the Lepton Number Conservation (LFC), a new symmetry that protects 
phenomenology from these dangerous processes. However, if we believe that 
this new symmetry is not exact and we expect to find out FCNC in near 
future experiments, SM provides a hopeless framework since predictions from 
it, are by far out of the scope of next generation colliders [7]. This is because 
in the SM, FCNC are absent in the lepton sector, and in the quark sector 
they are prohibited at tree level and further supressed at one loop by the GIM 
mechanism [6]. So detection of these kind of events would imply the presence 
of new Physics effects. Perhaps the simplest (but not the unique) framework 
to look for these rare processes is the Two Higgs Doublet Model (2HDM). 
Owing to the addition of the second doublet, the Yukawa interactions lead 
naturally to tree level FCNC unless we make additional assumptions. On the 
other hand, because of the increasing evidence about neutrino oscillations, 
especial attention has been addressed to the Lepton Flavor Violation (LFV) 
in the neutral leptonic sector by experiments with atmospheric neutrinos [8]. 
Further, since neutrino oscillations imply LFV in the neutral lepton sector, it 
is generally expected to find out LFV processes involving the charged lepton 
sector as well, such fact encourage us to study the 2HDM as a possible source 
for these FCNC. More details in section 3.2.4 and in chapter 4. 

Finally, another motivation to study the 2IIDM is the fact that some mod- 
els have a low energy limit with a non minimal Higgs sector. For instance, 
at least two Higgs doublets are necessary in supersymmetric models and 
the so called 2HDM type II has the same Yukawa couplings as the Minimal 
Supersymmetric Standard Model (MSSM). In particular, if the supersym- 
metric particles arc heavy enough, the Higgs sector of the MSSM becomes 
a constrained 2HDM type II at low energies. SUSY models with two Higgs 
doublets could provide solutions to some problems of the SM such as the 
Higgs mass behavior at very high scales, the Planck and Elcctroweak scale 
hierarchy, the mass hierarchy among fermion families and the existence of 
masses for neutrinos and neutrino oscillations. 
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With these motivations in mind, let us go to see in what way the terms 
in the Higgs sector are modified by the introduction of a new Higgs doublet 



2.2 The contribution of the Higgs sector in 
the 2HDM 

As explained above, we introduce a new Higgs doublet that is a replication of 
the first one, so the Higgs sector includes two Higgs doublets with the same 
quantum numbers 

with hypercharges Yi = F2 = 1, in general, both doublets could acquire VEV 

so it is more convenient to parametrize the doublets in the following way 

/ 0+ \ f <f>2 \ 

^1= y hi+vj+igi j ; $2 = I h2+V2ei^+ig2 j (2-2) 

Now, we examine each part of the Lagrangian that couples to the Higgs 
doublets 



2.2.1 The Higgs potential 

Since the Higgs potential is the sector that determines the SSB structure as 
well as the Higgs masses, Higgs mass eigenstates and Higgs self-interactions, 
we start examining this part of the Higgs sector. Unlike the SM case, the 
Higgs potential corresponding to the 2HDM is not unique, and each potential 
leads to different Feynman rules, details below. 

In order to write the most general renormalizable Higgs potential compat- 
ible with gauge invariance, it is convenient to introduce a basis of hermitian, 
gauge invariant operators 

A = , B = ^1^2, C=^ ($I$2 + = Re ($1*2) , 

D = ~ (^$1$2 - *2*i) = Im (^$1$2) 

13 



and write down all possible hermitian bilinear and quartic interactions com- 
patible with gauge invariance 

Vg{^i,^2) = -iilA-i4B-iild-tJ,lD + XiA:' + X2B^ + X3d'' + X^D'' 
+X5AB + XqAC + XsAD + X7BC + XgBD + XwCD (2.3) 

This Lagrangian is much more complex than the SM one given by Eq. (1.6), 
since in the potential (2.3) we have fourteen free parameters. As we shall 
sec below, four new Higgses will be generated from it. The most general 
2HDM potential (2.3), contains interaction vertices that are independent on 
the mass matrix and the VEV's. Nevertheless, these interactions vanish if 
we assume that the Higgs potential holds a charge conjugation invariance 
(C— invariance), and the number of parameters reduces to ten 

V^($i,$2) = -lilA- lilB- lilC + \iA^ + X2P + XzC^ + X^D^ 

+X^AB + X(iAC + XjBC (2.4) 

It is important to say that at this step, charge conjugation invariance is equiv- 
alent to CP invariance since all fields are scalars^. Under charge conjugation, 
a Higgs doublet of hypercharge 1, transforms as $i — > e*"*$* where the 
parameters ccj are arbitrary. Consequently, under charge conjugation we ob- 
tain — > e*^"-'~"'^$j<l>j. In particular, if we choose ctj = aj the operator 

D reverse sign under C— conjugation, while the other ones are invariant^, 
leading to the Lagrangian (2.4). 

However, Lagrangian (2.4) could induce spontaneous CP violation [9]. 
There are two ways of naturally imposing for the minimum of the potential 

^The appearing of pseudoscalar couplings come from the introduction of parity violation 
in the theory. It is carried out by introducing left-handed fermion doublets and right 
handed fermion singlets which in turn produces scalar and pseudoscalar couplings, the 
latter are proportional to 75 and are responsible for parity violation. Since the Higgs 
doublets do not have chirality, their self couplings respect parity and therefore they behave 
as scalars. However as we shall see in section (2.2.3) when the Higgs doublets couple to 
fermions some pseudoscalar couplings arise and one Higgs boson behaves as a pseudoscalar. 

^Of course, we could have chosen a, — aj = ±7r, in whose case the operator C = 
Re ^$|$2^ if' the one that violates charge conservation. Additionally, any other choice 
for ai — aj is possible, and in general none parameter vanishes. However, taking into 
account that these phases must be fixed (though arbitrary), C— invariance would impose 
relations among the coefficients so that the number of free parameters is always the same 
(for instance (13 and (14 would not be independent any more). 
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to be CP invariant [10]. The first one consists of demanding invariance under 
a Z2 symmetry where $1 ^ $1, $2 ^ — ^2- The resulting potential, that 
will be denoted is 

= -i4A - lilB + XiA^ + AsS^ + AgC^ + + X^AB (2.5) 

and correspond to setting = Ae = A7 = in Eq. (2.4). If we permit a soft 
breaking term of the form —nlC, spontaneous CP violation occurs [29], in 
that case the potential reads 

Va = V^- ixlC (2.6) 

The other potential without spontaneous CP violation, results from im- 
posing the global symmetry $2 — e"^$2- This potential (called V^) reads 

V'^ = -filA - filB + XiA^ + A252 + A3 (C^ + + X5AB (2.7) 

and is obtained by using //| = Ae = Ay = and A3 = A4 in Eq. (2.4). Since 
we have a global broken symmetry for the potential V^, there is an extra 
Goldstone boson in the theory. 

Additionally, it is customary to allow a soft breaking term —ii\C in La- 
grangian (2.7), obtaining 

Vb^V^- ixlC (2.8) 

Neither or have spontaneous CP violation. In the case of Vb, we end 
up in the CP conserving case of the scalar potential considered in the Higgs 
Hunter's guide (see appendix A), which in turn is the one of the MSSM. 
Both of them contain seven parameters and lead to different phenomenology, 
however the kinetic sector and the Yukawa sector are identical for both po- 
tentials. The soft breaking term in Vb is a quadratic term and consequently 
it does not affect the renormalizability of the model, the complete renormal- 
ization scheme for the potential Va has been accomplished in Rcf. [11], the 
results for Vb are similar but changing appropiately the cubic and quartic 
scalar vertices. 

The Higgs masses and Higgs eigenstates are defined in terms of the pa- 
rameters /ij, Aj from the potential, and consequently depend on the potential 
chosen (see appendix A). When the mass matrix described in appendix A 
is properly diagonalized, we get the Higgs masses and Higgs mass eigen- 
states. From now on, we will consider the CP conserving case for which 
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both VEV can be taken real, in that case the Higgs sector consists of the 
following spectrum, two Higgs CP-even scalars {H^,h^), one CP-odd scalar 
(^4°), two charged Higgs bosons (i?^), and the Goldstone bosons 
corresponding to W^, Z respectively. The mass eigenstastes described above 
are obtained from the gauge eigenstates defined in (2.2) by the following 
transformations 



where 



cos/3 sin/3 \ ( (Pt \ _ ( G+ 

-sin 13 cos 13 )\<pt ) ~ \ H+ 

cos a sin a \ ( hi\ _ / 

— sin a cos a J y h2 J y 

cos P sin P \ f gi \ f 



^1 



— sin/3 cos/3 J \ g2 J \ A' 



tan [3 — — , sm /? 



) ^2.9) 



2 



and a is the mixing angle for the CP— even Higgs bosons, which is different 
for each potential (see below), tan/3 is a new parameter that clearly arises 
from the fact that both Higgs doublets could acquire VEV. In most of the 
2HDM's we shall see that new physics contributions are very sensitive to it. 

In the following we summarize the results obtained in appendix A, which 
are in agreement with [12, 17]. 

For the potential V in Eq. (2.4) the minimum conditions (Tadpoles at 
tree level) are 



U - J-b--^2.^ 2~ 



the Higgs masses and the mixing angle a, are written as 
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'>T^H+ = -lA + ' '^A" = -lA + ^ (-^4 + A5) vl 

h = ^J^Xlvf {Xivf + vlX+) + {Xlvf + Xlvf - 2i4X+) vf + /ij 

tan 2a = -— ; — r-^^ ^ 

(2Ai - A+) vl + 1^1 

For the potential Eq. (2.5), the minimum conditions are 

= Ti = vi {-fil + Xivl + X+vl) 

= T2 = V2{-iJ,l + X2vi + X+vl) (2.10) 

where A+ = | (A3 + A5) . Relations (2.10) lead to the following solutions 
i) 

2 X2^i — X+^2 2 A1/X2 — A+z^i 



A1A2 — A^ ^ A1A2 — A^ 



or ii) 



2 n 2/^1 

^2 = ; = T- 



for i) the masses of the Higgs bosons and the mixing angle a are given by 
mH± = -A3 {vl + vl) ; m\o = ^ (A4 - A3) {vl + v^) 



'^HO,hO = Xivl + X2V2 ± y {Xivf - X2vif + Avfv^Xl 

tan 2a = ^^T^ (2-11) 



and for ii) they are 



= -III + ^X^vl ; m\o = -l4 + ^ (A4 + A5) vl ; m^o = 2Ai^.l2) 
mlo = -/^2 + ^ (A3 + A5) ; tan2a = (2.13) 
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Finally, for the potential Vb Eq. (2.8), the minimum conditions are 



= Ti 



whose solutions are 



Ax — A2 i Zi 



2 (Ai - A+) (A2 - A+) 

A2 — Ai ih Z2 

2 (Ai - A+) (A2 - A+) 



(2.14) 



Zi = J(Ai-A2)'-4(Ai-A+)(A2-A+) 



(A+'f;2 - /if) (A+v^ - fj,l) - 



Z2 = W(Ai-A2)'-4(A2-A+)(Ai-A+) 



(A+V2 - ;xi) (Aiv2 - ;xf ) - 



The masses and the mixing angle a are given by 



mH± = -A3 {vl + vl) + 1^1 



vi + vi 



2 "^1 



V1V2 



2 _ 1 2 ^1+^2 



1 



Aif X + A2f2 + ^A'-s (tan (3 + cot /3) 



tan2Q; 



±W Ai^i — A2t'| + ^/ii (tan/? — cot/9) 
2t;it;2A+ - 



+ ( 2V1V2A+ - -nl 



Ai-yf - \2vl + (tan /? - cot /?) \ ■ ) 

Observe that a solution with one of the VEV's equal to zero is not pos- 
sible for the potential Vb- As explained before, each potential has different 
Feynman rules and consequently, leads to different phenomenology [25], a 
complete set of Feynman rules for these potentials could be found in Refs. 
[25, 46]. 

The potentials V'a^ and Vb are different because they differ in some cubic 
and quartic interactions. For example, the coupling h^H^H^ reveals some 
subtle aspects of the phenomenology of the potential; in terms of the A's it 
is given by 

gh°H+H- — 2t>2A2 cos^ /9 cosa + f2A3 sin a COS/9 sin/9 — f 1A5 cos^ /9 sin a 
— 2^1 Ai sin^ /3 sin a + t'2A5 sin^ (3 cos a — ViX^ cos a cos /3 sin j3 
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and coincides for both potentials and Vb', it is because this interaction 
does not involve A4 nor /ia which are the factors that make the difference 
between them. However, by writing the coupling in terms of the Higgs boson 
masses the result is different for each potential 

^mloj sin {a-P) 



we can resolve the puzzle by remembering that what really matters in per- 
turbative calculations is the position of the minimum of the potential and 
the values of the derivatives at that point. Now, the position of the vacuum 
(minimum) is different for each potential and the generation of masses come 
from the second derivative of the potential evaluated at this minimum, so the 
relation among the X's and the masses are different for and Vb, thus the 
coupling h^H^H~ in terms of physical quantities differs for each potential. 
It should be pointed out that even before being able to test the Higgs bosons 
self-couplings in order to discriminate among the potentials, we might see a 
signature of them in processes with Higgs boson loops, for example, the decay 
77 could be very sensitive to the difference among the self couplings of 
and Vb- 

Finally, it is also important to notice that the symmetry of the potential 
ought to be extended to the other Higgs sectors, this fact is particularly 
important to write the Yukawa Lagrangian to be discussed in section 2.2.3. 

2.2.2 The kinetic sector 

The kinetic Lagrangian (1.7) of the SM is extended to become 

£kin = (L'^$i)+(L"'$i) + {D^^2y{D^^2) (2.16) 

where the covariant derivative is defined by Eq. (1.7). This Lagrangian 
endows the gauge bosons with mass and provides the interactions among 
gauge and Higgs bosons. In contrast to the Higgs potential and the Yukawa 
Lagrangian (see sections 2.2.1, 2.2.3), the kinetic sector is basically unique 



{ghOH+H-)A 
{ghPH+H-)B 



9 



mw 

9 
mw 



2 cos(a + /3) / 2 

^ " ^ ^ sm2/5 
- \^hoj sin (a -/?) 
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because of the gauge invariancc''. Indeed, we can easily check that the ki- 
netic Lagrangian described by Eq. (2.16) is aheady invariant under charge 
conjugation as well as under the discrete and global symmetries described in 
Sec. (2.2.1), therefore the imposition of these symmetries does not produce 
any difference in the kinetic sector unlike the case of the potential. 

In order to expand the Lagrangian, it is convenient to work on a real repre- 
sentation for the generators, for which we make the assigment La = —iTa and 
double the dimension of the representation by means of the definition 



( Rec^t 
Im(j)'l 
Recj>l 



(2.17) 



From which we find the real representation by looking at the action of the 
initial generators (multiplied by —i), over the two dimensional representation 
of the doublet defined in Eq. (2.17). For example, for ti 



Li$i = -ITi^i = 



$1 = - 

2 2 



Im4>i 



iRe(f)i 



(2.18) 



And we extend with the same correspondence rule 

Li$i ^ I 



(2.19) 



Im(j)f 
-Re(j)\ J 

Finally, we look for a matrix Li that acting on (2.17) reproduces (2.19), 
we get 



(2.20) 





/ 











1 \ 


1 










-1 





2 







1 










V 


-1 








0/ 



Similarly we obtain for the other generators 



^However, for some potentials and Yukawa Lagraiigians it is possible to rotate the 
Doublets such that only one of them acquire VEV, as wc will sec later. In that case the 
kinetic term is basically the same but taking (say) V2 =Q i.e. tan /3 = 0. 
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(2.21) 



Prom which the covariant derivative in (1.7) reads 



d^ + gL,Wl + g'UW^ 



(2.22) 



Thus, we can expand the Lagrangian (2.16) by using the four dimen- 
sional representation for each doublet Eq. (2.17), and the four dimensional 
representation of the generators Eqs. (2.20), (2.21). 

Gauge fields 

The mass terms are obtained from the VEV's, the resulting mass matrix is 
given by 

1 ^ 

-MlW;Wl ■ Ml, = 2Y,{9aLaVk)\9bUvk) (2.23) 

k=l 

Where a,b = 1,2,3,4 indicates the gauge bosons (gauge eigenstates) 
corresponding to the generators r', Y respectively. After the diagonalization 
of (2.23) the mass terms and eigenstates for the Gauge bosons read 



± 



W. 



M' = 



" " ■ M? 



V2 



,2^ _ 



/ cos 9w — sin 9w 
y sin^^ cos^^ 



■ 2/ 2 I 2\ 



cos^ 9w 



(2.24) 
(2.25) 
(2.26) 
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1— 




1+ 


Z 


1+- 




0+ 




0++ 








0++ 








0+- 







Table 2.1: Assigments to the parity (P) and charge conjugation (C), quan- 
tum numbers for the Higgs and vector bosons before the introduction of 
fermions. When fermions are absent, P and C are conserved separately, (i.e. 
in the Higgs potential and kinetic Lagrangian). 

where 9w is the Weinberg mixing angle. We can realize that the expressions 
for the vector bosons masses coincide with the ones in SM if vf + V2 = 
v"^ (where v is the VEV of the Higgs doublet of the SM), and since v"^ = 
4M^/gf^ is a known parameter, we have at tree level the constraint v\-\-V2 = 

Since the kinetic sector provides the interactions of Higgs bosons with 
gauge bosons, let us discuss some interesting features of these interactions. 

Interactions in the kinetic sector 

We can obtain the interactions in the kinetic sector by expanding the La- 
grangian (2.16) in terms of the mass eigenstates of Higgs bosons and Vector 
bosons Eqs. (2.24, 2.26, 2.9). On the other hand, it worths to note that 
some interactions that at first glance should appear, are absent in this ex- 
pansion as is the case of A°W+H^~, and A^ZZ. Let us discuss shortly the 
origin of these missings. First of all we should emphasize that in the SM, 
the discrete symmetries C and P are preserved separately when fermions are 
absent. Then, we can assign a unique set of quantum numbers J^^ to all the 
bosons of the theory when fermions are ignored. The same argument holds 
for the two Higgs doublet model. These facts dictate the presence or missing 
of some interactions [17]. 

First of all, let us examine the assigment of these quantum numbers for 
the scalar and vector bosons in the table (2.1). 

The existence of the vertex ZH^H~ indicates that J^^ (Z) = l+~, the 



22 



existence of H^hPhP in the Higgs potential says that J^*^ {H^) = 0"^+; sym- 
metrically the existence of h^H^H^ tells us that J^*^ {h^) = 0^+. From 
these quantum number assigments we deduce from the vertex ZH^A^ that 
J^^ (AO) = 0+-. 

Furthermore, in the 2HDM the imaginary parts combine to produce 
and G° so both must have the same quantum numbers, consequently 
A^ and have the quantum numbers 0"*"". Owing to the assigment J^^ (A^) 
— 0~^~ the vertices A^li^+VT" and A^ZZ are forbidden. We should remark 
that despite the Higgs A'^ is usually called a pseudoscalar Higgs, at this step 
^4° is a scalar particle. However, we shall see later that ^4° behaves as a pseu- 
doscalar when it couples to matter (fermions), see section (2.2.3). Thus, the 
term pseudoscalar is not so proper for the Higgs A^ and the denomination 
as a CP— odd Higgs is more appropiate. 

Another way to explain the absence of vertices A^l^+l^" and A^ZZ is 
the following: These interactions come from the kinetic term (D^^) (D'*$)^ 
after replacing one of the 0° fields in Eq. (2.1) by its VEV, but in a CP 
conserving model the VEV is real while A° comes from the imaginary parts 
of the 0° fields as we can see from Eqs. (2.2) and (2.9). Therefore A^ cannot 
be coupled to any massive vector boson. By the same token, the couplings 
(if°, /i°) 77 are forbidden at tree level, since /i° cannot be coupled at 
tree level to massless vector bosons. Finally, another interesting argument is 
that A^W~^W~ and A^ZZ are prohibited by C— invariance, however, when 
fermions are introduced (by means of loops) C— invariance is no longer valid 
(though CP— invariance still holds in good aproximation) . Consequently, 
these vertices can appear at one loop level by inserting fermions into the 
loops. 

So in general, vertices with neutral particles only and one or two photons 
clearly vanish at tree level but can be generated at one loop. And same for 
a pair of gluons coupled to any of the neutral Higgs bosons. Gluon gluon 
fusion generated by one loop couplings (A°, /i°) gg, is a very important 
source for Higgs production in Hadron colliders. 

Additionally, there are other couplings whose absence can be explained by 
the quantum number assigments in table (2.1). For instance, the coupling of 
Z to non identical Higgs bosons is allowed only if such scalars have opposite 
CP numbers. Consequently, the couphngs ZA^H^, ZA^hP are allowed while 
ZH^hP is forbidden. On the other hand, the couplings of Z to identical Higgs 
bosons arc forbidden by Bose-Einstein symmetry. Moreover, H^W^^ is 
prohibited by conservation of electromagnetic current. 
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As for the vertices H^W^Z, they do not appear at tree level in the 
2HDM. We should note however, that they are not forbidden by these quan- 
tum numbers (they appear in some triplet representations). In the case of 
multi-Higgs doublet models (including the 2HDM) they are prohibited at 
tree level because of the isospin symmetry. Notwithstanding, such couplings 
are allowed at loop levels owing to the breaking of the isospin symmetry by 
the loop particles. One loop contributions to the H^W~Z vertex are studied 
in [18, 19, 20] and will be discussed in section (3.2). 

Finally, we should emphasize that when CP is violated, many of the 
missing vertices described above appear, though they are expected to be 
quite suppressed. 

Furthermore, it is interesting to notice that the couplings of the kinetic 
sector satisfy automatically some tree level unitarity bounds. Since partial 
amplitudes cannot grow with energy, cancellations to avoid violation of uni- 
tarity are necessary. For example, in the case of VlVl — > VlYl scattering, 
cancelation is possible in SM because the vertex is of the form gcj,"ww = 
gmw- When we have more than one doublet, this work does not have to be 
done by only one Higgs boson, instead we have the sum rule 

''^dlyv = dlovv (2.27) 

i 

where i labels all the neutral Higgs bosons of the extended Higgs sector, 
and 0° denotes the SM Higgs. The sum rule (2.27) ensures the unitarity 
of VlVl — > VlVl. Moreover, to ensure // — > VlVl unitarity we ought to 
demand 

9hyv9h1 fj = 9<t>ovvg<j,o fj (2.28) 

i 

in such a way that only the contribution of all Higgs bosons cancels the effect. 
In the particular case of the 2HDM, Eq. (2.27) becomes 

9hovv + 9hovv = 9I0VV (2.29) 

where (p^VV denotes the coupling of the SM Higgs with two vector bosons^. 
We can check that in terms of the mixing angles a and (3 we get 

9h°VV . ( a \ 9hOVV f a \ fr, on\ 

= sm [P — a) ; = cos (p — a) (2.30) 

9<t>ovv g<t>ovv 

^Remember that the couphng A°VV does not appear at tree level. 
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from which the constraint (2.29) is obviously accomphshed. 

A very important phenomenological consequence of (2.29) is that the 
couphngs of H'^ to VV are supressed respect to SM ones. In particular, 
if one of the scalar Higgs bosons decouples at tree level from VV then the 
Higgs coupling at tree level of the other scalar Higgs is SM-like. In that case 
we say that one of the Higgs bosons "exhausts" or "saturates" the sum rule, 
this is a natural scenario in the Minimal Supersymmctric Standard Model 
(MSSM) in which cos (/3 — a) ~ is expected, so the coupling h^VV almost 
saturates the sum rule and is SM hke, while the coupling H'^VV is highly 
supressed. 

It is worthwhile to point out that the existence of V— Higgs— Higgs cou- 
plings is a new feature of the 2HDM respect to SM. As we shall discuss in 
section (3.1), tree level unitarity constraints involving those kind of couplings 
arc also satisfied. 

As for the sum rule (2.28), it involves also the couplings of the Higgs 
bosons to fermions, which are model dependent, we shall see in next section 
that such rule is also automatically accomplished at tree level by some sets 
of Yukawa couplings. 

2.2.3 The Yukawa Lagrangian 

The Most general gauge invariant Lagrangian that couples the Higgs fields 
to fermions reads 

-£y - V^'^^L^iUf^ + Vri'^.L^iD% + ^''^,^^^^^ 

V^t^iE]^ + ^S^L^^E]^ + h.c, (2.31) 

where $1,2 represent the Higgs doublets, $1,2 = i(T2$i,2, Vij ^^^d ^ij 
diagonal 3x3 matrices and i,j denote family indices. refers to the three 
down-type weak isospin quark singlets = (rf^j, s^^, iP^) . U refers to the 
three up-type weak isospin quark singlets = ("^^/j, c^j, t^j)^ -^i? to the 
three charged leptons. Finally, Q^j^, l^^ denote the quark and lepton weak 
isospin left-handed doublets respectively. The superscript "0" indicates that 
the fields are not mass eigenstates yet. 

Prom now on, I shall restrict the discussion to the quark sector only, 
since the extension of the arguments to the lepton sector is straightforward. 
What we can see from the model described by Eq. (2.31) is that in the most 
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general case, both Higgs bosons couple (and consequently give masses) to 
the up and down sectors simultaneously. However, this general case leads 
to processes with Flavor Changing Neutral Currents (FCNC) at tree level. 
This is basically because by rotating the fermion gauge eigenstates of the 
down sector to get the mass eigenstates we are not able to diagonalize both 
coupling matrices fj^fi^^^fi simultaneously^, the situation is similar for the 
up and lepton sectors. 

Now, since processes containing FCNC are very strongly supressed ex- 
perimentally, especially due to the small Kl — Ks mass difference, these 
processes were considered dangerous. Consequently, several mechanisms to 
avoid them at the tree level were proposed, for instance one possibility is 
to consider the exchange of heavy scalar or pseudoscalar Higgs Fields [4] 
or by cancellation of large contributions with opposite sign. Perhaps the 
most "elegant" or "natural" supression mechanism is the one suggested by 
Glashow and Weinberg [30] , who implemented a discrete symmetry that au- 
tomatically forbids the couplings that generate such rare decays. Another 
mechanism was proposed by Cheng and Sher arguing that a natural value 
for the FC couplings from different families ought to be of the order of the 
geometric average of their Yukawa couplings [5], this assumption has been 
tested by Ref. [31] in the quark sector. Remarkably, absence of FCNC at 
tree level is possible even with both Higgs doublets coupled simultaneously 
to the up and down sectors; Ref. [22] studies the possibility of suppressing 
FCNC in a general 2HDM, by demanding a 5*3 permutation symmetry among 
the three fermion families, this 5*3 permutation symmetry would be broken 
spontaneously when the SSB of the gauge symmetry occurs, more details in 
section 3.1. On the other hand, Ref. [46] shows that the decoupling limit (to 
be discussed later on) is a natural scenario to get rid of tree level FCNC in 
the general 2HDM. 

Let us first discuss the mechanism proposed by Weinberg and Glashow. 
We have said that the matrices rf'^, ^^'^ cannot be diagonahzed simultane- 
ously, and the same applies to the couple of matrices 7^^'°, So it is im- 
mediate to realize that we can supress the FCNC at tree level in Lagrangian 
(2.31), if we manage to get rid of one of the pair of matrices {jf^'^ that 
couples the up sector to the Higgs doublets, and same for the down sector. 

^Of course, we still have the possibility of a fine tuning that allows us to diagonalize 
both matrices with the same transformation. Although it is not a very natural assumption, 
it leads to a model without FCNC different from the ones obtained from arguments of 
symmetry. 
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We can achieve that, by implementing the following discrete symmetry 

$1 $1 and $2 -$2 

DjR TDjR and UjR -U^r 

so by demanding invariance under this discrete symmetry we have two cases 

• When we use Djr — > —DjR we should drop out r/^^'^and t]I^'^- So, $i 
decouples from the Yukawa sector and only $2 couples and gives masses 
to the up and down sectors. This case is known as the 2HDM type I 

• When we use Djr —>■ Djr then 77^^'° and vanish and therefore 
$1 couples and gives masses to the down sector while $2 couples and 
gives masses to the up sector. In this case we call it, the 2HDM type 
II. 

In the most general framework of multi-Higgs doublet models, this su- 
pression mechanism acquires the character of a theorem [30], the theorem 
of Glashow and Weinberg says that the FCNC processes mediated by Higgs 
exchange are absent at tree level in a general multi-Higgs doublet model if 
all fermions of a given electric charge couple to no more than one Higgs dou- 
blet. The Lagrangians type I and II discussed above clearly accomplish the 
theorem. It is important to say that we can use the same type of couplings 
for both the quark and lepton sectors or on the other hand, we can treat 
them asymmetrically giving a total of four different Lagrangians^. 

Furthermore, the Yukawa Lagrangians type I and II can also be generated 
from a continuous global symmetry. The set of transformations 

$1 $1 and $2 — > e*'^<l>2 

Djr ^ e-'^DjR and Ujr ^ e-'^^UjR 

with uj = ip,0; leads to the models type I and type II, respectively. 

The discrete symmetry arises as a special case of this continuous symme- 
try when we fix (/9 = ir. Nevertheless, we should remember that the discrete 
symmetry leads to a Higgs potential which is phenomenologically different 
from the one generated by the continuous global symmetry as we saw in Sec. 

^For instance, Ref. [11] classified these Yukawa Lagrangians as models 1, 11, 111, and 
IV. 1 prefer to use the most common notation taking into account that for instance, the 
leptonic sector could be of type 1 and the quark sector of type 11. 
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particle 




particle 


J"- 


7 


1~ 






Z 


!+-,!-+ 




o+, 0- 




0++, 0~ 








0++, 0~ 








0+-, 0-+ 







Table 2.2: Assigments to the parity (P) and charge conjugation (C), quan- 
tum numbers for the Higgs and vector bosons when fermions are introduced. 
After the introduction of fermions, P and C are not conserved separately 
anymore. However, in good aproximation CP is still a good quantum num- 
ber. 

(2.2.1). Consequently, we should be careful in choosing the symmetry. For 
instance, as we saw in Sec. (2.2.1), the potential Vb in Eq. (2.8) that coin- 
cides with the potential of the Higgs Hunter's Guide and the MSSM, comes 
from the imposition of the continuous global symmetry (plus a soft breaking 
term) and not from the discrete symmetry^, while the potential in Eq. 
(2.5) comes from the discrete symmetry. 

However, the discrete (or global) symmetry is not compulsory and we 
should explore the possibility of having the whole model Eq. (2.31) including 
the FCNC, and look for the constraints that the experimental measurements 
impose on the region of parameters of the model. When we take into account 
all terms in (2.31), the two Higgs doublets couple and yield masses to both 
up-type and down-type fermions, in that case we call it the 2HDM type 111. 

From the discussion above, we conclude that the Feynman rules of the 
Yukawa Lagrangian are highly model dependent, and could be enhanced or 
supressed respect to SM in contrast to the kinetic couplings which are always 
supressed. Consequently, the phenomenology is very sensitive to the choice 
in the Yukawa sector. 

There are however, some general characteristics related to the quantum 
spectrum J^^ for the particles involved in the Yukawa interactions. When we 

''We could argue that the potential Vb i.e. the potential in the Higgs Hunter's Guide 
and the MSSM, conies from the discrete symmetry but adding a fine tuning such that A3 = 
A4. Notwithstanding, the global continuous symmetry generates this potential without the 
necessity of any fine tuning. 
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introduce fermions, C and P arc not conserved separately any more, however 
CP is still a good symmetry (but not exact). If we assume CP conservation, 
the Higgs bosons and vectors that were considered above could be thought 
as admixtures of two states with well defined C and P (appropiate linear 
combinations of them are CP eigenstates, and even we can introduce a CP 
deviation e), we show in table (2.2) the J^*^ quantum numbers for Higgs 
and vector bosons when fermions are introduced i.e. when C and P are not 
conserved separately but CP is. 

On the other hand, it is well known that for a // system, the quan- 
tum C and P numbers are P = (-1)^"^^ , C = (-1)^+^; where L, S_ 
denote the total orbital angular momentum and the total spin of the // 
system, respectively. Then, // cannot be coupled to states of the type 
,0"*"". Therefore, in the couplings // (i7°, /i°) the Higgs bosons behave 
as purely O"*""*", while in f fA^ the Higgs behaves as purely 0~+. Hence, 
hP couples with parity +1 to matter (fermions) while ^4° couples with 
parity —1 to fermions, it is because of this fact that the Higgs boson IS usu- 
ally called a pseudoscalar boson. Obviously, ^ continue being CP— even 
and A^ continues being CP— odd as long as CP invariance is demanded. 

Additionally, there are two interesting gauge frameworks to generate the 
Yukawa couplings, the [/—gauge and the P— gauge. Though calculations 
in the unitary gauge are simpler because the would be Goldstone bosons 
are removed, the usage of the P— gauge in which the would be Goldstone 
bosons appears explicitly, could be useful if we are interested in calculations 
involving longitudinal vector bosons at high energies, since the equivalence 
theorem could be applied. 

Now we shall examine the most interesting choices for the Yukawa La- 
grangian and their features. Once again, for the sake of simplicity, 1 will 
restrict the equations and discussions to the quark sector, nevertheless most 
of the results hold for the leptonic sector too, with the proper changes. 

The 2HDM type I 

In the 2HDM type I, only one Higgs doublet (say $2), couples to the fermions, 
so the Yukawa Lagrangian becomes 



When we expand in terms of mass eigenstates the Lagrangian reads 



£y (type I) 



+leptonic sector + h.c, 



(2.32) 
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£y (type I) 




2Mw 




+leptonic sector + h.c. 



(2.33) 



The interactions involving would be Goldstone bosons given above, ap- 
pears in the i?— gauge and are identical in all Yukawa Lagrangians, however 
they vanish in the unitary gauge. 

Note that since only one Higgs doublet couples to the fermions, the La- 
grangian (2.32) is SM-like. However, it does not mean that the Higgs fermion 
interactions are SM-like, in Eq. (2.33) we see that all five Higgs bosons and 
the mixing angles from the 2HDM appears. It is because when we replace 
the scalar gauge cigcnstates contained in the doublet by the corresponding 
mass eigenstates, we are taking into account the mixing of both doublets 
coming from the potential to obtain such mass eigenstates. 

On the other hand, the Yukawa couphngs in SM are given by 



in both the up and down sectors (/ denotes any quark or charged lepton). It 
is interesting to analyze the deviations of these couphngs respect to the ones 
in SM. In order to do it, wc calculate the relative couplings (the quotient 
between the Yukawa couplings of the new physics and the SM Yukawa cou- 
plings). I will follow the notation defined in [72] for the relative couplings 




2Mw 




SM 



with "i" denoting a fermion antifermion (or VV) pair 
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and "/i" represents a generic Higgs boson. For the Lagrangian type I we get 
in third family notation 





sm a 

= cos 

sm p 


(/? — a) — cot /? sin 


w- 


a) 


xf : 


sin a 

= cos 

sm p 


(/3 — a) — cot /3 sin 


{13- 


a) 


xf ■ 


cos a 
sm p 


(/3 — a) + cot /3 cos 




a) 


xf : 


cos a 

: = sm 
sm p 


{P — a) + cot /? cos 




a) 


xf : 


—i'js cot /3 


; xf ■ ilb cot p 







There are several interesting features that could be seen above. First, the 
relative couplings of a certain Higgs boson to the down-type quarks are equal 
to the ones to the up type quarks, only the couplings to the up and down 
type quarks differs from a relative sign. Second, from Eqs. (2.30) and (2.34) 
we sec that the model type I satisfies the tree level unitarity bound given 
in Eq. (2.28). In the case of minimal mixing in the CP even sector i.e. 
a = 0, is decoupled from the Yukawa sector while the couplings of 

are maximal (for a given value of tan/3). In the case of a = 7r/2, the 
behavior is the same but interchanging the roles of both CP even Higgses. 
Additionally, the Yukawa couplings for a certain Higgs can be enhanced or 
supressed depending on the values of the mixing angles a, p. 

On the other hand, from (2.34) we can realize that in model I the heaviest 
(lightest) CP even Higgs becomes totally fermiophobic (i.e. all (/i") // 
couplings vanish) if a = (7r/2). We should bear in mind however, that de- 
cays of H° (/i°) to fermions are still pos_sible via H° (/i°) W*W {Z*Z) 
2/7 or H° (h^) W*W* {Z*Z*) 2fJ. If the fermiophobic scenario were 
accomplished it would be a signature for physics beyond the SM, and the 
2HDM type I would be an interesting candidate. In addition, we see that 
if we further assume that P — a = 71/2(0) then (/i°) becomes also boso- 
phobic (respect to vector bosons, see Eq. (2.30) and ghostphobic, while 

(H^) acquires SM-like couplings to fermions and vector bosons. In the lat- 
ter case, (/i*^) always needs another scalar particle to decay into fermions 
or vector bosons. 
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The 2HDM type II 

In this case, one Higgs doublet (e.g. $i) couples to the down sector of 
fermions while the other Higgs doublet (e.g. $2) couples to the up sector. 



-£y (type II) = r?J'°g;'^$iD5« + i'^f^,L^2U% + leptonic sector + h.c. 

(2.35) 

In this case the expanded Lagrangian becomes 

-i^y (type II) = — DU'^^'^'D (cos aH"^ -smahP) 

^■^^ ' 2Mw cos f3 ^ ^ ' 

+^T^ nUMff'^'U (sinai/° + cosa/i°) 

2Mwsm/3 ^ ^ ' 

ZIVl\Y 

+ tan 

— U \(kM^'''Pr - M^^"^XPi)l DG 
v 2Mw L V / . 

+leptonic sector + h.c. (2.36) 



w 



this is the Lagrangian that is required in the MSSM. Notice that since the 
down fermion sector receives its masses on the first Higgs doublet while the 
up sector receives its masses from the second one; the hierarchy of the Yukawa 
couplings between the top and bottom quarks of the third generation becomes 
more natural if tan/? ~ 35 and tana ~ 1. On the other hand, we can see 
that the couplings of up type fermions to Higgs bosons are the same for both 
models type 1 and II. It is because in both lagrangians (2.32) and (2.35) the 
up sector couples to the same Higgs doublet while the down sector couples 
to different ones^. 



we had chosen in the model type II to couple the doublet $1 to the up sector and 

$2 to the down sector, the couplings of the down type fermions would had been the ones 
that coincide for models type I and II. It is a fact of convention. 



32 



Once again, comparison with SM is useful, calculating the relative cou- 
plings in the third generation notation, we get 



Yf : = cos (p — q; — cotp sm (p — q; 

* sm(3 / A- VA- / 

itO cos a n ■ /n \ 

Xh : — cos [jj — a) + tan jj sm [jj — a) 

COS L/ 

1,0 COS a . / r, \ r, I r, \ 

V: : = sm p — a + cot p cos p — a 

'■ — ^^^^ = sin (3 — a) — tan 3 cos (3 — a) 

" cos (3 / ^ \^ ; 

xf ■ -i75Cot/3 ; x^" : -i75tan/? (2.37) 

from this results and Eqs. (2.30), we can sec the following interesting rela- 
tions among the relative couplings in the 2HDM type II 

(x;o^(x•?T = 1 

{xf'"' - xt"'"") (xj:-"" - xf'"') = 1 - {x'S'-'y 

(x;:"'"" + xf'"") x'^-"° = 1 + xr'xf 

2 \ 2 2 

From Eqs. (2.30) and (2.37) wc see that the model type II satisfies the 
tree level unitarity bound given in Eq. (2.28). The relative couplings of a 
certain Higgs to the down-type quarks are different from the ones to the up 
type quarks unlike the case of Lagrangian type I . 

On the other hand, if cos {P — a) ~ 0, the couphngs of the lightest 
CP even Higgs boson /i° are almost identical to the SM Higgs couplings 
not only in the Yukawa potential but in the kinetic sector as well, if we ad- 
ditionally take H^, sufficiently heavy while keeping the quartic Higgs 
self couplings < (9 (1), we obtain the decoupling limit in which the low energy 
effective theory is the SM [46]. In general, if one saturates the sume rules 
for VV with one scalar, the unitarity sum rules can be achieved by assuming 
SM-likc couplings of the same Higgs with fcrmions. This is a natural scenario 
in the MSSM. However, it is not the only way to satisfy the sum rules, for 
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instance the other Higgs (that is obviously weakly coupled to VV) could be 
strongly coupled to fcrmions so that the product ghwOhfJ ^^^^^ significant. 

It is worthwhile to point out that there is a very strong difference between 
models I and II, in model II when tan/5 > 1 some couphngs (hke Att) 
are supressed, and others like Abb are enhanced, for instance in H'^tb the 
contributions of the top (bottom) are suppressed (enhanced). In contrast, 
model I gives a uniform supression or enhancement pattern (for fixed a). 

Moreover, from (2.37) we can check that model type II does not exhibit 
a totally fermiophobic limit for any Higgs boson, unlike the model type I. 
For instance, if a = n/2 then /i° (if°) become fermiophobic to fermions of 
up(down)-type, while couplings to the down(up)-type fermions are maximal, 
the opposite occurs when a = 0. This partially fermiophobic behavior of the 
model type II is another important difference with model type I and with 
SM. 

The 2HDM type III 

The model type III consists of taking into account all terms in Lagrangian 
(2.31) 

+lepton sector + h.c. (2.38) 

in the case of the model type III, we are able to make a rotation of the 
doublets in such a way that only one of the doublets acquire VEV. Therefore, 
we can assume without any loss of generality that ($1) = f/v^, ($2) = 
(see appendices B,C). After expanding the Lagrangian (2.38) we obtain 

-i^y (type III) = -^DM^^W ( cos aH"^ -sm ah"^) 

+^DC^D (sin aH^ + cos ah^) 
v2 

+JLd^d^,DA'' + -^UMi'^'U (cosai/° - sina/i") 
v2 2Mw 

+^UC^U [sin aH^ + cosah^] - -^UC^-f^UA^ 
V 2 V 2 

+U {K^^Pr - ^^KPl) DH+ 

+Goldstone interactions + leptonic sector + h.c. (2.39) 
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where the Goldstone interactions are the same as in the case of (2.33) and 
(2.36). 

There are several interesting characteristics: the diagonal couplings of the 
pseudoscalar Higgs boson and the charged Higgs bosons to fermions are not 
proportional to the fermion mass as is the case of models type I and II; 
they are proportional to ^ff i.e. the diagonal elements of the mixing matrix 
^. Additionally, the diagonal couplings involving CP— even Higgs bosons 
have one term proportional to the fermion mass and other term proportional 
to ^ff. Once again we shall calculate the relative couplings of the model 



Xt '■ cos a H 7= V 

jjo , 6b sin a 

Xb : cosa+ v 

h" ■ , itt cos a 
Xt '■ —sma-\ ^= V 

hp . , ^fccoso; 

Xb : - sm a + — = v 

y Ivrib 

Xt ■ n — ^ ' Xb ■ -1= — ^ (2.40) 

y2mt y2mb 



The presence of diagonal mixing vertices could play a crucial role in 
looking for FCNC. At this respect, it is worthwhile to point out that the 
relative couplings in model type III are not universal because of the contri- 
bution of the factor ^ff/mj^ as a manner of example the relative couplings 
X? , x^ ! Xu in general different^, a deviation from the universal be- 
havior could be a clear signature of FCNC at the tree level; even if only 
diagonal processes are observed. In particular, we cannot obtain a totally 
fermiophobic limit of (say) the Higgs hP unless that a very precise pattern of 
the quotient C///"^/ is demanded, we can have for instance the possibility of 

we can easily see from Eqs. (2.34, 2.37) in model type I, the relative couplings for 
a certain Higgs h, are the same for all fermions. In model type II, these relative couplings 
are equal for all up type fermions and for all down type fermions. In contrast, in the model 
type III all these relative couplings arc in general, different each other. This breaking of 

the universality of the relative couplings in the model type III owes to the presence of the 
flavor diagonal contributions ^u/mi, unless we assume them to be universal as well. 
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/i° to be top-phobic^'' but not charm-phobic or up-phobic, this is due to the 
lack of universahty discussed previously [38]. In contrast, model type I, can 
only be fermiophobic to all fermions simultaneously [11], while model type II 
can only be fermiophobic to all down-type or all up-type fermions simultane- 
ously. Moreover, as it was explained above, A° couples to a pair of fermions 
// through the matrix element ^//. Therefore in model III, if ^// = for 
a certain specific fcrmion, then becomes fermiophobic to it (and only 
to it). It is also interesting to see that the matrix elements modify the 
charged Higgs couplings as weU: gn+fifj = {Kikikj^R- iikKkiPh), where 
K is the Kobayashi Maskawa matrix and Pl{r) are the left(right) projec- 
tion operators. At this respect there arc two features to point out i) the 
flavor changing charged currents (FCCC) in the quark sector are modified 
by the same matrix that produces FCNC, ii) in the lepton sector FCCC are 
generated by the same matrix that generates FCNC [39] . 

Finally, we can see that the tree level unitarity constraints (2.28) impose 
a condition on the diagonal elements ^// of the mixing matrix, preventing 
them to be arbitrarily large. 

Since in model III we can find a basis in which we are able to get rid 
of one of the VEV's, the mass eigenstates become simpler. By using say 
($2) = 0, we obtain the mass eigenstates by setting /3 = in Eq. (2.9). We 
shall call it the "fundamental parametrization" . In that case, only the 
doublet $1 generates the spontaneous symmetry breaking, thus it is the only 
one that provides masses for the fermions and vector bosons. The second 
doublet $2 only provides interactions with them. 

Additionally, taking ($2) = and using the basis of states (/ii, /i2, 92-1 H^) 
we can see that $1 corresponds to the SM doublet and hi to the SM Higgs 
field since it has the same couplings and no interactions with /i2, fl'2[32]. Con- 
sequently, all the new scalar fields belong to the second doublet $2- More- 
over, h2,g2 do not have couplings to gauge bosons of the form h2{g2)ZZ or 
h2{g2)W'^W~ . Nevertheless, we ought to remember that hi, /12 are not mass 
eigenstates, they mix through the angle a (see Eq. (2.9)). However, as soon 
as we consider a scenario with a ~ these features become important since 
the set {hi, h2, g2, H"^) becomes a mass eigenstates basis. As we will see 
later on, when a = in the fundamental parametrization, we arrive to the 

^°With the top-phobic hmit for certain Higgs, we mean that the coupling of that Higgs 
to a pair of top quarks vanishes. Nevertheless, the couplings of such Higgs to a top and 
any other quark could exist. 
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decoupling limit of the 2HDM. 

On the other hand, we can parametrize the Yukawa Lagrangian type III 
by preserving both VEV's different from zero, we shall see below that the 
relation between the model type III with the models type I or type II is 
clearer with this parametrization[36]. In order to distinguish between the 
parametrization with only one VEV Eq. (2.38) from the parametrization 
with two VEV, we write down the latter in the following way 



-A 



+lepton sector + h.c. (2-41) 



U.D ^U.D 



where {^[) = vi, ($2) = ^2- The relations among and y,^^ 

are calculated in appendix (B). First we split the Lagrangian (2.41) in two 
parts 



-£ 
-£ 



Y(U) 
Y(D) 



-£y (type III) 



Vf/^,L^,D% + g'°g-i$2/^J« + lepton sector + h.c. 



-£ 



Y(U) 



- £ 



Y(D) 



In order to convert this lagrangian into mass eigenstates we make the unitary 
transformations 



Dl,r = {Vl,r)DIj, 
Ul,r = (Tl^r) U^ r 



(2.42) 



from which we obtain the mass matrices. In the parametrization with both 
VEV different from zero we get 



j^diag 



Vl 
Tl 



V2^ V2 . 



(2.43) 



We can solve for obtaining 
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^^0 ^ ^yt^^^as^^ _ !^~D,0 

^■0 ^ ^jIm^^^T^ - 



i-z) (2-44) 

V2 V2 



Let us call the Eqs. (2.44), parametrization of type I. Replacing them 
into (2.41) the expanded Lagrangians for up and down sectors are [36] 

-£vnn (type HI) = —r^ -UMfrU (sin a'H"" + cos 

^ 2Mwsmf3 ^ ^ ' 

^ VlfU [sin {ol -p)H° + cos (a' - (3) 



V2sm(3 



ig cot f3— diag^, jr AO 



I 



2M^ ^ V2sin/5 
V2Mw sin/3 



+/I.C. (2.45) 



-i:Sn^ (type ni) = — -DW^'^'D (sin + cos q;'/i°) 

) 2MvKsin/3 ^ ^ ' 

j}- DlpD [sin {ol -(3)H^ + cos (a' - (3) h^] 

v2 sin/3 

i^COt^-^ d,a, _ ^ j^j^O 

2Mw ^ V2 sin f3 ' ^ 

ocot^-^ + - ^UK^''PrDH+ 

V2Mw sin/3 ' 

^-^DM^'^'-f^DG'^ + -j^ — UKM^'''PrDG+, 



2Mw V2Mi 



+leptonic sector + h.c. (2.46) 

where K is the CKM matrix, ^^(^) = TL{VL)rj^^^^'^Tl{VR)^ and similarly 
for The superindex (/) refers to the parametrization type I. In order 

to make the expansion above, we have used the Eqs. (2.9) but replacing 
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a ^ a' since the notation a will be reserved for the mixing angle between 
CP even scalars in the fundamental parametrization (see appendix B). The 
leptonic sector is obtained from the equation (2.46) replacing the down-type 
(up-type) quarks by the charged leptons (neutrinos). 

It is easy to check that if we add (2.45) and (2.46) we obtain a Lagrangian 
consisting of the one in the 2HDM type I [17], plus some FC interactions, 
i.e., — <£y([/) — "^yp)- Therefore, we obtain the Lagrangian of type I from 
Eqs. (2.45) and (2.46) by setting = ^ = 0. 

On the other hand, from (2.43), we can also solve for 77^'°, rf^'^ instead of 



which we call parametrization of type II. Replacing them into (2.41) the 
expanded Lagrangians for up and down sectors become [36] 



to get 




(2.47) 



£'^\l. (type III) 




+/I.C. (2.48) 
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-4'2)) (type III) = ^ DM'^^^D (cos - sin 



2MvK cos /5 
1 ^ 

4 



D^^D [sin (a' -p)H^ + cos (a' - /?) h^] 

v2 cos/3 

ZMw V 2 cos p 



V2Mw ^ cos/3 



+leptonic sector + h.c. (2.49) 



The superindex (//) refers to the parametrization type II. Moreover, if we 
add the Lagrangians (2.45) and (2.49) we find the Lagrangian of the 2HDM 
type II [17] plus some FC interactions, —£^y(ij) ~~ ■^y{d)- Similarly like before, 
Lagrangian type II is obtained setting = = 0, in the total Lagrangian 
— — £Y^ij)y Therefore, Lagrangian type III can be written as the model 
type II plus FC interactions if we use the parametrization type I in the up 
sector and the parametrization type II in the down sector, it is valid since 

and are independent each other and same to rj^'^ . 

Moreover, we can build two additional parametrizations by adding —£y^^y 

^Y{D) ~'^Y{u) ~ '^y\d)- former corresponds to a model type I plus 
FC interactions but interchanging the role of the doublets $i ^ $2; the lat- 
ter gives a model type II plus FC interactions with the same interchange. On 
the other hand, terms involving would-be Goldstone bosons are the same in 
all parametrizations in the R-gauge, while in the unitary gauge they vanish 
[17, 36]. 

Finally, as it is demostrated in appendix (B) the FC matrices rf'^, ^^'^ are 
related to the FC matrices rf^'^ , ^^'^ in the fundamental parametrization Eq. 
(2.39) through the following formulae (see Eqs. (B.14)), here we shall write 
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them in terms of tan (3 




Mu,D I 2 tan /3 



V Vl + tan2/5 ^l + tan2/? 



U,D 




U,D 



V 



y/1 + tan2 /3 



1 



(2.50) 



The development of these new parametrizations of the 2HDM type III, 
facihtates the comparison of it with the models type 1 and type II. In addition, 
one of these parametrizations has been used by Haber and Gunion [46] to 
study the decoupling limit of the general two Higgs doublet model. Further, 
in Sec. (4.5) we shall obtain some limits on the free parameters of the model 
in the framework of one of these parametrizations. 
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Chapter 3 



Present Status of the general 
Two Higgs Doublet Models 

3.1 Theoretical constraints 

The most important theoretical constraints come from the p parameter and 
the unitarity hmits. In SM, p = 1 at tree level. On the other hand, in 
models beyond the standard model values around p ~ 1 are expected too. In 
arbitrary Higgs multiplet representations the value of p at tree level is given 

by[l7] 

_ _ Et^.Y, m{T, + 1) - 1-^] \VT^yf Cr^y^ 

where V^.y. = < $(Tjyj) >o defines the VEV of the neutral Higgses, Yi is 
the hyperchargc value of the multiplet number i. T-i is the SU{2)l isospin 
i.e. Tj(Tj + 1) are the eigenvalues of the generators: E^=i Finally, the 

coefficients Cr y are defined as 

\ |, if Tj is integer and Yi = 0. 

in the case of the 2HDM, we have Yi ^ Y2 ^ 1, Ct^Yi = CT2Y2 = 1, = 
T2 = 1/2. Prom which we can see that the relation p = 1 is maintained at 
tree level. 

In SM the one loop correction to the p parameter [33] is dominated by 
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the top quark loop 

Ap= — ^ 

by contrast, the contributions coming from the first family are given by 
mild logarithmic terms. This fact is related to the breaking of a global 
SU {2)y symmetry [34]. This isospin symmetry "protects" the relation p — 
1 at next to leading order. The first generation of quarks preserves in good 
approximation such "custodial symmetry". In contrast, the third generation 
of quarks breaks the symmetry strongly producing a very sizeable one loop 
contribution [34] . Constraints coming from radiative corrections to p in the 
2HDM will be discussed later on. 

Other important theoretical constraints come from unitarity. As it was 
mentioned in section (2.2.2), the unitarity bound coming from VlVl — > 
VlVl scatterings is accomplished if 

dlyv ^ dlovv (3-1) 

i 

where i labels all the neutral Higgs bosons of the 2HDM and 0° refers to the 
SM Higgs. Additionally, we saw that the requirement of unitarity for the 
// — > VlVl scattering leads to 

y^.9h?vv9h9fl = 9<ifivv94>of-f (3-2) 

i 

As we saw in Sec. (2.2.2), the condition (3.1) is accomplished for the 
general 2HDM since the couplings of any Higgs to VV are universal. In 
contrast, the second condition (3.2) is accomphshed automatically by the 
2HDM type I and II, but not by the type HI one, in the latter case this 
condition provides a constraint for the diagonal terms of the mixing matrix 
as it was already mentioned. 

On the other hand, unitarity bounds of processes involving V — Higgs — 
Higgs couplings determines the form of such couplings. For instance, unitar- 
ity constraints in H^W~ — > H^W~ and A^Z — > A^Z leads to the following 
sum rule [17]. 

9H+w-ho + 9H+W-H0 = {9h"A"z + 9hoaoz) COS^ 6*1^ = ^ (3.3) 

Moreover, unitarity of A^Z W^W^ requires a precise ratio between 
9hoaoz and 5'^o^o^[17]. Similarly, unitarity of H~^W~ —>■ W^W~ determines 
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a relation among the couplings Qn+w-h" Qh+w-h^- Combining the sum 
rule (3.3) with the unitarity bounds from A^Z W^W~ and H^W^ 
W~^W~, the following values for the V — Higgs — Higgs couplings are ob- 
tained 

Qh^A^z = ^ ^ a (/5 - «) ; 9H+w-ho = ^ COS (/? - a) 

Z COS uw z 

Qhoaoz = ^ Q COS {(5 - a) 
2 COS 9w 

They coincide with the values gotten by direct expansion of the kinetic 
Lagrangian, and tree level unitarity is also automatically satisfied by y — 

Higgs — Higgs couplings. 

On the other hand, tree level unitarity bounds in the Higgs potential can 
also be derived from the two body scattering of scalars S1S2 8^84^. Tree 
level unitarity constraints for the potential in Eq. (2.5) have been cal- 
culated [43] obtaining strong upper limits on some free parameters of the 
Higgs potential with exact Z2 invariance. For instance, Ref. [43] obtains the 
following useful bound 



rriho < \ \ — cos^ (3 - mj^o cot^ /? 



Notwithstanding, Akeroyd et. al. [44] showed that by introducing the 
soft breaking term (from which we arrive to the potential Va Eq.(2.6)) 
such bounds are weaken considerably. Additionally, Ref. [44] also shows that 
the inclusion of charged states not included in [43] changes significantly the 
bounds on mfj± and rriAO- In the absence of the soft breaking term, general 
unitarity bounds are m^± < 691 GeV, m^o < 695 GeV, m/jO < 435 GeV, 
"^^0 < 638 GeV. These bounds are obtained for relatively small values of 
tan/3 (tan/3 Ri 0.5). For large values of it the upper bound is smaller, but 
only is very sensitive to the variation in tan j3, especially for small values 
of the soft breaking term. For = 0, the sensitivity is very strong, but for 
large values of e.g. /i| = 15, the plot rriho vs tan (3 becomes a horizontal 
line predicting m/jO < 670 GeV. The relaxation of the strong correlation 
between m/jO and tan/3 could be a mechanism to distinguish between the 
potential with the exact symmetry and with the softly broken symmetry. 

Unfortunately, no general unitarity bounds have been calculated for the 
potentials V Eq. (2.4) and Vb Eq. (2.8). However in the fermiophobic limit 
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of model I, a tree level unitarity bound can be easily obtained for Vb [25] 

m^o = rriAO — 2 (A+ — Ai) f ^ cos^ (3 (3.4) 

wc sec that if wc assume /3 ~ 7r/2 then m^o ~ tri^o- As for the potential 
V, Ref. [46] estimates that in the decoupling limit the introduction of uni- 
tarity implies for the heavy Higgs bosons to be nearly degenerate. 

However, it should be pointed out that tree level unitarity is only a guide, 
its violation could imply a strongly interacting Higgs sector. Nevertheless, 
in that case perturbative arguments become suspect. 

It has also been mentioned that flavor changing neutral currents impose 
strong constraints on the 2HDM, since the introduction of the second doublet 
produces them automatically. These rare processes can be suppressed on ei- 
ther theoretical or phenomenological grounds. The most common theoretical 
mechanism consists of the imposition of a Z2 symmetry (or a global U (1) 
symmetry) as dicussed in section (2.2.3) generating the models type I and II. 
Nevertheless, there are also some theoretical considerations that permits us 
to avoid FCNC in the type HI model. Two interesting scenarios for a flavor 
conserving type HI model deserves special attention. 1) A model type HI 
with a permutation symmetry among the fermion families [22], and 2) the 
decoupling limit of the 2HDM type HI [46]. Let us discuss them briefly. 

The motivation for introducing the S3 permutation symmetry, comes from 
the similarity among the fermion families. Since there are not fermion masses 
before the SSB of the gauge symmetry, it is logical to assume that before the 
SSB the three generations of fermions are indistinguishable i.e. they hold 
an S3 symmetry. After the SSB of the gauge symmetry, fermions acquire 
masses i.e. distinguibility, thus the 5*3 symmetry should be spontaneously 
broken as well. The implementation of such permutation symmetry conduces 
automatically to the supression of FCNC at tree level [22], but also leads 
to a 3x3 unit CKM matrix. Notwithstanding, Ref. [22] showed that by 
introducing an appropiate soft breaking term, we could still get a flavor 
conserving type HI model while obtaining a reasonable CKM matrix. 

As for the decoupling limit, Ref. [46] shows that this particular scenario 
automatically suppress both FCNC and CP violating couphngs in the 2HDM 
type HI. More details about the decouphng limit will be discussed in Sec. 
(3.2.2). 

Furthermore, renormalization group equations could be utilized to obtain 
some useful constraints. For instance, we can demand that the self interacting 
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parameters in the Higgs potential remain perturbative up to a large scale 
(typically the grand unification scale i.e. Aqut ~ lO^^GeV), and also that 
the vacuum of the potential remains stable up to such large scale. In SM 
such bounds have been studied at first by Cabibbo et. al. [13] and Lindner 
[14]. The criterium of perturbativity provides an upper limit of the SM 
Higgs mass of about 180 GeV when the cut-off scale is taken to be Aqut ~ 
10^*'GeV. Meanwhile, vacuum stability provides a lower limit of about 130 
GeV [15], for the same cut-off. On the other hand, if the cut-off scale is 
extended up to the Planck scale Apianck ~ lO"*^^ GeV, the allowed interval 
becomes 145 GeV< nih < 175 GeV [16], therefore perturbativity and vacuum 
stability requirements gives a rather strong allowed interval for the Higgs 
mass in SM. As for the 2HDM, Ref. [15] has examined the constraints coming 
from perturbativity and vacuum stability up to the grand unification scale 
in the type II version without the soft breaking term in the potential. Such 
reference found that these criteria impose an upper limit for the charged 
Higgs of mH± ^ 150 GeV which is in conflict with the lower bound for 
it > 500 GeV^ Consequently, the 2HDM type II (without the soft 

breaking term) cannot be valid up to the unification scale. In contrast, 
model I is not excluded by these criteria, since the experimental lower limit 
for mH± is significantly milder. 

We should remark however that bounds from perturbativity and vacuum 
stability from the electroweak (EW) scale up to any cut-off scale is only a 
guidance, if they are not accomplished it merely means that the model is only 
valid up to a lower scale (or that it belongs to a non-perturbative regime). 
Indeed, we can reverse the problem and check until what scale might our 
model be valid. For instance, Ref. [16] have calculated lower and upper 
bounds of m/jo in the context of the 2HDM types I and II but including 
the soft breaking term in the potential. These bounds are estimated as a 
function of the cutoff scale and also of the soft breaking term of the discrete 
symmetry, requiring that the running coupling constants neither blow up nor 
fall down below A. Comparing with the corresponding SM constraints, the 
upper bound docs not change so much but the lower bound is substantially 
reduced. In particular, in the decoupling limit with a large soft breaking 
term, the upper bound for the 2HDM for Apianck = 10^^ GeV, and rrit = 175 
GeV, is about 175 GeV, which is basically the same as the corresponding SM 



^At the time of ref [15], the lower bound was significantly smaller m^f+ > 165 GeV, 
but enough for the conflict to arise. 
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upper bound. By contrast, the corresponding lower bound for the 2HDM is 
about 100 GeV while the SM bound is about 145 GeV. In the region of small 
soft breaking term the lower and upper bounds depend on the soft-breaking 
mass and when the soft breaking term vanishes, the lower bound no longer 
appears [16]. 

Finally, perturbativity could be also examined at a fixed scale (typically 
the electroweak one) to bound the free parameters of the model at that scale. 

For example, perturbative constraints at the EW scale [45] on the Aj restrict 
the allowed values of tan P, from these grounds it is found that tan /9 > 30 is 
strongly disfavoured [44]. 



3.2 Phenomenological constraints 

Let us discuss radiative corrections to p first. There are different renor- 
malization schemes, Marciano and Sirlin [48] maintain the relation myy = 
mzCOsOw and radiative corrections appear in the relation between Gp and 
mw In the Veltman scheme [49], the free parameters of the electroweak the- 
ory are g, sin 9w and mw- The effective one loop value for m| is extracted 
from i^^e and ■j^^e at one loop, and writing them in the tree level form with a 
corrected niz value. From the Veltman scheme the one loop correction Ap is 
written as 



m^Ap = Aww {k^ = 0) - cos^ OwAzz {k^ = O) (3.5) 

Avv are the propagator corrections at zero momentum transfer. 

In the i?— gauge, the terms that could be quadratically dependent on the 

Higgs boson masses come from diagrams with either two Higgs bosons or 
a Higgs boson and a Goldstone boson, as well as associated tadpole type 
diagrams. 

The contribution of the Higgs sector in the R — gauge in the 2HDM type 
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II is given by [17] 



A^%, (0) - cos' OwA^i (0) = ^ {Fap {mj,+ ,m\o) + [Fap (m^+,m^o) 

-Fap {m^, mjjo)] sin^ (/3 - a) 

-Fap {rn^Ao^mlo)] cos^ {P - a)} 

-Fap (m^, m^o) - (m|, m^o) 

-Fap (m|, m|o) + Fap m^o)] (3.6) 

where Fap {ml, mf) is defined as 

2 2 2 2 

Fa, {ml ml) ^^{ml + ml) - -^^^ In ^ (3.7) 

Zi 1 1 1-^ '"'2 '"'2 

Notwithstanding, we should be careful in removing the contribution from 
the SM Higgs, a natural approach could be to take /i° as the SM Higgs with 
modified couplings, and extract the contribution of hP with SM couplings. 

Since my/ ~ mz the HG (Higgs-Goldstone) contribution is very small. 
While the HH contribution could be large depending upon the mixing angles 
and masses. From expressions (3.5-3.7) it is easy to check that Ap depends 
basically on the splitting between mH+ and m^io as well as from the splitting 
of mH+ with the CP even Higgses. Combined experimental measurements 
on neutral currents of sin^ 6w , and mw , imply that the charged Higgs 
mass cannot be much larger than the neutral Higgs mass (except if the Higgs 
is SM like). There are some bounds about the splitting between mj^o and 
mH+- In SUSY a very natural scenario consists of m^o ~ m^o ~ mH+ with 
a SM hke hP Higgs. The term Fa, (m^+, m^o) vanishes in this case and the 
other terms tend to be cancelled by the small splitting of mw ~ ^z, then 
these contributions are expected to be small in SUSY scenarios. Notice that 
Azz (0) tends to cancel the otherwise large Aww (0) contribution. 

For the model type I, the contribution to Ap has been estimated in the 
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unitary gauge [50] 



Ap = iQT^2y2 [^^^^ {p-a)F {mH± , m^o, m^o) 

+ cos^ (,(? — a) F (mj^±, m^o, m^o)] 

be b ab a ac a 

F (a, b,c) = a + m m m - 

b — c c a — b b a — c c 

In the case of model type I, if a light /i° Higgs still exists, the variation 
of Ap with is soft, but it depends on m^o and m^o — mH± (or mH± and 
mAO—TnH±). Once again, Ap can be used to restrict the splitting between the 
charged Higgs and neutral Higgs bosons. Ref [51] shows that in the general 
2HDM when m^o is small, the contribution to Ap can grow quadratically 
with the other Higgs masses, so in order to keep Ap small, a correlation 
among the Higgs boson masses to cancel this large loop contributions should 
be assumed. Ref [52] presents constraints in the framework of the general 
2HDM with a very light CP-odd scalar, these bounds are shown in a rnn — 
mH+ plane for three different values of sin^ (/3 — a) , constraining strongly 
the splitting between such scalars. 

Other observables useful to constrain the 2HDM are Br {b sj) , the Z — > 
bb hadronic decay branchig ratio {Rb) , the forward backward asymmetry of 
the bottom quark in Z decays (Aft); and the muon anomalous magnetic mo- 
ment of the muon, a^. The effectiveness of each observable to constrain the 
model depend on the specific scenario, for instance, when tan /3 is of the order 
of 1 the asymmetry is less effective to constrain in the m//+ — tan f3 plane 
than Rh [52]. The anomalous magnetic moment of the muon is a very 
important parameter to restrict any physics beyond the SM, in particular 
the 2HDM. Since much of our analysis is based on this observable, we shall 
revise its present status in section 4.1. For now, I just say that its present 
world average experimental value is given by [68] 

/ r)\cxp 

^exp ^ y9 - )^ ^ (8) X 10"^° 

and the SM estimation has been calculated recently by a variety of authors 
[70], though their estimations are rather different, all of them coincide in the 
fact that the SM value based on e^e~ — ^hadrons data, is roughly Scr below 
the present world average measurement (see details in section 4.1). Thus, the 
theoretical and experimental current estimations seem to show the necessity 
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of new physics. The fact that Aa^ should be positive at 95% CL, impose 
strong restrictions to any physics beyond the SM. 

On the other hand, bounds on the 2HDM might also be extracted from 
the non-observation of direct Higgs production and/or decays. For instance, 
in scenarios with a very light CP odd Higgs boson constraints can be derived 
from Br(T-> Z AAA, f {f')J Z (W^) AA -> Z (W^) + 

"77" and // ^ H+H~ Z {W+W-) + "77" [52]. Similar decays with 
yl° replaced by can be used to bound scenarios with a very light CP even 
Higgs boson. For intermediate Higgs bosons, limits could be derived from 
the non observations of processes like e+e" Z* ^ ZhP, and in the heavy 
regime from processes like — > ZZ — >^leptons. 

Further, as well as constraining the model it is necessary to explore the 
possibility to produce the Higgs bosons from colliders. The expected observed 
production rate is usually very model dependent and, of course depend also 
on the experimental luminosity and detector efficiency. 

At e^e~ machines production mechanisms depend upon substantial VV 
couplings. The most advantageous scenario for these machines would be a 
roughly equally VV couplings for both scalars, if one of the Higgses tend 
to saturate these couphngs (as in MSSM) then the discovery of the weak 
coupled Higgs would be problematical in e+e" machines. ^4° is particularly 
problematical in these colliders since it is not coupled to VV, the main pro- 
duction mode that is available is e^e~ — > Z* — > A^h^ or A^H^ . The cross 
sections of these processes are proportional to cos^ (/5 — a) , sin^ (/5 — a) for 
hP, respectively. Detection might be possible if m^o -|- m/jo (or m^o -|- 
mHo) is not too large compared to the machine energy, and the vertex 
ZA^{h^)H^ is a saturating coupling. Further, if hP,A^ arc sufficiently 
light, they may be produced with a Z on shell and the production rate 
would be higher. With sufficiently high energy, the heavy states can be 
produced via e+e^ H^H^ or e+e~ — > We should notice that 

modes in which one or both Higgses go to ee are difficult to detect be- 
cause of the QED background. Prospects for detection of a single Higgs in 
e~^e~ colliders are e+e" Z Zh^ (A°) where one of the Z's is off shell. 
If m^o > 2mz "the gold-plated" detection mode is ZZ ^Icptons (for 

rriho > 130 GeV, /i° ZZ* could be a signature). In the intermediate mass 
Higgs regimes other decays are required [66] . With increasing center of mass 
energy, W~^W~ fusion begins to be the dominant mechanism. 

In hadron colliders, the dominant mechanism for Higgs production is the 
gg fusion through a top quark loop. Important prospects for single Higgs 
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production are gg A'^, H'^ and gb ^ H t. An interesting method could 
be from gg ^ QQ , , H"^) where Q is a heavy quark (6 or t) and the 
Higgs decays into a Heavy quark pair. Another interesting source is the 
quarkonium decay to which could be enhanced or supressed according 
to the model. Indirect signatures could be given by ^Si (tt) — > bb, the SM 
contributions come from W exchange, but in two doublets a could be 
placed wherever the is. So alterations of the rate could be detected 
especially if tan (3 < 1 so that the tbH~^ coupling is enhanced. Many other 
quarkonium decays can be a direct source for Higgs boson pairs. Some of 
them take into account exotic quarks with new CKM elements. 

Another possibility is the 77 collider mode via 77 fusion (dominated by 
pW" or tt loop), depending on the mass of the Higgs in the final state decay. 

On the other hand, muon muon colliders offer interesting perspectives 
for Higgs discovery, because Higgs bosons s— channels are enhanced respect 
to e+e~ colliders, though still weak (proportional to the muon mass) but 
production can be significant if the collider is run on the Higgs resonance 
^/s = rriH- Further, they could be made feasible to reach energy in the 
multi-TeV regime. Many references have address the perpective of looking 
for the Higgs bosons in this kind of machines [21]. In particular, the fourth 
of Refs. [21] shows that in the large tan/3 regime, production of a single 
charged or pseudoscalar Higgs in association with a gauge boson is possible 
in muon colliders (i.e. H^W^,A^Z) with sizeable cross sections 

whose anologies at e"'"e~ would be very small. 

As for e"'"e~ machines, LEP have obtained limits on the charged Higgs 
mass from e+e~ — > 7*, Z* — > H^H~ . It worths to point out that far from 
the Z mass shell electromagnetic interaction become dominant, thus the 
bound on mH+ obtained from this process is nearly model independent, LEP 
excludes a region of 2mT- < m//+ < SO.SGeK The region 711^+ < 2mr is 
excluded by the non observation of i? — > H^Xc- In this region they look for 
H'^ T-dr when mH+ > rur- 

Some other limits arc gotten from some rare decays like b — > S7, b — >■ 
sg, K — > Tidd. In particular, for the model type II a lower bound of ^ 
500 GeV at 95% CL is reported based on NLO analysis of the b s^j data 
[67]. Further, as it was mentioned above constraints on the niH —TnH+ plane 
are obtained based on Ap and on the mH+ — tan (5 plane from R}, and A}, [51, 
52]. 

The top and charged Higgs phenomenologies are intimately related. Since 
the present upper bound of mH+ > 500GeV for model II seems to show 



52 



that mH+ > rrit + rrib we shall examine this assumption first. In that case 
two important scenarios emerge. In the first, mH+ < mw + iriho so that 
i7+ — > W^h^ is forbidden, thus — > tb dominates fermion modes. But if 
iJ+ — > W^hP is allowed there is a competition. We get 

BR{H+ W+h^) _ 2cos2 {(5 - a)p^m^+ 
BR [H+ tb) ~ %M 

M = [{ml cot^ (3 + ml tan^ P) (m^+ - m] - ml) - Amlml] 

for we replace cos {(3 — a) sin (/3 — a) . So for large mH+, — 
W^h^ could be important or even dominant for cos {j3 — a) not too small, 
because of the availability of longitudinal W^'s, if we assume saturation of 
/i° i.e. cos^ (/5 — a) = 1, we see that for large mH+ the mechanism —>■ 
W'^h^ would be dominant. However for masses accesible to the LEP II 
iJ+ — > W^h^ is smaller than i7+ — > tb. Additionally, since the condition 
mt < + mi), is obviously accomplished and taking into account that 

t W'^b is the dominant top decay, the off shell contribution of the charged 
Higgs is very small even for light Higgses; in that case it could be interesting 
to examine branchings of the type BR{t — > H*'^b — > t+'j?^^). On the other 
hand, if mH+ < mt + m^ and mH+ < mw + m^o (as it is still possible in 
models I and III) only — > T^dr and cs are relevant. So we 

examine the ratio of their branchings, in model I this ratio is independent of 
tan /3, while in model type II it yields 

BR {H+ T^d^) _ pr ml tan^ /?cos2 Oc (m|+ - ml) 
BR {H+ cs) ^ 3pc K 

K = (ml tan^ P + ml cot^ P) + — ml — ml) — Amlml tan /3 cot P 

in model I it is ~ 30% and independent of tan /3. We should remember that 
small values of tan/5 with small values of mH+ are probably excluded by 
B^ — B^ mixing. Now, if t — > H~^b is allowed, it is very likely that 
mass is sufficiently large to produce W^h^. Since iJ+ tb is not al- 

lowed then W^h^ will be dominant because of the H'^W^h^ coupling 

{{g/2) cos {P — a)), as long as it is not very supressed. So a possible sce- 
nario for H~ discovery could be tt production followed by the conventional 
t — > W'^b and t — >• H~b — > W^h%. Since the is most likely to decay 
to bb, we would have to distinguish between a W~^bbb final state from a 
W~^b final state. 
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Furthermore, the loop induced — > W^Z decay, though highly sup- 
press in most of the parameter space, can be significantly enhanced owing 
to non- decoupling efi^ects of Heavy Higgs bosons on the H^W^Z vertices. 
As we explained in Sec. (2.2.2), these vertices are forbidden at tree level in 
multi-Higgs doublet models and in particular in the 2HDM because of the 
isospin symmetry, but might arise at loop levels due to the breaking of the 
isospin symmetry through the loop particles. The quark loop contributions 
for the H^W^Z couplings have been studied in [18]. The full one loop cal- 
culation in the 2HDM was carried out in [19], finding that the inclusion of 
non-decoupling heavy Higgs modes with large mass splitting between and 

could yield a substantial enhancement for the decay width W^Z 
due to the strong breakdown of the custodial SU (2)y invariance [34] in the 
Higgs sector. Ref. [19], found that the branching ratio for this decay could 
reach values up to 10~^ ~ 10~^ for ■mH+ = 300 GeV, giving a detectable 
mode at LHC or future e'^e~ colliders. On the other hand, Diaz-Cruz et. al. 
[20] studied the same decay with an effective Lagrangian approach, finding 
a branching ratio enhancement up to about 10~^. Additionally, Ref. [20] 
also studies the rare decays — > W'^j, and — > W~^h^ with an effective 
Lagrangian approach, finding that they can have ratios of order 10~^ and 
O (1) respectively. Moreover, the relative behavior among the three decays 
if + — > (7, Z, hP) could give clues about the underlying Higgs structure 
[20]. 

Since detection of a charged Higgs would be an unambiguous signal of 
physics beyond the SM, many estimations of the potential of either hadron 
or electron positron colliders to discover such Higgs boson and its properties 
have been carried out. The production of charged Higgs at the Tevatron and 
the LHC proceeds via the partonic processes: gh — > tH^ and gg thH^ [53]. 
Other partonic processes could generate pair production of charged Higgs 
bosons: gg{qq) — >■ H^H^[54:] and gg{bb) H'^W^ [55], notwithstanding 
H'^W^ production lead to a much smaller rate at Hadron colliders. At 
e'^e~ colliders, the main source is e~^e~ — > 'y*,Z* H^H~. If the collider 
energy is not enough for pair production the associated H^W^ production 
becomes dominant if kinematically accesible. This latter possibility has been 
studied in Ref. [56], where H^W^ production is analized in both hadron 
and e^e' colliders in the framework of the MSSM and the general 2HDM. At 
hadron colliders H^W^ final states arise via gg fusion and bb annihilation, 
Ref. [56] found that hadronic cross section in the general 2HDM type II, 
can be much larger than in the MSSM (about a factor 500-1000!), thus, its 
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discovery at hadron colliders could in addition provide a distinctive signature 
of the Higgs sector it belongs to. In contrast, for e^e~ colliders the MSSM 
signals are considerably higher than for the 2HDM type II (up to two orders 
of magnitude). 

On the other hand, Kanemura et. al. [57], showed that there are only two 
channels that offers an opportunity for single H"^ detection at e^e~ colliders 
for 'mij± > when ^/s = 500 GeV, (1) the one loop H^W^ production 

discussed above, and (2) the tree level scattering e+e^ —>■ t~VtH^ ,t^i/tH~ . 
The latter is relevant in the large tan (5 region while the former is relevant in 
the small tan/3 regime. Since the rate of process (1) is rather poor [57], Ref. 
[58] resort to the dominant decay channel H'^ tb. Such process has as the 
main background source, the SM process e~^e~ ^ ti ^ T'Vrtb. It was found 
that for 1 and 5 ab^^ of accumulated luminosity, neither evidence (> 3cr) 
nor discovery (> 5a) of charged Higgs bosons is possible if mH± ^ \fs/'^-i 
whereas if mH± ^ V^/^ the signal should be easily observed. 

Further, many of the technics to look for the SM Higgs boson, can be 
extrapolated to the exploration for neutral Higgs bosons of extended Higgs 
sectors, if a neutral CP even scalar were discovered a precision test of its 
couplings to fermions and vector bosons, as well as its self couplings, could 
give us information on the underlying Higgs structure. Linear colliders (LC) 
will explore the Higgs structure in case of discovery by precise measurements 
of the Higgs couplings to fermions and vector bosons. However, the structure 
of the Higgs potential can only be revealed by measuring the Higgs self 
couplings. The trilinear Higgs coupling XhOhOho can be directly measured 
from e+e- ^ Z* ^ ZhPhP and e+e" ^ W+*W-*vV hPhPvV, when the 
Higgs boson is light [59]. Ref. [60] shows that at e^er colliders with the 
energy of 500 GeV (3TeV) and the integrated luminosity of 1 ab^^ (5 a6^^), 
XhOhOh'-'' can be measured by about 20% (7%) accuracy for rrih ~ 120 GeV. It 
worths to remark that even in the decoupling limit in which the couplings of 
hP to fermions and vector bosons are SM-like, the one loop corrections to the 
trilinear self coupling Xh'-nfiyfi can undergo a significant deviation from the 
SM behavior in both the MSSM and the general 2HDM [61]. In particular, 
they find that quantum corrections can be of the order of 100% in the general 
2HDM due to the quartic mass terms of heavier Higgses. In the case of the 
MSSM they find that the deviation can exceed 5% in light stop scenarios. 

The rr decay mode [62] is a promising discovery channel for A° and 
at the CERN LHC for large tan/?. On the other hand, the potential 
of the muon pair channel for MSSM Higgs bosons was studied in [62, 63]. 
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Ref. [64] focus on the potentiality to look for neutral Higgs bosons of the 
MSSM via pp bbH bbfiJI + X at the LHC. He estimated that such 
process could serve to discover and at the LHC with an integrated 
luminosity of 30 fb~^ if m^o < 300 GeV. At a higher luminosity of 300 fb~^, 
the discovery reach is not expanded much. In addition, the inclusive final 
state H lip could allow the discovery of and at the LHC with an 
integrated luminosity of 30 fb^^ if rrij^ < 450 GeV. At a higher luminosity of 
300 fb~^, the discovery region in m^o increases significantly to < 650 GeV 
for tan /3 = 50. 

In this section we have examined the most hopeful channels for possible 
direct and indirect detection of the 2HDM Higgs sector. However, the im- 
portance of this channels depend strongly on the specific scenario, in what 
follows we describe some interesting scenarios of the general 2HDM. 

3.2.1 The fermiophobic limit 

If we take into account the relative Yukawa couphngs given in Eqs. (2.34) 
we can see that in model I the heaviest (hghtest) CP even Higgs becomes to- 
tally fermiophobic (i.e. all hff couplings vanish) if o; = (vr/2). In contrast, 
model type II does not exhibit a totally fermiophobic limit for any Higgs 
boson, see Eqs. (2.37); for instance, if a = 7r/2 then /i° (i?°) become fermio- 
phobic to fermions of up(down)-type, while couplings to the down(up)-type 
fermions are maximal, the opposite occurs when a — 0. On the other hand, 
in model HI, there is still a possible fermiophobic scenario, nevertheless in 
order for a certain Higgs boson to be totally fermiophobic, a non natural fine 
tuning must be accomplished^; see details in section (2.2.3). 

Consequently, only the model type I provides a natural totally fermio- 
phobic framework. The fermiophobic limit for (a — 7r/2) in the model I 
can be obtained in the potential in two ways: A+ = 0, or = ; the 
latter assumption leads to a massless h^. In potential Vb there is only one 
way: 2viV2X+ = We have assumed Vi < V2 but of course we can in- 
terchange the role of each VEV. Despite /i° becomes totally fermiophobic 
by setting a — 7r/2; hP can still decay into fermions by the channels /i° — > 
W*W {Z*Z) 2ff or by means of decays of the type /i° ff through 

^There are other scenarios with fermiophobic Higgs bosons such as the Higgs triplet 
model of Rcfs. [23], and some models with extradimensions [24]. Ref. [24] shows that the 
latter can exhibit a more natural fermiophobic limit with less fine tuning. 
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scalar and vector boson loops; the dominant fermionic decay in the fermio- 
phobic limit is bb. 

If we additionally assume a — P = 6 = then becomes bosophobic and 
ghostphobic as well. Besides, in this scenario are also fermiophobic. 

In that case, /i° needs another scalar particle to be able to decay. Moreover, 
in the 2IIDM type I with these assumptions, H'^ acquires the couplings of 
the SM Higgs. In contrast, the signature of is model dependent. 

Briicher and Santos [25] shows that in the fermiophobic limit the most 
important three level decays for h° are 

/i° ^ WW, ZZ, ZA, WH, AA, H+H' 

and the most important one loop decays are 

— >• 77, Zj, bb 

and decays into fermions via virtual vector bosons. The dominance of any 
of them depend on the parameters of the model, especially on S = a — 
(5 — 71/2 — 13 see details in Ref. [25]. The fact to emphasize is that the 
branching fractions depend strongly on the potential chosen and that hP — > 
U> is the fermionic dominant decay. In addition, the decay hP — > 77 is 
always important for light hP (and sometimes for intermediate hP), but for 
hP sufficiently heavy it is suppressed. 

As for the pseudoscalar, below the ZhP and the W^H^ thresholds A^ 
decays mainly into fermions, but if 5 ^ 0, ^4^^ becomes stable or decay 
outside the detector; thus the only signature would be missing energy or 
momentum. Over the threshold of Zh^ or W^H"^ the pseudoscalar Higgs 
boson decays inside the detector. Finally, the signature of A^ is the same for 
both and Vb- 

Now for the charged Higgs boson, if mj^± is below mw± the decays are 
fermionic and independent of 5, but if 6 decreases the branching is unchanged 
while the total decay width decreases (so becomes stable), if mH± > 
mi4/± then W^'~f becomes dominant for tiny S but for large values the 

signature is again fermionic, but as soon as H"^ W^A'^^h^) is possible the 
sum of both becomes dominant, except in the case in which m^f± > mt+nib in 
whose case if+ — > tb will be dominant. In the case of potential Vb below 
the W^j and above the W^h^ or PF^A^ thresholds, the situation is the 
same as in the potential Va- But in the regime of H"^ W^j dominance, 
the situation depends strongly on the choice of parameters; for instance. 
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when 5 — > the unitarity constraint for Vb given by Eq. (3.4) predicts 
a degeneracy between m^o and mho; such degeneracy can suppress — >■ 
Vl^^7 to a few percent respect to the fermionic decays even above the mj^± > 
mw± threshold unhke the case of V^. 

Finally, is not fermiophobic in this scenario so the one loop decays 
are not important, the dominant decay is obviously bb and — > 

WW above the WW threshold like in SM, this situation is the same in both 
potentials. A difference in signature can only be detected by purely scalar 
decay modes. 

From the experimental point of view, LEP and Tevatron have looked for 
fermiophobic Higgs bosons. The main channel at LEP has been e+e~ hPZ, 
hP 77 obtaining a lower bound of about 100 GeV [26], the channels 
e^e~ — * h^A^^ 77 and WW* have also been considered [26]. As 

for Tevatron run 1, the mechanism used is qq' ^ V* ^ h^V, — * 77 and 
a lower limit of about 80 GeV was obtained [27]. Nevertheless, such limits 
were obtained assuming that the couplings of the type QhPvv a-re of the same 
order as the ones in SM. Notwithstanding, in the limit of fermiophobia we 
can see that the relative couplings Xv behave hke 

therefore, in scenarios with large values of tan j3 the couplings QhPvv might 
be highly suppress and the bounds described above would not be valid. If 
it is the case, a light Higgs boson {h^ « 100 GeV) could have eluded the 
LEP and Tevatron constraints and may also escape detection at Tevatron 
run II [28]. Owing to it, Akeroyd and M. A. Diaz [28], have proposed new 
production mechanisms at Tevatron run II, they are: (1) — > 7*,^"* — > 
H+H-- (2) qq' ^ W* ^ H^hP] (3) qq' ^ W* ^ H^A^; (4) qq ^ Z* ^ 
A^'h''. 

The subsequent decays h^W*, A° ^ Z*h^, ^ '-ff would give 

rise to 7777 final states very easy to measure accurately. In the fermio- 
phobic scenarios, these new channels are effective even when h^VV are very 
suppressed [28]. These mechanisms are complementary to the traditional one 
qq' — > Vh^ which could suffer a strong suppression in realistic models. Note 
that gluon gluon fussion is not relevant in this framework and that the four 
additional channels proposed are based on the couplings Qhhv- 
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3.2.2 The decoupling limit 

The discussion of the decouphng hmit will be given based on the following 
parametrization of the most general CP invariant potential[46] 



V = mii$l$i + 77122 *2*2 - 771i2*l*2 + h.C. + ^Ai ^$1*1^ 

+ \\2 ($^$2)' + A3 ($^$2) + A4 ($I<f2) ($^$1) 

^As (^$1$2) V Ae + A7 (^\^-^ $1$2 + h.c 



+ 



In principle 771^2! ^5; ^6 and A7 can be complex, but for simplicity all of 
them are chosen real. The exploration of the decoupling limit is inspired in 
the idea that the standard model could be an effective theory embedded in 
a more fundamental structure characterized by a scale A which is supposed 
to be much larger than the Electro Weak Symmetry Breaking (EWSB) scale, 
V = 246GeV. In particular, many models end up in a low energy model with 
a non-minimal Higgs sector, as is the case of the MSSM whose low energy 
Higgs sector correspond to that of a 2HDM. 

Based on this idea, there are two important scenarios to differenciate 
in the 2HDM, in the first one there is not a low energy theory containing 
only one light Higgs boson i.e. no decoupling limit exists. In the second 
case, the lightest CP even Higgs boson is much lighter than the other Higgs 
bosons whose masses are of the order of a new scale h.2HDM- In particular, 
if A2HDM >> V and all dimensionless Higgs self-coupling parameters ac- 
complish the condition Aj < O (1) , then all couplings of the light Higgs are 
SM-like, the deviations lie on the order of O {v"^ / A\hdm) •• t'^i^ called the 
decoupling limit. What we do in an effective theory is to integrate out the 
heavy modes; if we believe that the SM is an effective theory of a 2HDM then 
one of the doublets should be integrated out. Since only one Higgs boson 
must exist in the low energy theory then the mass of the light Higgs should 
be of the order of EWSB scale while the heavy scalars masses should be of 
the order of A2hdm- And as the low energy theory only contain one doublet, 
then h^^ should be indistinguishable from the SM Higgs. For some choices 
of the scalar potential, no decoupling limit exists. Haber and Gunion [46] 
has established that no decoupling limit exists if and only if \q = = in 
the basis where m^2 = 0- Thus, the absence of a decoupling limit implies a 
discrete symmetry for the scalar potential (this symmetry could be hidden 
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in other bases). The low energy effective theory reduces the potential to the 
potential of the SM with all the SM constraints. 

When we assume that m/jO << m^o , m^io , m//± one can simply impose 
tree level unitarity constraints to the low energy effective scalar theory i.e. to 
the SM in this case. However, at one loop, the heavier scalars can contribute 
by virtual exchanges so the restrictions on the self couplings involves all Higgs 
bosons, to maintain unitarity and perturbativity the splitting of the squared 
masses among the heavier states should be of order O (v^) so they should be 
rather degenerate [46]. 

In MSSM the decouphng limit is achieved by setting m^o » mz- However, 
one loop effects mediated by loops of SUSY particles can generate a devi- 
ation from the SM expectations (e.g. light squarks contributions to /i° — > 
111 99)1 so to avoid significant SUSY contributions we should assume heavy 
SUSY particles (say of order ITeV). If the latter condition is fulfilled, is 
SM-hke even at one loop level. The leading one loop radiative correction to 
9h%h ^s order O (m|tan/9/ (?7^^o)) and formally decouples in the regime 
in which m^o >> m|tan/9 which is more demanding in the case of large 
tan/3 regime. This behavior is called delay decoupling in Ref. [47], this 
phenomenon can also occur in the general 2HDM, with tree level couplings. 

In general, the amplitudes of all loop induced processes which involve 
hP and SM particles as external states approaches the SM values whenever 
■^j ~ O iX) m^o 00, this is guaranteed by the applequist-Carazone 
decoupling theorem. 

After describing the framework of the decoupling limit, phenomenological 
consequences must be discussed. For instance, in the most general Yukawa 
Lagrangian (model type III) FCNC are generated through the mixing matri- 
ces ^l^'^^and 'rfij^ ■ These matrices are in general complex and non-diagonal; 
thus, they could be also a source for CP violating couplings between the 
neutral Higgs bosons and fermions. Notwithstanding, in the decoupling 
limit (which correspond to the condition cos (/3 — a) = in the Yukawa 
Lagrangian), both non-diagonal and CP violating couplings of vanish 
(but not for A^) and the couplings to fermions reduces to the SM 
ones. However, FCNC and CP violating processes mediated by or are 
supressed by their square masses (relative to w), due to the propagator su- 
pression, and since cos (/3 — a) ~ O (v^/m^) the suppression factor is roughly 
the same for hP, H^, A^ . Therefore, the decoupling limit is a natural mecha- 
nism to suppress Higgs mediated FCNC and to suppress Higgs mediated CP 
violating couplings in the most general 2HDM. Additionally, all vertices con- 
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taining at least one vector boson and exactly one of the heavy Higgs states 
are proportional to cos {(3 — a) and hence vanish in the decoupling limit. 

Nevertheless, Ref. [46] shows that there are some alternative scenarios in 
which cos {13 — a) = 0, but all Higgs bosons masses are of O {v) , although 
hP is SM-like, there is not an effective theory with only one light Higgs. 
Radiative corrections introduce in this case significant deviation respect to 
the trcc-lcvcl behavior. 

The scenario in which sin (/3 — a) = 0, is usually called the non-decoupling 
regime, since hP couplings deviate maximally from SM, while ones are 
SM-hke. 

Now the following question arise, if only one CP even Higgs boson were 
discovered, how could wc distinguish whether this scalar belongs to the min- 
imal Higgs sector of the SM or to a non minimal Higgs sector?. Far from 
the decoupling limit, the deviation from the coupling of the Higgs boson to 
fermions and gauge bosons (respect to the ones in SM) could provide the 
answer. However, in the decoupling regime couplings of this Higgs boson are 
similar to the SM ones and the discrimination becomes problematic; in that 
case there are two alternatives: (1) Direct production of heavy states and 
(2) High precision measurements of the couplings of the light Higgs to look 
even for tiny deviations. 

As for the first alternative, at the LHC gg and gb H~t signa- 

tures could be suppressed in the decoupling scenario. In addition, the branch- 
ing — i> ZZ which is the gold plated in SM is almost absent in this regime. 
An interesting method could be from gg — > QQ (if °, A^, H^) where Q is a 
heavy quark {h or t) and the Higgs decays into a Heavy quark pair because 
these couplings are not suppressed in the decoupling limit. In e+e^ colliders 
with sufficiently high energy the heavy states can be produced via e^e~ — > 
H^H~ or e^e~ — > H^A^ without a rate suppression, the mechanism e^e~ 
Z h^A^ discussed in section (3.2) is supressed by cos (/3 — a). These pro- 
duction mechanisms requires very high ^/s since they involve a pair of heavy 
states. Another possibility is to consider the production of one heavy state in 
association with light states (the light Higgs and the top quark are considered 
light states), however this option seems to be hopeless [46]. 

On the other hand, for the second alternative, one of the most promising 
channels seems to be Br{hP — > hh). Additionally, 77 coUiders could play an 
important role by either measuring the — > 77 coupling with sufficient 
precision or directly producing A^ and/or in 77 fusion, only the latter is 
viable in the decoupling regime. Finally, as it was pointed out by Kanemura 
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et. al. [61], another promising coupling to distinghish the underlying Higgs 
structure could be the trilinear self coupling XhOhOho which has been estimated 
to be very sensitive to the introduction of new physics in both the 2HDM 
type II and the MSSM, see discussion in Sec. (3.2). 



3.2.3 Scenarios with a light Higgs boson 

As we saw above, SM estimation of the muon anomalous magnetic mo- 
ment differs from the experimental measurement by about 3a deviations. 
In the SM the Higgs contribution to a^u is neghgible since it is proportional 
to the factor m^/m^, and the present SM Higgs mass limit provided by LEP 
is iTifi ^ 114.4 GcV. If we believe that an extended Higgs sector provides 
the necessary enhancement to get a value of the {g — 2)^ compatible with 
experimental measurement, we should examine possible factors to generate 
such enhancement. In a general 2HDM with no FCNC, it might occur due 
to two factors. First, an increment on the h^fi~^fi~ coupling proportional to 
tan/9, and second owing to the possibility of a light Higgs boson arising from 
the suppression of the vertex h^ZZ from which the LEP bound could be 
avoided. The total enhancement at one loop level would be of the order of 
[71] 

tan'f3ln(^] (3.8) 



The couplings of h^ZZ and for the 2HDM read 

gmz sin {(3 — a) gcos{P — a) 

COS Uw 2 cos &w 

so we can evade the LEP constraint by assuming sin (/5 — a) ~ 0, i.e. a non 
decoupling scenario, in that case, a light CP even scalar can be considered. 
Exploration of light scalar and pseudoscalar modes has been studied by a 
variety of authors [24, 52, 71, 72]. In general, one light Higgs boson is neces- 
sary in the 2HDM type 1 and 11, in order to get a sufficient value of Aa^ to 
be compatible with present experimental measurement. M. Krawczyk [72] 
explores the possibility of having a light scalar or pseudoscalar Higgs boson 
in the framework of the 2HDM (II), considering only the case in which the 
soft breaking term yu| vanishes. The LEP data allows the existence of a light 
neutral Higgs in the 2HDM, but not two of them could be light simultane- 
ously [73, 74, 75, 76], the constraints are given in a (m/iO,m^o) plane and 
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they indicate roughly that rrifio + m^o > 90 GeV. On the other hand, the 
Bjorken process gives an upper hmit on the relative coupling Xv ['''2], for 
instance, it is found that at 95% CL, << 1 for 171^0 < 50 GeV. Moreover, 
the Yukawa couplings Xd constrained by low energy data coming from 

T ^ /i° 7 [77, 78] while LEP do it for masses > 4 GeV. Additionally, 
the analysis of the decay Z — > /i° 7 at LEP [73] gives both the upper 
and lower limits for \xd\ ■ 

Assuming that the light Higgs or A"), is the only one contributing at 
one loop to the muon anomalous magnetic moment, the one loop diagram 
gives a positive (negative) contribution for /i° (A") . However, calculation of 
the two loops diagram with a charged Higgs in one of the loops produce 
a contribution that is positive (negative) for (hP) . Taking into account 
both contributions Ref. [72] found that the scalar Higgs gives a positive 
(negative) contribution for m^o below (above) 5 GeV, while the pseudoscalar 
gives a positive (negative) contribution for niho above (below) 3 GeV. 

Constraints from — [70], current upper limits of LEP from the 
Yukawa processes [76], lower hmits from Z — > (A°) 7 [73], upper 90% 
CL limits from T — > (A°) 7 and upper limits from the Tevatron [79] are 
combined to obtain exclusion regions in rriho — tan f3 and m^o — tan f3 planes. 
The two loop analysis of [72] shows that light scalars are excluded at 95% 
CL, but pseudoscalar masses between ~ 25 GeV and 70 GeV with 25 < 
tan/3 < 115 are still allowed. 

An alternative study of very light pseudoscalar scenarios (of the order of 
0.2 GeV) has been carried out in Ref. [52], they show that in the 2HDM type 
I and II this very light CP-odd scalar can be compatible with the p parameter, 
Br(6 57) , Ah, {g — 2)^ Br (T A^^) , and direct searches via the 
Yukawa process at LEP. For a mass of around m^o ~ 0.2 GeV, they got that 
tan/3 ~ 1 and that behaves as a CP odd fermiophobic scalar, decaying 
predominantly into a 77 pair. The second of Refs. [52] found that a very light 
CP odd Scalar can be cither allowed or excluded depending on the statistic 
used for the Hadronic contributions for a^. Therefore, a better stabilization 
of the Hadronic contributions for a'^^ is necessary to elucidate this point. 

Recently, Ref. [24] showed that additional constraints on light pseu- 
doscalars belonging to a 2HDM are obtained from K and B decays. They 
showed that masses between 100 MeV and 200 MeV might evade such strin- 
gent constraints owing to a possible cancellation in the decay amplitude. In 
that case, the coupling of ^4° to quarks is strong and may produce sufficient 
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pscudoscalar states via photoproduction. This kind of test could be per- 
formed by the Jefferson laboratory and Ref. [24] shows the perspective for 
detection of this modes in this laboratory. 

3.2.4 The 2HDM with FCNC 
FCNC, general framework 

Flavor Changing Neutral Currents (FCNC) are processes highly suppressed 
by some underlying principle still unknown, despite they seem not to violate 
any fundamental law of nature. On the other hand, Standard Model (SM) 
observables are compatible with experimental constraints on FCNC so far, 
with the remarkable exception of neutrino oscillations [8]. In the case of 
the lepton sector, the high supression of FCNC suggested by direct search 
is "explained" by the implementation of the Lepton Number Conservation 
(LFC), a new symmetry that protects the phenomenology from these danger- 
ous processes. Notwithstanding, the observation of oscillation of neutrinos 
[81] coming from the sun and the atmosphere, seem to indicate the existence 
of such rare couplings. Neutrino oscillations can be explained by introduc- 
ing mass terms for these particles, in whose case the mass eigenstates are 
different from the interaction eigenstates. It should be emphasized that the 
oscillation phenomenon implies that the lepton family number is violated, 
and such fact leads us in turn to consider the existence of physics beyond the 
standard model (SM), because in SM neutrinos are predicted to be massless 
and lepton flavor violating mechanisms are basically absent. These consider- 
ations motivate the stiidy of scenarios with Lepton Flavor Violation (LFV). 

The original motivation for the introduction of neutrino oscillations comes 
from the first experiment designed to measure the flux of solar neutrinos [82] , 
such measurement was several times smaller than the value expected from 
the standard solar model, so Ref. [83] suggested the neutrino oscillation 
mechanism as a possible explanation of the neutrino deficit problem. In 
addition, models of neutrino oscillations in matter [84] arose to solve the 
neutrino deficit confirmed by SuperKamiokande [85]. Since then, further 
evidence of solar neutrino oscillations has been found by SuperKamiokande 
and SNO [86]. Besides, this phenomenon can be inferred from experiments 
with atmospheric neutrinos as well [8] . 

On the other hand, since neutrino oscillations imply LFV in the neutral 
lepton sector, it is generally expected to find out LFV processes involving 
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the charged lepton sector as well. Experiments to look directly for LFV have 
been performed for many years, all with null results so far, those searches for 
these processes have only provided some upper limits, some examples are n — 
e conversion in nuclei [87], /j, — > eee [88], /x — 67 [89], and /i~ z/ee~X7^[90]. 
In addition, the search for LFV can also be made by analyzing semileptonic 
decays, upper bounds from LFV meson decays have been estimated, some 
examples arc /i+e" [91], and Kl vr" (7r+) /i+e^ [92]. 

As for experimental perspectives, a muon muon collider could provide 
very interesting new constraints on FCNC, for example /x/x fj.T{eT) medi- 
ated by Higgs exchange [21] which test the mixing between the second and 
third generations. Additionally, the muon collider could be a Higgs factory 
and it is well known that the Higgs sector is crucial for FCNC. Several au- 
thors have studied the potentiality of n/j, colliders as a Higgs factory and as 
a source of Higgs mediated FCNC [21]. 

Other interesting perspectives consist of improving the sensitivity of ex- 
periments that look for LFV processes. For example, in the case of — 67, 
the present upper bound of its branching ratio comes from MEGA collab- 
oration (1.2 X 10~^^) [98]. A new experiment is under construction at PSI 
whose aim is to improve present bound by about three orders of magnitude 
[99]. On the other hand, the SINDRUM II experiment at PSI provides the 
current upper bound of /i — e conversion branching ratio (6.1 x 10~^^) [100]. 
The MECO experiment at BNL [101] is planned to analyse such conversion 
in aluminum nuclei with a sensitivity below 10^^^. In addition, the PRISM 
project [102] aims to improve upper limits for /j, ej, and /j, — e conversion 
in nuclei, by one or two orders of magnitude. 

Another posiblc source to study LFV is by means of lepton flavor violating 
decays of a Higgs, in Higgs factories. In particular, the Fermilab Tevatron 
and the CERN Large Hadron Collider have the potentiality to study LFV 
Higgs boson decays such as /i — > r/x [64, 65]. 

From the theoretical point of view, the phenomenon of LFV has been 
widely studied in different scenarios such as Two Higgs Doublet Models, 
Supersymmetry, Grand Unification, effective theories, leptoquark models, 
technicolor, Supcrstrings and SM extensions with heavy neutrinos [93]. For 
instance, LFV processes in SU(5) SUSY models with right-handed neutrinos 
have been examined based on recent results of neutrino oscillation experi- 
ments [94]. On the other hand, Ref. [95] discuss how the supersymmctric 
SM with a see-saw mechanism could face the flavor problem; in particu- 
lar, current constraints on the leptonic soft supersymmetry breaking terms 
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and possible improvement on those constraints are studied in this reference. 
Moreover, Grimus et. al. [96] explores the generation of LFV by using a 
multi-Higgs doublet model with additional right-handed neutrinos for each 
lepton generation, finding a non-decoupling behavior of some LFV ampli- 
tudes respect to the right-handed neutrino masses. 

Further, Babu and Pakvasa [97], studied the AL = 2 process /i"*" — > 
e^dci% {i = e, fi, r), in an effective Lagrangian approach showing that it can 
explain the neutrino anomaly reported by the LSND experiment. Ref. [97], 
found two effective operators in SM that lead to such decays and no other 
processes. These operators arise from integrating out scalar fields with LFV 
interactions. 

In addition, Cvetic et. al. [42] , studied lepton flavor violation in tau decays 
induced by heavy Majorana neutrinos within two models (1) the SM with 
additional right-handed heavy Majorana neutrinos which is a typical see-saw 
model; and (2) the SM with left-handed and right-handed neutral singlets. 
The first of these models predicts very small branching ratios of the LFV 
tau decays considered, in most of the parameter space. Unlike the second 
of these models, which might show large enough branching ratios for such 
decays to be tested in near future experiments. 

Now, since there are interesting perspectives to improve the experimental 
sensitivity for the LFV processes fi ej, and n ^ e conversion in nuclei. It 
might also be interesting to check for theoretical predictions for them com- 
ing from different scenarios. Both processes might occur in many models 
with LFV. Nevertheless, the simple see-saw neutrino model does not induce 
experimentally observable rates for the // — > 67 decay [103]. But in the 
case of SUSY models the rate for this decay can be significant owing to the 
production of LFV processes by means of one loop diagrams of sleptons. Par- 
ticularly interesting are SUSY models with i?— parity broken because LFV 
interactions could arise even at tree level [104]. The branching ratios of LFV 
processes has been calculated for a variety of models such as SUSY-GUTs 
and SUSY models with right-handed neutrinos [105]. These references show 
that the branching ratios for fi 67, and fi — e conversion could get values 
near to the observable threshold in forthcoming experiments. 

The n — e conversion rate in nuclei, has been calculated by a variety of 
authors with a variety of approaches [106]. In particular, Kitano et. al. [103] 
have evaluated the n — e conversion by using the method of Gzarnecki et. 
al. [106]. These calculations were carried out for nuclei in a wide range of 
atomic numbers with an effective theory approach. In general they found 
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that the conversion branching ratios grow along with the atomic number (Z) 
up to Z ~ 30, those branching ratios become the largest for Z ~ 30 — 60, 
and decrease for heavy nuclei Z > 60. Though this tendency is general, 
the Z dependence on conversion rates is significantly different for several 
LFV couphngs. Therefore, the atomic number dependence of the conversion 
branching ratio could be useful to discriminate among the theoretical models 
with LFV. Additionally, the conversion rate also depends on proton and 
neutron densities for each nucleus, Ref. [103] uses several models for proton 
and neutron density distributions (see appendix A on that reference) , finding 
a reasonable agreement among them. 

In general, many extensions of the SM permit FCNC at tree level. The 
introduction of new representations of fermions different from doublets pro- 
duce them by means of the Z-coupling [107]. In addition, they are generated 
at tree level in the Yukawa sector by adding a second doublet to the SM [9] . 
Such couplings also appear in SUSY theories with R-parity broken [108], be- 
cause FCNC coming from R-parity violation generates massive neutrinos and 
neutrino oscillations [109]. In this work we shall concentrate on the FCNC 
arising from the general 2IIDM. For other scenarios with FCNC, I refer the 
reader to the literature [80, 107, 108]. 

On the other hand, scenarios with FCNC in the lepton sector automat- 
ically generates FCNC in the quark sector as well, experimental bounds on 
FCNC in the quark sector come from AF = 2 processes, rare B-decays, 
Z ^bb and the p-parameter [31]. Reference [31] also explored the implica- 
tions of FCNC at tree level for e+e-(//+/x-) tc + tc, t c-f{Z,g), L>° 
and B^ — Bg. Moreover, there are other important processes involving FCNC. 
For instance, the decay B^[D~) — > K~ ij,'^t~ which depends on /x — r mixing 
and vanishes in the SM. Another one is B~{D~) K~iJi^e~ whose form 
factors have been calculated in [4], [110]. 

As explained above, FCNC in the lepton sector are basically absent in 
the SM. On the other hand, in the quark sector FCNC are very tiny because 
they are prohibited at the tree level and arc further suppressed at one loop by 
the GIM mechanism [6]. Consequently, SM provides a hopeless framework 
for this rare processes since predictions from it, are by far out of the scope 
of next generation coUiders [31]. So detection of these kind of events in near 
future experiments in either the lepton or quark sector, would imply the 
existence of new Physics. 

Since FCNC are strongly constrained by experimental data, several mech- 
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anism to avoid potentially larger contribution to these exotic processes have 
been developed. For instance, Glashow and Weinberg [30] proposed a dis- 
crete symmetry in the Two Higgs Doublet Model (2HDM) which forbids the 
couplings that generate such rare decays, hence they do not appear at tree 
level. This discrete symmetry led to the so called 2HDM type I and II dis- 
cussed in section (2.2.3). Another possibility is to consider heavy exchange 
of scalar or pseudoscalar Higgs fields [4] or by cancellation of large contribu- 
tions with opposite sign. Another mechanism was proposed by Cheng and 
Sher arguing that a natural value for the FC couplings from different families 
should be of the order of the geometric average of their Yukawa couphngs 
[5] . On the other hand, the possibility of having a flavor conserving type III 
model at the tree level has been considered in Ref. [22] by demanding an 
5*3 permutation symmetry on the fermion families, see details in Sec. (3.1). 
Finally, as dicussed in Sec. (3.1), the decouphng hmit provides another nat- 
ural scenario for suppression of both FCNC and CP violating couphngs in 
the 2HDM type III. 

FCNC in the 2HDM 

Although FCNC in SM are very suppressed because of the GIM mecha- 
nism, FC transitions involving down-type quarks get enhanced because of 
the presence of a top quark in the loop. In contrast, transitions involv- 
ing up-type quarks are usually very small. Therefore in SM processes like 

Q Q 

6— > (s, d)7, K — K , B — B ; are more interesting than for example 
t — > c (7, Z) , D — D . However, by introducing FCNC at tree level (as is the 
case of the 2HDM type HI) transitions like t — >^ c (7, Z) could be enhanced 
too, giving a clear signal of new physics [31, 35]. 

In e+e~ colliders a very interesting FCNC reaction is e+e~ — > tc. Its sig- 
nature would be very clear and could compensate the lower statistics respect 
to hadron coUiders. This process can proceeds by a Higgs in the s— channel 
at the tree level, and at one loop via corrections to the Ztc, jtc vertices. 
The former is negligible and only the latter have to be considered. Defining 



and using the following parametrization for the mixing matrix elements 




tc + tc) 



a (e+e 



7* ii+ii ) 




V 



68 



Atwood et. al. [31] found that assuming A^^ = = A, values of i?*^/A^ up 
to 10~^ are reached at yfs — 400 — 500 GeV and that not much is gained for 
larger values of a/s. 

The decay i ^ 07 is highly supress in SM since the GIM mechanism 
supresses it strongly owing to the smallness of the quarks running into the 
loop (down type quarks) and also because of the large tree level rate t — > 
hW . This process as well as the processes t ^ c {Z, g) can be substantially 
enhanced in model type III respect to SM and also respect to type I and type 
II [31, 35]. 

Due to the great agreement between SM and experimental determination 
of the mass difference in the K — K , and Bd — systems, any tree level 
contribution from elementary flavor changing couplings needs to be strongly 
suppressed. This was the original motivation to impose the discrete symme- 
try given in Ref. [30] on the 2HDM. If we were to work on a model with 
FCNC a very usual assumption is that ^.y = \ij [5]. if ^e assume 

that all \ij are equal in the quark sector then AF = 2 mixing imposes se- 
vere constraints to the common A. If on the other hand, we assume that 
the mixings involving the first family are negligible then these processes do 
not place further constraints on the remaining FC couplings, allowing quite 
strong mixing between the second and third family [31]. Finally, further 
constraints coming from seems to indicate that A{,{, >> l,Att << 1 and 
\sb » l,Act << 1. i.e. only the down sector of the second and third 
family is strongly mixed. Disregarding constraints coming from i?^ the lat- 
ter assumptions are not neccesary and the assumption of negligible mixing 
with the first family is compatible with constraints given by processes like 
Br {B — > Xs'j) and the p parameter. In general, the mixing between the sec- 
ond and third generation of down-type quarks stays almost unconstrained, 
for example the contribution of new Physics for F {b ^ sec) is appreciably 
smaller than SM one for Xgf, < 40, i.e. in a large range of values. 

As for the lepton sector, FCNC (i.e. LFV processes) in extended Higgs 
sectors have been analized by several authors [4, 15, 37, 38, 39, 121]. Nie and 
Shcr examined the contribution on from Higgs mediated FCNC coming 
from a very general model with extended Higgs sector. They found that 
those contributions can be significantly enhanced if the mass of the scalar is 
light enough [15], and obtained a bound on from by assuming that the 
heavier generations have larger flavor changing couplings [5]. Kang and Lee 
[121], got constraints on some LFV couphngs (^^^, ^er, C^r, Cmm, Ce/., Cee), in 
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the framework of a general extended Higgs sector by means of the g — 2 muon 
factor and several LFV decays. These limits are gotten by setting the lightest 
scalar mass as m/jO = 100 GeV, 1000 GeV, and assuming that the other Higgs 
bosons are decoupled. On the other hand, Ref. [38] provides bounds on the 
mixing vertices {^tt, Cer, Cut, Cmm) framework of the 2HDM type III, 

by a similar approach but sweeping a wide region of the free parameters 
of the model. For instance, since the contribution of the scalars and the 
pseudoscalar have opposite signs, those limits could change significantly if 
the pseudoscalar Higgs is not decoupled. Specifically, Ref. [38] shows that 
lower limits on are obtained by using light values of the pseudoscalar 
Higgs boson, while upper limits are found by using scenarios with very heavy 
pseudoscalar, details in next chapter. Finally, Diaz et. al. [39] have found 
the constraint |C/ie/"^_ff+ 1 < 3.8 x lO^'^Gel^"^, based on the fiavor changing 
charged current decay /x v^evy,. Despite the bound is not so strong, it 
should be pointed out that it is independent on the other parameters of the 
model, so improvements on the upper limit of — > Ve'^v^ could provide very 
useful information on ^^g. 
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Chapter 4 



Flavor Changing Neutral 
Currents in the Lepton Sector 
of the 2HDM (III) 



Having motivated the study of Higgs mediated FCNC in the previous chapter, 
we shall concentrate on the lepton sector of the 2HDM type III henceforth. 
Now, since neutrino oscillations seem to indicate the existence of Lepton 
Flavor Violation (LEV) in the neutral lepton sector; we shall consider the 
possibility of having LEV couplings in the form of either ECNC in the charged 
sector, or Elavor Changing Charged Currents (ECCC) between the charged 
and neutral sector. 

For easy reference, I write the Yukawa Lagrangian type III Eq. (2.39) in 
the fundamental parametrization, for the lepton sector only 



-A 



E 



9 



-M 



diag 



E {cos aH^ - sin ah^) 



+ —=E^^E (sin aH^ + cos 
v2 

+^i''PREH+ + ^Ee^75^A° + h.c. 
v2 



(4.1) 



In order to get information on LEV couplings in the 2HDM type III, we 
shall use the following observables: the anomalous magnetic moment of the 
muon {g — 2)^ and the LEV decays: /j, — > 67, /j, — > eee, r — > eee, r — > 
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In general, the muon anomalous magnetic moment a^, is one of the most 
important observables to constrain contributions coming from new Physics. 
Since, it is perhaps the most important source to obtain our bounds, we shall 
make a brief survey on the present status of its experimental and theoretical 
estimations. 



4.1 Brief survey on (g — 2) 



The first study of (7— factors of subatomic particles was realized by Stern 
in 1921 [111], the famous Stern Gerlach experiment was done in 1924 [112] 
the conclusion was that the magnetic moment of the electron was one Bohr 
magneton i.e. that the gyromagnetic ratio g was 2. This gyromagnetic factor 
is the proportionality constant between the magnetic moment and the spin, 

the discovery that fife 7^ 2 [113] and the calculation by Swinger [114] was one 
of the great success of quantum electrodynamics. Additionally the long life- 
time of the muon allows a precision measurement of its anomalous magnetic 
moment at ppm level. The anomalous muon magnetic moment is given by 



//^ = (1 + a^) 



eh 

2^ '''' 



A series of experiments to evaluate were carried out at CERN [115]. Nowa- 
days, an ongoing experiment at BNL has been running with much higher 
statistics and a very stable magnetic field in its storage ring [116]. By com- 
bining CERN and BNL statistics the present world average result is [68] 

( 2 ) 

^exp ^ y9 - ^ (8) X 10-'° 

the ultimate goal of the E821 experiment at BNL is to reduce present uncer- 
tainty by a factor of two. 

In the framework of the SM, theoretical predictions of comes from 
QED, electroweak, and hadronic contributions 

a. = + «r + c 
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^QED constitute the bulk of the SM contributions (11 658 470.57 x 10"^°) 
and its uncertainty is small (2.9 x 10^^^) [117], the EW contribution is small 
(152.0 X 10~^^) and it is known with a similar accuracy (4.0 x 10~^^) [118], 
noteworthy, the EW contribution has the same order of magnitude as the 
current experimental uncertainty. Hadronic contributions are the second 
largest ones (~ 7000 x 10^^^) [69, 70] and their uncertainty constitutes the 
bulk of the uncertainty contributions (~ 70.7 x 10^^^). 

As for the hadronic contributions, the most important ones in terms of 
both contributions to and its uncertainties are the leading vacuum po- 
larization terms, they have been estimated recently by a set of authors [70] . 
The first of Refs. [70] has used the most recent data on hadron production 
from e+e" colhsions from BES, CMD-2, and SND [119]. The HMNT group 
(see fourth of Refs. [70]) has made estimations based on both inclusive and 
exclusive analysis of e+e" ^hadrons. Combining such exclusive and inclu- 
sive sets of data, they obtain a more optimum estimate, and evaluate QCD 
sum rules in order to solve a discrepancy between the inclusive and exclusive 
analysis. In addition, the second and third Refs. [70] presents an analysis 
based on e^e~ data as well as r hadronic decays data. All of Refs. [70] found 
that analysis based on e^e^ — >hadrons, predicts a SM value which is roughly 
3(7 below the present world average measurement, however, the second and 
third of Refs. [70] have showed that data based on hadronic r decays only 
differs with the experimental world average by about I.60". Stabilization of 
the predictions coming from e+e^ based data and r decays data is neces- 
sary in order to say something definitive. For our purposes we shall use the 
estimation given by Jegerlehner [70] whose analysis is based on e'^e~ data 
by taking the most recent results from CMD-2 [119]. Furthermore, other 
important contribution is the light by light scattering whose sign has been 
corrected recently [69]. 

On the other hand, it is important to emphasize that despite the mea- 
surement of the muon anomalous magnetic moment of the muon is about 
180 times less accurate than the electron anomalous magnetic moment; the 
former is more sensitive to new physics by a factor of ~ (m^/mg) , this fac- 
tor more than compensates the difference in precision. Therefore, is more 
important than in looking for new physics. 

In constraining the 2HDM type III we shall use the estimation for Aa^ at 
95% CL reported by [72] 

9.38 X 10"^° < Aa^ < 51.28 x 10"^° (4.2) 
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which in turn are based on prehminary results by Jegerlehner [70] . 



4.2 The {g - 2)^ factor in the 2HDM type III 

The first one loop electrowcak corrections for coming from new Physics 
(Aa^) were calculated in Refs. [120]. For Aa^ the integral expression [17] is 
given by 



'5(P) = 



r2 ^/^ 



P 
2 ^1 



X 



Stt^ Jq m^x^ + (m^ — m^) x + M'^^p^ (1 — 



X] 



-dx (4.3) 



where Is{p) is an integral involving an Scalar (Pseudoscalar) Higgs boson 
with mass Ms(p), and Cs{p) is the corresponding coefficient in the Yukawa 
Lagrangian Eq. (4.1). If we assume that << and << m^o jjo y^o in 
the calculation of (4.3), the one loop contribution from all neutral Higgs 
bosons reads [37] 



Aa,. 
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^Cll^FiMs)-^G{Ms 
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J2cl(FiM,) + ^G{M, 



(4.4) 



where 



G{M) = 



F{M) 



M 



2 + 3M2 + 6M2 In 
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We have neglected the contribution of the charged Higgs boson because of 
two reasons [37]: on one hand this contribution involves the neutrino mass 
and on the other hand, the CERN e'^e~ collider LEP bound on its mass is 
mH+ > 80.5 GeV. 



4.3 Lepton Flavor Violating decays 

As we explained below, some lepton flavor violating decays are useful to 
constrain the 2HDM type 111. The relevant expressions of the lepton decays 
involved at leading order are given by 



r (/X ^ 67) = (mj/O , m/,0 , m^o , a, ^rr 



(4.6) 



where I = e, /i denotes a light charged lepton. In addition, we have defined 

v^^rr J^iii'^ a \ \n[rii jjii / iiIt] 



R {tuho , niho , tuao , a, Crr) 



TTir sin 2a + 



m 



HO 



— I uIt sin 2a — 
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2^irr ln[m^o/ 
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sin q; ^5 h cos a ^ 
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In [m^o / m^] 



A" 



(4.7) 



On the other hand, the expression for a lepton L going to three leptons of 
the same flavor I is given by 



r (L ^ III) 



m. 



20 487r3 
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sin2Q;^ 
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V2\ 



mi 



sin^ a cos^ a 
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HO 



H 2 



(4.8) 
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Figure 4.1: Leading order diagrams for the processes: {g — 2)^ (top left), 
L ^ Ij (top right), L — >• III (bottom left), and /i~ — >• 'de'df^.e'^ (bottom right). 
In the two former, we have neglected the contribution of the charged Higgs 
and a neutrino into the loop. 

Finally, the decay width for /i~ — > z/gC"!/^, reads 



r U- ^ v,e-V^ = ^ (4.9) 

' 24 5767r3 Vmjy+y ^ ' 

the latter is a process with Flavor Changing Charged Currents (FCCC) me- 
diated by a charged Higgs boson as we see in figure (4.1). 

In calculating Aa^ and the leptonic decays of the type L Z7, we neglect 
the contribution of light leptons into the loop, so only the contribution with 
a tao into the loop is considered, see figure (4.1). 



4.4 Obtaining the bounds for the flavor chang- 
ing vertices 

From the previous section, we see that the free parameters that we are in- 
volved with, are: the three neutral Higgs boson masses (m/jO, mji/o, m^o), the 
mixing angle a, and some flavor changing vertices ^ij. Based on the present 
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Bounds on 


Bounds on 


Bounds on 


case 1 


7.62 X 10-^ < il, < 8.31 X 10-3 


ee ^ 7.33 X 10-18 


^er < 4.82 X 10-1^ 


case 2 


1.29 X 10-3 < ^ 4.42 X 10-2 


^ 2.24 X 10-16 


e < 2.77 X 10-14 


case 3 


1.53 X 10-3 <e^. 


< 2.24 X 10-16 


e < 2.76 X 10-14 


case 4 


9.57 X 10-4 < < 1.40 X 10-2 


ea< 2.10x10-17 


< 8.22 X 10-15 


case 5 


1.02 X 10-3 < < 1.66 X 10-2 


< 2.93 X 10-17 


e < 9.65 X 10-15 



Table 4.1: Constraints on the mixing parameters ^gr^^r Cer ^'^^ ^^e 
five cases mentioned in the text. The two former are generated from Aa^ and 
r (/X — > 67) respectively, while the latter comes from the combination of the 
lower hmit on and the upper bound on Cer^^T- 

bounds from LEP2 we shall assume that 171^0 ^115 GcV. In addition, we 
shall assume that m^o > rriho, both assumptions will be held throughout this 
section. Now, since we are going to consider plots in the ^ij — ttiao plane, we 
should manage to use appropiate values of (m^/o, a) in order to sweep a wide 
region of parameters. To sweep a reasonable set of this pair of parameters 
we utilize for rnno values of the order of 115 GeV (light), 300 GeV (interme- 
diate), and very large masses (heavy). As for the angle a, we consider values 
of q; = (minimal mixing), a = n/A (intermediate mixing), and a = 7r/2 
(maximal mixing). It can be seen that all possible combinations of (mjjo, a) 
could be done by considering five cases 1) when tuho — 115 GeV; 2) when 
rriHO — 300 GcV and a = n/2; 3) when rriH" is very large and a = 7r/2; 
4) when ttiho ~ 300 GeV and a = 7r/4; 5) when 171^0 is very large, and 
a — 7r/4. In all these cases the mass of the pseudoscalar will be considered 
in the range 115 GeV< m^o. 

It is worthwhile to remark that these five cases cover all possible combina- 
tions of {ttiho, a), with the three different values for uiho (115 GeV, 300 GeV, 
and very heavy) and the three different values for a (0, 7r/4, 7r/2), though 
at the first glance some cases seem to miss. For example, in the first case 
we have uiho — rriho. In this scenario, all the processes that we calculate be- 
come independent on a and therefore no distinction is neccesary among the 
three different values of it. The cases a) tjiho = 300 GeV, with a = 0; and 
b) rrifjo very large, with a = 0; are also included in case 1, because when 
a = the processes in the previous section does not depend on ttiho. 
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The first bounds come from Aa^. We use tlie estimated value of it, given 
by [72] at 95% C.L. Eq. (4.2). Since Aa^ in Eq. (4.2) is positive, lower and 
upper bounds for the PC vertex can be gotten at 95% C.L. The results 
are indicated in table (4.1) column 1. The lower bounds in each case are 
obtained when m^o 115 GeV, and using the minimum value of Aa^ in Eq. 
(4.2), while the upper bounds arc obtained when A'^ is very heavy, and using 
the maximum value of Aa^ in Eq. (4.2). From these results a quite general 
and conservative allowed interval can be extracted [38] 

7.62 X 10-^ < < 4.44 x 10"^ (4.10) 

Furthermore, upper bounds for the product obtained from the 

expression of the decay width F {/j, —>■ ej) in Eq. (4.6) and from the experi- 
mental upper hmit F (// ^ 67) < 3.6 x 10"^° GeV [122]. The most general 
upper bounds arc obtained for very heavy. The results are shown in table 
(4.1) column 2. From this table we infer that quite generally the bound is 
[38] 

^ 2.24 X 10-16 (4.11) 

Moreover, combining these upper limits with the lower bounds on 
given in the first column of table (4.1), we find upper limits on ^g^. The 
results appear on table (4.1) third column, and the general bound can be 
written as 

< 2.77 X 10-1^ . (4.12) 

Noteworthy, these constraints predict a strong hierarchy between the mixing 
elements and ^er i-e. |^er| << iCt^rl and they differ by at least five orders 
of magnitude. 

On the other hand, from Eq. (4.6) we see that the decay widths F (r — > jj.'j) 
and F (r ^ A*A*A*) depend on two mixing vertices and ^^^^ respec- 

tively. Then, we can get conservative constraints on the diagonal mixing 
vertices ^TT, Cii/i by using once again the lower bounds on obtained from 
Aa^. Since F (r — fij), and F (r — ^ /i/x/x) are rather complicate functions of 
^rr, and ^fj_fj_ respectively, we present these constraints in the form of contour- 
plots in the m^o — ^^-^ plane and the m^o — plane, figures 4.2 and 4.3, 
each one for the five cases. 

We can check that for each contourplot there is a vahic of m^o for which 
^T-^ or stays unconstrained, and they are shown in tables 4.2, and 4.3 
respectively. Additionally, the bounds for ^^-t, when m^o is very large and 
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Figure 4.2: Contourplots in the — "Iao plane for each of the five cases 
cited in the text. On left: case 1 corresponds to the long-dashed line, case 
4 to the short-dashed line, and case 5 to the solid line. On right: case 2 
corresponds to dashed line, and case 3 is solid line. 




Figure 4.3: Contourplots in the — m^o plane for each of the five cases 
cited in the text. On left: case 1 corresponds to the long-dashed line, case 
4 to the short-dashed line, and case 5 to the solid line. On right: case 2 
corresponds to dashed line, and case 3 is solid line. 
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Values of m^o for 
unconstrained 


^T-T intervals for 
m^o very large 


^T-T intervals 

for m^o f« 115 GeV 


case 1 


m^o ^ 170 GeV 


-0.0072 < < 0.0072 


-0.010 < ^rr < 0.010 


case 2 


m^o ~ 440 GeV 


-0.018 < irr < 0.018 


-0.0043 < Crr < 0.0043 


case 3 




unconstrained 


-0.0036 < < 0.0036 


case 4 


m^o « 228 GeV 


-0.0075 < < 0.022 


-0.0094 < < 0.0036 


case 5 


m^o 250 GeV 


0.00024 < Crr < 0.021 


-0.0093 < Crr < 0.0026 



Table 4.2: Bounds extracted from the contourplots shown in Fig. 4.2. The 
first column indicates the values of m^o for which stays unconstrained for 
each of the five cases. The second and third columns show the allowed inter- 
vals on when m^o is very large and when m^o fa 115 GeV respectively. 



when rriA" ~ 115 GeV are also included in tables 4.2, 4.3 respectively. We 
see that the general constraints for read 

-1.8 X 10~^ < ^ 2.2 X 10~^ , for niAo very large 

-1.0 X 10~^ < ^rr < 1-0 X 10"^ , for niAO ^ 115 GeV (4.13) 

These constraints are valid for all cases, except for the third case with 
tuao very large, since in that scenario remains unconstrained. Now, 
for the general bounds read 

I^mmI ~ 0.12, for niAO very large 

\^^^\ < 0.13, for m^o f« 115 GeV (4.14) 

Once again these constraints are not valid for the third case when m^o is 
very large, but are valid in all the other cases. 

We should remember that the diagonal mixing vertices generate differ- 
ences among the relative couplings in model type 111, and that this lack of 
universality could be a distinctive signature for it. This fact might leads us 
to scenarios in which a certain Higgs boson (say /i°) may be tau-phobic-*^ but 
not muon- phobic or electron- phobic (see section 2.2.3). 

^With the tau-phobic limit for certain Higgs, we mean that the coupling of that Higgs 
to a pair of tau Icptons vanishes. Nevertheless, the couplings of such Higgs to a tau and 
any other lepton could exist. 
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Values of m^o for 
^fjfj unconstrained 


^^p, intervals 
for iriA » rriho 


intervals 
for m^o 115 GeV 


case 1 


m^yo ^ 115 GeV 




< 0.0 13 


unconstrained 


case 2 


TTiAo ^ 300 GeV 




<0.12 




< 0.055 


case 3 




unconstrained 




< 0.043 


case 4 


niAO ^ 152 GeV 




< 0.058 




<0.13 


case 5 


TTiAO ~ 163 GeV 




< 0.061 




<o.ii 



Tabic 4.3: Bounds extracted from the contourplots shown in Fig. 4.3. The 
first column indicates the values of m^o for which stays unconstrained for 
each one of the five cases. The second and third columns show the allowed in- 
tervals on when m^o is very large and when m^o 115 GeV respectively. 

We also ought to remember that in model type III, the pseudoscalar Higgs 
boson couples to fermions only by means of the matrix element ^jj, while in 
models I and II it couples through the mass of the corresponding fermion. 
Therefore in model III, if C// = for a certain specific fermion, then A'^ becomes 
fermiophobic to it (and only to it). 

At this step, it is important to remark that our bounds are compatible 
with all tau-phobic limits for each Higgs boson. It is easy to check that 

becomes tau-phobic if = for a = 0, and if C,tt — 0.01 for a = n/A; for 
a = 7r/2 there is not tau-phobic limit. On the other hand, becomes tau- 
phobic if ^rr = for a = 7r/2, and if ^rr — —0.01 for a — 7r/4; and there is 
no tau-phobic limit for a — 0. The tau-phobic limit for the pseudoscalar is 
the simplest one, = 0. Finally, if = and a = {■k/2) then and 

(H^) become tau-phobic simultaneously. 

On the other hand, bounds given above could in contrast produce a sig- 
nificant enhancement to the Yukawa couplings. As a matter of example, in 
the case 5 with a — n/A using 0.1 ^ we see that the contribution com- 

ing from the term proportional to ^rr, is about 10 times larger in magnitude 
than the contribution coming from the term proportional to the tan mass, 
for both the and couplings. So it is even possible for the couplings of 

^In the case 5 both values = 0.01 and = 0.1 arc allowed at least for m^o around 
250 GeV, as it is shown in Fig. (4.2). However, for a very light (m^o « 115 GeV), or a 
very heavy pseudoscalar Higgs boson (m^o >> m^jo), constraints on ^^-^ are considerably 
stronger in all cases, see table 4.2. 
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the CP— even Higgs bosons to be overwhelmed by the ^^-^ contribution. Of 
course, since the bounds for are weaker than the constraints on E,tt as 
we can see from Fig. 4.3 and table 4.3, we can also have highly supressed or 
highly enhanced Yukawa couplings to a pair of muon fermions. 

Based on the bounds obtained above, we are able to estimate upper limits 
for some leptonic decays by using the expressions (4.6,4.7,4.8). We shall 
assume that l^rA ^ 0.1, and that |^ee| i$ 0.1; from these assumptions we can 
infer the following upper bounds [38] 

r(r^e7) < 1.5xl0-2^ 

r (r ^ eee) < 5 x 10"^^ (4.15) 

If we compare with the current experimental upper bounds, [122] 

r(r ^ 67) < 6.12 X IQ-^^ 
r (r ^ eee) < 6.57 x 10"^^ 

We see that these decays are predicted to be very far from the reach of next 
generation experiments in the context of the 2HDM type III, unless that 
^rr, and/or ^^e acquire unexpectedly large values, this remarkably strong 
supression owes to the dependence of these decays on the mixing vertex ^er- 
For the sake of comparison, Ref. [42] has obtained upper limits for these 
decays in the context of two models with heavy Majorana neutrinos: (I) the 
SM with additional right-handed heavy Majorana neutrinos i.e. a typical 
see-saw type model and (II) the standard model with right-handed and left- 
handed neutral singlets. The first of these models predicts very small decay 
widths for LFV processes in most of the parameter space, but the second 
one might show large enough decay widths. Namely, the upper limits in the 
context of the second of these models read [42] 

r (r ^ e7) < 2.27 x 10"^° GeV , 
r (r ^ eee) < 2. 27 x 10^^^ GeV. (4.16) 

which arc not so far from the present experimental threshold and could be 
tested by near future experiments, unlike the case of the 2HDM (III). Such 
fact could be useful to discriminate between these models. 

On the other hand, the decay /x~ i^ee"!?^ provides information about 
the mixing between the first and second lepton family. The leptonic pro- 
cesses analyzed so far, involves FCNC mediated by neutral Higgs bosons. 
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Instead, the process n~ — >■ z/ge"!/^, involves Flavor Changing Charged Cur- 
rents mediated by a charged Higgs boson. Nevertheless, from Lagrangian 
(4.1) we see that the matrix elements ^fj that generates FCNC automat- 
ically generates FCCC which are strongly suppress in the leptonic sector. 
Consequently, by constraining flavor changing charged currents we are indi- 
rectly constraining flavor changing neutral currents as well. In the case of 
the 2HDM (III), interactions involving FCCC at tree level only contains the 
contribution from the charged Higgs boson, reducing the free parameters to 
manage. Thus, from the decay width F (/x ^ Uf.e~Vn) we are able to extract 
a bound for the quotient ien/i^H+ independent on the other free parameters 
of the model. Taking the analytical expression for the decay width Eq. (4.9) 
and the current experimental upper bound for it [122] 

F ^ T^ee'V^,) < 3. 6 X 10-^^ GeV, (4.17) 

the following constraint is gotten 

< 3.8 X 10^^ GeV-^ (4.18) 

Despite this constraint is not so strong, it is interesting since it does not 
depend on the other free parameters of the model, because the calculation 
does not involve neutral Higgs bosons nor mixing angles. This is a good 
motivation to improve the experimental upper limit for processes involving 
FCCC in the leptonic sector. 

4.5 Bounds on the Higgs boson masses 

In the previous section, we obtained constraints for the mixing vertices based 
on some sets of values for the Higgs masses. Conversely, we can assume a 
set of values for the flavor changing vertices and try to get constraints on the 
Higgs boson masses. We shall do it in the framework of the new parametriza- 
tion developed in Ref. [36] and in section (2.2.3) of this work. We should 
be very careful in interpreting such results because a different basis has been 
used for the Yukawa Lagrangian. For instance, we shall show allowed regions 
in the rriAQ — tan/3 plane by using flxed values for mj^o, m/jo. However, 
we ought to bear in mind that S,^t depends explicitly on tan/5 as we see in 
Eqs. (2.50), checking Eqs. (2.50) we see that the only way to keep fixed 
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with changing tan/3, is by varying the ^^j- parameter of the fundamental 
parametrization; applying the last of Eqs. (2.50) to the vertex we get 

= (Vl + tan^/^) (4.19) 

so by using a fixed value of and varying tan /3 (say from OtoN), what we 
are really doing from the point of view of the fundamental parametrization, 
is sweeping from to Vl + N^^^r ■ In addition, we should take into 
account that a' = a + (3; Therefore, if we want to keep a' constant as well, 
a should be varied accordingly. Consequently, to vary tan/3 while keeping 
fixed and a', it is necessary to vary the parameters and a in the 
fundamental parametrization^. It should be pointed out then, that changing 
tan/5 would imply a correlation among the variation of a and C.^^- which is 
not a necessary condition. A more appropiate picture would be fixing tan (3 
(the basis) and use the other parameters as variables. Notwithstanding, we 
shall show some plots in the m^o — tan P plane, assuming that correlation. 

The possibility of plotting in a m^o — tan f3 plane for the model type III, 
facilitates the comparison of the model type III with the models type I or II. 

With this clarification we proceed to constrain the FC vertex involving 
the second and third charged leptonic sector by using the estimated value for 
Aa^^. Additionally, we get lower bounds on the Pseudoscalar Higgs mass 
by taking into account the lower experimental value of Aa^^ at 95% CL 
reported in [72] and making reasonable assumptions on the FC vertex. 

For easy reference, I include the Yukawa Lagrangian type III for the 
leptonic sector by using parametrizations of type I and II explained in section 
(2.2.3). 

-£y(f^ = TTVy^ -EM^'^'E(sma'H^ + cosa'h^) 

+ _^EMf'-i^EG^ + ':^^^EMf'-i^EA^ 
^ -E^^E [sin {a' - P) + cos {a' - (3) 



V2sin/5 

—^E^^^.EA^ + h.c. (4.20) 
V 2 sm (i 



^We can assume another point of view in which tan /? and a' are changed in such a way 
that only varies in the fundamental parametrization. But it implies to work with an 
additional variable in the non-trivial parametrization. 
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where the superindex (/) refers to the parametrization type I. It is easy to 
check that Lagrangian (4.20) is just the one in the 2HDM type I, plus some 
FC interactions. Therefore, we obtain the Lagrangian of the 2HDM type 
I from Eq. (4.20) by setting rj^ — 0. In this case, it is clear that when 
tan /3 — > then rj^ should go to zero, in order to have a finite contribution 
for FCNC at the tree level. 

In the parametrization of type II, the Yukawa Lagrangian reads 

2Mw cos 13 ^ ^ ^ 

+ eI'^E [sin {a' -(3)H'^ + cos {a' - (5) hP] 

V2 cos/3 

+ ' E^'^^.EA' + h.c. (4.21) 
v2cos/5 

The Lagrangian (4.21) coincides with the one of the 2HDM type II [17], plus 
some FC interactions. So, the Lagrangian of the 2HDM type II is obtained 
setting = 0. In this case it is clear that when tan (3 — > oo, then C,^ should 
go to zero, in order to have a finite contribution for FCNC at the tree level. 

We now use the expression for Aa^ from either Eq. (4.3) or Eq. (4.4) to 
get bounds on the Higgs boson masses. Additionally, we can notice that in 
Eq. (4.5) for Higgs bosons heavier than 100 GeV we find 

777 



3m 



r 



and the contribution of the factor involving G (niHi) is negligible [37]. How- 
ever, we put it for completeness. 

In Ref. [37] the estimated value for Aa^ has been taken from [116]. In 
this section we are going to use the basic procedure developed in Ref. [37], 
but using a more updated value for Aa^ i.e. the one reported in Ref. [72], 
instead of the one reported in [116]. We get some lower and upper bounds 

on the mixing vertex rj for the parametrizations of type I (II). In Fig. 

(4.4), we display lower and upper bounds for the FC vertices as a function 
of tan/3 for both types of parametrizations with m/^o = mno = 150 GeV 
and m^o cxd. In the first case, with parametrization type I, the allowed 
region for rj^r is 0.045 < rjf^r ^ 0.12 for large values of tan/3. Meanwhile, 
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Figure 4.4: Lower and upper bounds for rj^^ tan/?, for parametriza- 

tions I and II using m/jO = rriHO = 150 GeV and m^o —>■ oo. 
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Figure 4.5: Lower and upper bounds for ^^t- vs m^/o, for parametrization of 
type II, taking ruho — mjjo and m^o — > oo, the pair of dotted lines correspond 
to tan /3 = 0.1, the dashed hnes are for tan (3—1, and the sohd hnes are for 
tan 13 = 30. 

for parametrization type II, the allowed region for small tan/5 is the same. 
From the first of Figs. (4.4) we can sec that when tan/? — > 0, rj^j. goes to 
zero as expected. Similarly, the second of these figures, shows that when 
tan/3 — > oo, ^^^r goes to zero as expected too. 

In addition, the Figs. (4.4) show us that the FC vertices are basis depen- 
dent. We can realize that by varying tan/5, what we are doing is changing 
the basis. Figs. (4.4) show that for different values of tan/3 (i.e. different 
bases) we get different values for the lower and upper limits on rj^j. ( j . 



In figure (4.5), wc show lower and upper bounds for the FC vertex as 
a function of m//o for parametrization of type II when m/jo = mj^o and 
m^o — > oo. We see that the larger value for tan/3 the smaller value of 
. We only consider the case of parametrization type II because there is a 
complementary behaviour between both parametrizations as could be seen 
in figure (4.4). In particular, for tan/3 = 1, the behaviour of the bounds for 
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both parametrizations is the same. Since in this case we are fixing tan/?, we 
are fixing the basis and the interpretation is clear. It should be emphasized 
that according to Eq. (4.19) for different values of tan/?, the same value of 

correspond to different values of .^^^ in the fundamental parametrization. 
Thus, we see again that the values of the FC vertices are basis dependent. 

Observe that according to its current estimated value, Aa^^ is positive 
definite, and it is a new feature from most updated results [70]. On the 
other hand, the Feynman rules from (4.20) and (4.21), show that the Scalar 
(Pseudoscalar) contribution to Aa^^ Eqs. (4.3, 4.4) is positive (negative). 
Such fact permits us to impose lower bounds on the pseudoscalar Higgs mass, 
by using the lower limit for Aa^^. However, it is very important to clarify 
that for very light Pseudoscalars (i.e. my^o << 100 GeV), the contributions 
up to two loops become important [72]. Further, the diagrams at two loops 
for the scalars and the pseudoscalar could have opposite signs respect to the 
contribution at one loop. Consequently, the lower bounds shown here are 
not valid in the very light pseudoscalar regime, and they are only valid if we 
assume that m^o > 100 GeV. 

Now, to take into account the experimental value, we should make a 
supposition about the value of the FC vertex. A reasonable assumption 
consists of taking the geometric average of the Yukawa couplings [5] i.e. 



2.5 X 10 ^ and V \^) ~ 2.5 x 10 ^ which are one order of magnitude larger 



and smaller than the former. Notwithstanding, we should remember that 
these values are basis dependent. Therefore, in this work we shall take the 
values yielded above for the FC vertex among the second and third family by 
using different bases, i.e. different values of tan /3. Using these suppositions 
and the experimental vahie for Aa^^ reported in Ref. [72] we get lower 
bounds for m^o and they are plotted in figures (4.6-4.8). 

Figure 4.6 displays m^o vs tan j3 utilizing the parametrization type II and 
using = 2.5 X 10~^, and ^p,- = 2.5 x 10~^. Additionally it is assumed 
that rriho = mjjo with m^jO = 115, 300 GeV*. According to the interpretation 
explained above, since we are using < tan/3 < 100, it is equivalent to 

''By using ^^j- = 2.5 x 10~^ we still have an allowed region but it appears only for 
tan/3 > 155 and for tan/3 > 360 when niho = 115 GeV and when m/jo = 300 GeV 
respectively. The behavior of those curves is similar to the others and is not shown in Fig. 




2.5 X 10~^ 



Additionally, we shall use also the values rj ( ^ 




4.6. 
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Figure 4.6: Contourplot of m^o vs tan/? using parametrization type II and 
assuming m/jO = m^^o. From left to right: (1) m/jO = 115 GeV, = 
2.5 X 10-2 (2) rriho = 300 GeV, ^i^r = 2.5 x lO'^ (3) m^o = 115 GeV, = 
2.5 X 10-3 (4) 300 GeV, = 2.5 x 10"=^. 
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use a range of from to 100.^^t- in the fundamental parametrization. 
It could be seen that in the limit of large tan/9, the lower limit reduces 
to m^o m/jO . The same behavior can be seen in parametrization type I but 
the bound m^o m/jO is gotten in the limit of small tan p. We also see that 
the smaller value of the stronger lower limit for m^o . 

As a proof of consistency, we can see that for instance in the case r7i/,() = 
niHO = 115 GeV with = 2.5 x 10^^ the allowed region starts at about 
tan/3 fti 13. If we calculate the value of in the fundamental parametriza- 
tion corresponding to tan/3 ^ 13 we obtain from Eq. (4.19) that ^"^^ a; 
1 X 10~^, in rough agreement with the minimum allowed value for re- 
ported in table (4.1) for the case 1^. 

Figure 4.7 shows m^o vs rriHO with = 2.5 x 10~^ and tan/3 = 2, 
using m^o = 111^0 and m^o = 115 GeV. With these settings, the values 

= 2.5 X 10~^ and = 2.5 x 10"^ are excluded. 

Once again, we can see the consistency with the analysis in Sec. (4.4) 
noting that for = 2.5 x 10~^ and = 2.5 x 10~^ with tan /3 = 2, we get 

= 3. 125 X 10~^ and = 3. 125 x 10"''' respectively. Looking at the case 
1 in table (4.1) we see that these values are not allowed, and hence they are 
excluded in Fig. 4.7. By contrast, the value C,^r = 2.5 x 10^^ with tan/? = 2, 
leads to = 3. 1 x lO^'^ which is clearly allowed according to the table (4.1) 
in the case 1. 

In figure 4.8 we suppose that m^^o = 115 GeV, mno = 300 GeV. The figure 
above shows the sensitivity of lower bounds on m^o with the mixing angle 
a', for parametrization type II, taking tan/3 = 20. The value ^^t- = 2.5 x 
10~^ is excluded again. The constraints are very sensitive to the a' mixing 
angle for ^^t- = 2.5 x 10~^ but less sensitive for = 2.5 x 10~^. The 
figure below shows ttia" vs tan/3 for rriho = 115 GeV, m^-o = 300 GeV, 
a' = 7r/6, for parametrization type II and considering the C,fj,T = 2.5 x 10~^ 
and ^^T- = 2.5 x 10^'^. The value C,fj,T = 2.5 x 10~^ is excluded. The m^o lower 
asymptotic limit for large tan /3 is approximately m/jO . 



^The difference owes to the fact that here we have used the approximate expression for 
Attfi given in Eq. (4.4), instead of the Eq. (4.3) which was the one used to obtain the 
table (4.1). Additional differences could be traced to the errors involved in the calculation 
of the contourplots. A more careful numerical analysis could improve the agreement, but 
it is beyond the purpose of the present work. 
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Figure 4.7: Contourplot of m^o vs tuho setting tan /3 = 2, the line above cor- 
respond to — 2.5 X 10~^ and m/jo = mijo, while the line below correspond 
to ^i,r = 2.5 X 10-2 for rriho = 115 GeV. 
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Figure 4.8: (top) Contourplot of m^o vs a', for parametrization type II, with 
m/jO = 115 GeV, ttiho = 300 GeV and tan /3 = 20. Line above corresponds to 
= 2.5 X 10~^ , and line below corresponds to = 2.5 X 10"^. (bottom) 
Contourplot of m^o vs tan (3 for m^o = 115 GeV, uiho — 300 GeV, a — n/6, and 
for parametrization type II. Line above corresponds to — 2.5 X 10^^, and line 
below corresponds to — 2.5 X 10^^. 
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Chapter 5 



Addendum: Top-squark 
searches at the Fermilab 
Tevatron in models of 
low-energy supersymmetry 
breaking 



5.1 General framework 



This chapter describes the work developed at the Fermi National Acceler- 
ator Laboratory, along with M. Carena, D. Choudhoury, H. E. Logan, and 
C.E.M. Wagner [41]. It is about the perspectives for detection of the lightest 
top squark at the Tevatron collider in the case in which such stop comes from 
a supcrsymmetric model with low energy supersymmetry breaking pattern. 
Since it is an addendum, I will provide only a very brief survey on the frame- 
work of supcrsymmetric models and supersymmetry breaking schemes. For 
more details 1 refer the reader to the literature [123, 124] 

One of the greatest problems of the SM is the so called "naturalness prob- 
lem" . The radiative corrections of the Higgs mass squared depend quadrat- 
ically on the cut off scale energy used. If one hopes to have a Higgs mass 
of the order of the EWSB scale, it is required a very precise fine tuning (of 
about one part in 10^^), in order to get the cancellation necessary at the 
Planck scale, the apparent necessity of such extremely accurate cancellation 
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is known as the naturalness problem. On the other hand, it is well known 
that such radiative corrections with fermions and bosons into the loops have 
opposite signs. Consequently, if we had a fermion of the same mass for each 
boson and vice versa, the necessary cancellation would occur. Supersymme- 
try provides a framework in which each SM particle posseses a superpartner 
i.e. a particle with the same mass and quantum numbers under the gauge 
symmetry but with different spin. Therefore, imposition of SUSY provides 
an elegant solution to the naturalness problem. 

Supersymmetry consists of an extension to the Poincare group by the 
introduction of new generators. The algebra of the group is determined by a 
set of conmuting and anticonmuting relations among the poincare generators 
and the new generators. The new supersymmetric generators permits to 
relate fermion fields with boson fields providing a unified description of them. 
In order to keep the supersymmetry, the spectrum of the SM should be 
doubled arising the so called superpartners. 

Nevertheless, since the superpartners have not been detected yet, they 
cannot have the same mass as their corresponding SM partners. Conse- 
quently, supersymmetry must be broken in such a way that the superpart- 
ners should be heavier than their partners. However, such breaking have to 
preserve the cancellation of quadratic divergencies that was one of the orig- 
inal motivations, this fact imposes constraints on the splittings among the 
partner and superpartner masses, in general such splitting cannot be larger 
than a few TeV. As a consequence, some of the lightest superpartners could 
have masses within the reach of accelerators like the Tevatron or the LHC. 
As we shall see later, the hghtest stop (i.e. one of the the superpartners of 
the top quark), is a good candidate to be one of the lightest supersymmetric 
particles. Before discussing that point, let us describe briefiy the particle 
content of the Minimal Supersymmetric Standard Model (MSSM) , and some 
breaking schemes of supersymmetric models. 

5.1.1 The MSSM 

In the MSSM, the particle spectrum of the SM is doubled. On the other hand, 
in order to get cancellation of anomalies and preserve SUSY, two Higgs dou- 
blets must be introduced. The Higgs spectrum is similar to the one of the 
general 2HDM i.e. two CP-even scalars {H^.h'^), one CP-odd scalar 
and two charged Higgs bosons {H"^). The Higgs bosons and the SM parti- 
cles have superpartners with the same quantum numbers under the gauge 
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symmetry but with different spin. The superpartners of the gau^ bosons 
(gauginos) are spin 1/2 particles denoted as photinos (7), Winos (W), Zinos 
(Z) and gluinos (g). The fermionic superpartners of the Higgs bosons (Hig- 
gsinos) are Hi, H2, H"^. Owing to the EWSB, the Higgsinos and the gaugi- 
nos of the electroweak symmetry can be mixed to give the mass eigenstates, 
arising two charged Dirac fermions called charginos x^, and four neutral ma- 
jorana fermions (the neutralinos Xi_4). By contrast, the gaugino of the color 
symmetry (gluino g) does not mix with Higgsinos nor electroweak gauginos, 
because it belongs to an unbroken symmetry. Finally, the spin-0 partners 
of the fermion fields (the sfermions) are the squarks q, charged sleptons 
I and sneutrinos u. There are two superpartners associated to each quark or 
charged slepton, one for each chirality, the left-handed sfermions transform 
as SU (2)^ doublets while right-handed ones are SU (2)^ singlets, there are 
only left-handed sneutrinos because of the lack of right-handed partners in 
the SM. The gluino is a SU (3)(^ octet while the squarks are SU (3)^ triplets. 

As well as the interactions of SUSY particles with SM particles, the La- 
grangian of Supersymmetry should contain some soft SUSY breaking terms. 
The term "soft" means that these terms breaks SUSY and generates split- 
tings among the SUSY particles and their partners but maintaining the 
gauge symmetry and the cancellation of the quadratic divergencies mentioned 
above. The soft SUSY breaking parameters are mass terms for gauginos and 
sfermions as well as trilinear scalar couplings. The number of those kind of 
terms depend on the specific SUSY breaking scheme. 

The chargino and neutralino masses and mixing angles are determined by 
the gauge boson masses M\y, Mz, the parameter tan/?, the SUSY Higgsino 
mass parameter /i and the two soft breaking parameters associated to the 
SU (2)^ gaugino mass (Mi) and the U (1) gaugino mass (M2), evaluated at 
the EW scale. The neutralino mass matrix in the B — — Hi — H2 basis 
is written as: 

2 _ / —Mz cos P sin 9w Mz sin (3 sin 9w 
y Mz cos (3 cos 9w —Mz sin j3 cos 9w 

The values of the parameters Mi, M2, /x, tan/3 determines the Higgsino 
and gaugino content of the neutralinos. There are some interesting limits in 
which the Higgsino and gaugino components become simpler. As a manner of 
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example, if >> Mz and Mi,M2^ Mz, with Mi < M2, the hghtest neu- 
trahnos are gaugino-hke and the heaviest are Higgsino hke, the eigenvalues 
of the matrix Mat in Eq. (5.1) yield the following spectrum 

M^o Ml ; Mrjo Ri Mo ; M^o M^o lul 

In particular, we shall work under the assumption that the lightest neutralino 
is a pure bino Xi ~ details later on. The Higgsino and gaugino compo- 
sition affects strongly the couplings of the neutrahnos to gauge bosons and 
sfermions leading to subtantial differences in production and decay rates. 
^ Similarly, the mass matrix for the two charginos can be written in the 
— basis, and the production and decay rates depend on their gaugino 
and Higgsino composition as well. 

As for the squarks and sleptons, their mass eigenstates are in general 
combinations of left-handed and right-handed components. The soft SUSY 
breaking sfermion mass parameters are strongly constrained by the suppres- 
sion of FCNC, requiring that (1) the soft SUSY breaking sfermion mass 
matrix be diagonal and degenerate or (2) the masses of the first and second 
generation sfermions be very large. 

The mixing among left and right-handed sfermions depends on the mass 
of the corresponding SM fermion. So those mixings are negligible for squarks 
and sleptons of the first and second generations but could be significant for 
the fermions of the third generation. In particular, the mass matrix for the 
stops in the tL, tu basis reads 

where rriQ^, At are SUSY breaking parameters and 

Dj^ = M|cos(2/?)(T3j-g/sin2^^) 
Dj^ = M| cos (2/3) g/sin^^v^ 

Unless there is a cancellation of the factor {^At — , the mixing will be 

significant because of the large top mass. This large mixing in turn produces 
a large splitting between the lightest (ti) and the heaviest {12) top squarks. 
This fact makes the lightest stop a good candidate to be one of the lightest 
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superparticles. In this work we shall assume that the stop is the next-to- 
next-to lightest sparticle, details below. 

Another important issue on supersymmetric models is the i?— parity, such 
discrete symmetry is defined as R = (^_i)25+3S+l^ where S is the particle 
spin, B the baryon number, and L the lepton number. It can be checked 
that this quantum number is i? = 1 for the SM particles and i? = ^1 for the 
sparticles, if this symmetry were conserved, a SUSY particle could not decay 
into just SM particles. Specifically, under i?— parity conservation each vertex 
should contain an even number of sparticles. An inmediate consequence is 
that the lightest supersymmetric particle (LSP) must be absolutely stable in 
the case of i?— parity conservation. Thus, from the point of view of collider 
physics the track of such LSP would be missing energy and momentum. 
The violation of i?— parity conduces to the violation of lepton and/or baryon 
number. We shall restrict to the case of i?— parity conservation. 

5.1.2 SUSY breaking schemes 

As we have explained above, supersymmetry must be broken and there is 
a variety of breaking mechanisms. We shall discuss briefiy two schemes: 
supergravity and gauge mediated susy breaking models. 

In supergravity (SUGRA) it is assumed the existence of additional su- 
perfields (the hidden sector) that couples to the MSSM particles by means 
of gravitational interactions. After the SUSY breaking, some components 
of the hidden sector acquire VEV, the soft SUSY breaking terms arise from 
the interaction among the MSSM superfields with the components of the 
hidden sector that acquire a VEV, from which it follows that the strength 
of these SUSY breaking terms is proportional to those VEV divided by the 
Planck scale. Despite the number of SUSY breaking terms is enormous in 
those kind of models, in minimal SUGRA they are reduced considerably be- 
cause the MSSM sparticles couples universally to the hidden sector. This 
universality is extended to the MSSM mass pattern at a scale of the order 
of Mpianck (~ lO^^GeV) or at a scale of the order of GUT (~ lO^^GeV); at 
such a high scale the scalars (Higgs bosons and sfermions) are assumed to 
have a common mass nio, all gauginos (Bino, Wino, and gluino) have a com- 
mon mass mi/2, and all trilinear couplings have a common strength as well 
(Aq). All these values at GUT or Planck scale can be run by Renormalization 
Group Equations (RGE) to obtain their values at EW scale. Evolution of 
the mass pattern by RGE shows that the lightest two neutralinos and the 
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lightest chargino tend to be gaugino like. In addition, squarks are in general 
heavier than sleptons. Finally, the lightest neutralino seems to be a good 
candidate to be the LSP, except in the cases in which mo or mi/2 are very 
small, in those cases the best candidates become the sneutrino and the gluino 
respectively. 

As for the scale of supersymmetry breaking, in SUGRA models it is ex- 
pected to be very high (~ lO^^GeV if we expect the SUSY masses to he at 
the TeV scale) as we shall discuss in next section. Thus, SUGRA are in 
general models with high energy supersymmetry breaking scale. 

Another interesting scheme is the gauge mediated SUSY breaking, in 
which the SUSY breaking terms are generated from gauge interactions. In 
these models the masses of sfermions with the same quantum numbers un- 
der the gauge group i.e. with the same gauge couplings are predicted to 
be degenerate, in this way FCNC are naturally suppressed. As we shall see 
in next section, in Gauge Mediated Spontaneously Breaking (GMSB) mod- 
els the scale of SUSY breaking is much smaller than the Planck or GUT 
scale avoiding corrections at those scales to the degeneracy. The existence 
of heavy messenger superfields is assumed. The breaking of SUSY occurs 
in a hidden sector which also couples to the messenger superfields, thus the 
fermion components of the messenger superfields acquire a commom mass 
M and the scalar components acquire a commom mass Ma/1 ± x where x is 
a dimensionless parameter related to the breaking scale. The gauginos and 
sfermions masses receive radiative contributions from the messenger fields 
that generates a splitting among them and their corresponding SM part- 
ners. Gaugino and scalar masses lie roughly on the same order of magnitude. 
Moreover, after evolving by RGE, sfermions with the same quantum numbers 
continue degenerate if we ignore the effect of Yukawa couplings^, therefore 
the mass hierarchy is directly related to the gauge coupling strength; the 
gauge coupling pattern also determines the gaugino mass hierarchy. In such 
models the gravitino might be very light, and it is a good candidate to be 
the LSP, in that case it plays a crucial role in low energy phenomenology. In 
this work the gravitino will be the LSP. 



^In this work we assume that there are ten degenerate squarks corresponding to the 
five light quarks. Only the top squarks arc; considered to have an sphtting between them, 
owing to the large Yukawa coupling of the top quark. 
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5.2 Introduction 



The Standard Model with a hght Higgs boson provides a very good descrip- 
tion of all experimental data. The consistency of the precision electroweak 
data with the predictions of the Standard Model suggests that, if new physics 
is present at the weak scale, it is most probably weakly interacting and con- 
sistent with the presence of a light Higgs boson in the spectrum. Extensions 
of the Standard Model based on softly broken low energy supersymmetry 
(SUSY) [125] provide the most attractive scenarios of physics beyond the 
Standard Model fulfilling these properties. If the supersymmetry breaking 
masses are < 0{1 TeV), supersymmetry stabilizes the hierarchy between the 
Planck scale Mp and the electroweak scale. Furthermore, the minimal super- 
symmetric extension of the Standard Model (MSSM) significantly improves 
the precision with which the three gauge couplings unify and leads to the 
presence of a light Higgs boson with a mass below about 128 GeV [126]^. 

Perhaps the most intriguing property of supersymmetry is that local su- 
persymmetry naturally leads to the presence of gravity (supcrgravity). In 
the case of local supersymmetry, the Goldstino provides the additional de- 
grees of freedom necessary to make the gravitino a massive particle [127]. In 
the simplest scenarios, the gravitino mass itIq is directly proportional to the 
square of the supersymmetry breaking scale V-^susy^ 

ma ~ Fsvsy/VsMp , (5.2) 

where Mp denotes the Planck mass. 

The relation between the supersymmetry breaking scale V -Fsusy and the 
masses of the supersymmetric partners depends on the specific supersym- 
metry breaking mechanism. In general, the superpartner masses Msusy ^'I'e 
directly proportional to -Fsusy and inversely proportional to the messenger 
scale at which the supersymmetry breaking is communicated to the vis- 
ible sector: 

Msusy Cm^^, (5.3) 

^This upper limit for the lightest CP even Higgs boson could be moved up in some 
non-minimal supersymmetric versions. For instance, Ref. [40] shows that in the super- 
symmetric version of the SU(3)xU(l) gauge model in which the Higgs triplets are included 
in the Icpton supcrficlds, the upper limit for the lightest CP even Higgs boson can be shifted 
up to about 140 GeV. 
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where Cm is the characteristic strength of the couphng between the messen- 
ger sector and the visible one. If the breakdown of supersymmetry is related 
to gravity effects, Mm is naturally of the order of the Planck scale and Cm is 
of order one; hence, for \/Fsusy ~ lO-*^^ GeV, Msusy is naturally at the TeV 
scale. In gauge mediated scenarios (GMSB) [128, 129], instead, the couphngs 
Cm are associated with the Standard Model gauge couplings (times a loop 
suppression factor), so that a FsusY/^m~100 TcV yields masses of the order 
of 100 GeV for the lighter Standard Model superpartners. 

When relevant at low energies, the gravitino interactions with matter 
are well described through the interactions of its spin 1/2 Goldstino compo- 
nent [127]. The Goldstino has derivative couplings with the visible sector 
with a strength proportional to l/FsusY- In scenarios with a high messenger 
scale, of order Mp, Eqs. 5.2 and 5.3 imply that the gravitino has a mass 
of the same order as the other SUSY particles, and its interactions are ex- 
tremely weak. In such scenarios, the gravitino plays no role in the low-energy 
phenomenology. However, in low energy supersymmetry breaking scenarios 
such as GMSB in which the messenger scale is significantly lower than the 
Planck scale, the supersymmetry breaking scale is much smaller. Typical val- 
ues in the GMSB case are Mm ~ 10^ — 10^ GeV, leading to a supersymmetry 
breaking scale V -Fsusy roughly between 10^ and a few times 10^ GeV. The 
gravitino then becomes significantly lighter than the superpartners of the 
quarks, leptons and gauge bosons, and its interaction strength is larger. As 
the lightest supersymmetric particle (LSP), the gravitino must ultimately be 
produced at the end of all superparticle decay chains if i?-parity is conserved 
(for an analysis of the case of i?-parity violation see, Ref. [130]). 

Depending on the strength of the gravitino coupling, the decay length 
of the next-to-lightest supersymmetric particle (NLSP) can be large (so 
that the NLSP is effectively stable from the point of view of collider phe- 
nomenology; this occurs when \/7susy ^ 1000 TeV), intermediate (so that 
the NLSP decays within the detector giving rise to spectacular displaced 
vertex signals; this occurs when 1000 TeV > V -Fsusy ^ 100 TeV), or micro- 
scopic (so that the NLSP decays promptly; this occurs when V Fsusy ^ 100 
TcV) [131, 132, 133]. The decay branching fractions of the Standard Model 
superpartners other than the NLSP into the gravitino are typically negligible. 
However, if the supersymmetry breaking scale is very low, V -Fsusy <S 100 
TeV (corresponding to a gravitino mass <^ 1 eV) , then the gravitino coupling 
strength can become large enough for superpartners other than the NLSP 
to decay directly into final states containing a gravitino [133, 134]. In any 
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case, the SUSY-breaking scale must be larger than the mass of the heaviest 
superparticle; an approximate lower bound of V i^suSY — 1 TeV corresponds 
to a gravitino mass of about 10^'^ eV. 

In many models, the lightest amongst the supersymmetric partners of 
the Standard Model particles is a neutralino, Xv The partial width for Xi 
decaying into the gravitino and an arbitrary SM particle X is given by 

r«;^XG).K.JV.^(^j (l--fj . (5.4, 

where Kx is a projection factor equal to the square of the component in the 
NLSP of the superpartner X, and Nx is the number of degrees of freedom 
of X. If X is a photon, for instance, Nx — 2 and 

Kx = l-ZVii cos Ow + Ni2 sin 6*^1^1^ , (5.5) 

where Nij is the mixing matrix connecting the neutralino mass eigenstates 
to the weak eigenstates in the basis B, W, Hi, H2. 

If the neutralino has a significant photino component, it will lead to ob- 
servable decays into photon and gravitino. Since the heavier supersymmetric 
particles decay into the NLSP, which subsequently decays into photon and 
gravitino, supersymmetric particle production will be characterized by events 
containing photons and missing energy. This is in contrast to supergravity 
scenarios, where, unless very specific conditions are fulfilled [135, 136], pho- 
tons do not represent a characteristic signature. The presence of two en- 
ergetic photons plus missing transverse energy provides a distinctive SUSY 
signature with very little Standard Model background. 

It might be argued that a Xit decaying with a large branching ratio into 
photons, would be severely constrained by LEP data. However, such bounds 
are extremely model-dependent. For example, there exists no tree-level cou- 
pling between a photon (or Z) and either two binos or two neutral winos. 
Since the bino is associated with the smallest of all gauge interactions and, in 
addition, its mass is more strongly renormalized downward at smaller scales 
compared to the wino mass, in many models the lightest neutralino has a 
significant bino component. Therefore, if the NLSP is approximately a pure 
bino, the bounds on its mass depend strongly on the selectron mass (since, 
at LEP, pair-production of binos or neutral winos could occur through t- 
channel exchange of selectrons). For selectron masses below 200 GeV, the 
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present bound on such a neutralino is approximately 90 GeV [137] (the bound 
weakens with increasing selectron mass). As emphasized before, once pro- 
duced, such a neutrahno would decay via Xi ~^ iG. If x5 is heavy enough, 
then the decays x? — > ZG and — > hPG are also allowed; however, the decay 
widths into these final states are kinematically suppressed compared to the 

final state and will be important only if either the photino component of 
X? is small or if x? itself is significantly heavier than Z and [133]. 

In most SUSY models, it is natural for the lighter top squark, ti, to be 
light compared to the other squarks. In general, due to the large top Yukawa 
coupling, there is a large mixing between the weak eigenstates and tn, 
which leads to a large splitting between the two stop mass eigenstates. In 
addition, even if all squarks have a common mass at the messenger scale, the 
large top Yukawa coupling typically results in the stop masses being driven 
(under renormalization group evolution) to smaller values at the weak scale. 
An extra motivation to consider light third generation squarks comes from 
the fact that light stops, with masses of about or smaller than the top quark 
mass, are demanded for the realization of the mechanism of electroweak 
baryogenesis within the context of the MSSM [138]. 

In this chapter, we examine, in detail, the production and decay of top 
squarks at Run II of the Tevatron collider in low-energy SUSY breaking 
scenarios wherein the lightest neutralino is the NLSP and decays promptly 
into 7(7. We also investigate the production and decay of the other squarks, 
and provide an estimate of the reach of Run II of the Tevatron in the heavy 
gluino limit. We work in the context of a general SUSY model in which the 
SUSY particle masses are not constrained by the relations predicted in the 
minimal GMSB models. We assume throughout that the gravitino coupling 
is strong enough (or, equivalently, that the scale of SUSY breaking is low 
enough) that the NLSP decays promptly. This implies an upper bound on 
the supersymmetry breaking scale of a few tens to a few hundred TeV [132, 
133, 131], depending on the mass of the NLSR Our analysis can be extended 
to higher supersymmetry breaking scales for which the NLSP has a finite 
decay length, although in this case some signal will be lost on account of the 
NLSP decaying outside the detector. At least 50% of the diphoton signal 
cross section remains though for NLSP decay lengths cr < 40 cm [131]; 
this corresponds to a supersymmetry breaking scale below a few hundred 
to about a thousand TeV, depending on the mass of the NLSP. Moreover, 
the displaced vertex associated with a finite decay length could be a very 
good additional discriminator for the signal. Thus, in totality, our choice is 
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certainly not an overly optimistic one. 

This chapter is organized as follows. In Sec. 5.3, we outline models of 
low-energy supersymmetry breaking wherein the lighter stop is the lightest 
sfermion and, moreover, is lighter than the charginos as well as the gluino. In 
the following section, we review the stop pair production cross section at the 
Tevatron. In Sec. 5.5, we summarize previous studies of stop production and 
decay at Run II of the Tevatron. This is followed, in Sec. 5.6, by a discussion 
of the SUSY parameter space and the relative partial widths of the various 
stop decay modes. In Sec. 5.7, we describe the signal for each of the stop 
decay modes considered. We describe the backgrounds and the cuts used to 
separate signal from background in each case, and give signal cross sections 
after cuts. This gives the reach at Run II. Wc also note the possibility that 
stops can be produced in the decays of top quarks. In Sec. 5.8 we consider the 
production and decay of 10 degenerate squarks that are the super symmetric 
partners of the five light quarks. Finally, we summarize our results (along 
with the results concerning the previous chapters of this Ph. D. thesis) in 
chapter 6. 



5.3 Light Stop in Low Energy Supersymme- 
try Breaking Models 

As discussed above, the mass of the gravitino as well as its interaction 
strength are governed by the supersymmetry breaking scale V-^susy- Low 
energy supersymmetry breaking models are defined as those obtained for 
low values of ^ -Fsusy (~10® GeV) and hence resulting in a gravitino lighter 
than a few keV. Apart from evading cosmological problems [139], a further 
striking consequence of such models is that sparticle decays into gravitinos 
may occur at scales that may be of interest for collider phenomenology. As 
is apparent from Eq. (5.3), in such models the messenger mass scale Myn 
should be smaller than 10^ GeV. 

In this chapter, wc are interested in the presence of light stops in low en- 
ergy supersymmetry breaking scenarios. The motivation is simple: assuming 
that the supersymmetry breaking masses are flavor independent, that is the 
left- and right-handed squark masses of the three generations are the same 
at the messenger scale, the lighter stop turns out to be the lightest of all 
the squarks. The reasons are twofold. On the one hand, there are renormal- 
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ization group effects induced by top-quark Yukawa interactions tfiat tend to 
reduce tlie stop mass scale compared to tlie otlier squark masses. On tlie 
otlier, tliere are non-trivial mixing effects that tend to push the lightest stop 
mass down compared to the overall left- and right-squark masses. For similar 
reasons, it is also natural to assume that the lightest stop will be lighter than 
the gluino. 

There is a further motivation behind our analysis, namely that of baryo- 
genesis. A crucial requirement for electroweak baryogenesis scenarios within 
the Minimal Supersymmetric Extension of the Standard Model is the pres- 
ence of a light stop with mass of the order of, or smaller than, the top quark 
mass. Since this scenario does not depend on the nature of supersymme- 
try breaking, it is very important to develop strategies to look for light top 
squarks in all their possible decay modes, in particular in those related to 
the possibility of low energy supersymmetry breaking. 

In general, light stops can induce large corrections to the precision elec- 
troweak parameter Ap [140, 141, 142], unless they are mainly right-handed 
or there is some correlation between the masses and mixing angles in the stop 
and sbottom sector [143]. The most natural way of suppressing potentially 
large contributions to the rho-paramctcr is to assume that there is a large 
hierarchy between the supersymmetry breaking mass parameters for the left- 
and right-handed stops. The simplest and most efficient way of doing this is 
to assume that sparticles which are charged under the weak gauge interac- 
tions acquire large supersymmetry breaking masses. Observe that, under this 
assumption, the charged wino will also be naturally heavier than the lightest 
stop and, for simplicity, we will assume that both charginos are heavier than 
the particle under study. We will further assume that the superpartner of 
the U{1)y gauge boson, the so-called bino, is the lightest standard model 
superpartner. 

The above mentioned properties are naturally obtained in simple exten- 
sions of the minimal gauge mediated model. Indeed, let us assume that there 
are copies of messengers, which transform as complete representations 
of SU (5) and hence do not spoil the unification relations. Under the Standard 
Model gauge group, some of these fields will then transform as left-handed 
lepton multiplets (W,, i — 1 . . . N) and their mirror partners (W^), while the 
others would transform as a right-handed down quark multiplet (Cj) and its 
mirror partner {Ci). We shall assume that < 4 in order to keep the gauge 
couplings weak up to the grand unification scale [144]. Let us further assume 
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that Wi and Cj couple to two different singlet fields ^'2,3, with 

{S,} = S, + e'F, , (5.6) 

and, for simplicity, assume that all couplings are of order one and that Fj ^ 
Sj. We will further assume that Si S2 — S, where S characterizes the 
messenger scale. 

In this simple case, the masses of the gauginos at the messenger scale will 
be given by 

Ms = iV^Aa, 
M2 = N^A2 , 

where A3 ^ and A2 ~ F2/S2 [144, 145]. Now, it is easy to see that for 

A2 > 3A3, the weak gaugino can have a mass similar to that of the gluino (or 
be even heavier), while the bino is still much lighter than the wino (about 
three times lighter than the wino at the messenger scale). 

The squared scalar masses at the messenger scale are also proportional 
to the number of messengers N 



m^Q = 2N 



a=2,3 ^ 



, (5. 



where, for a particle transforming in the fundamental representation of SU (n), 
eg = (n^ - l)/(2n), while = 3/5{Q - Ta)^. The above quoted masses 
should be renormalized to the weak scale. The details of this procedure in 
a general gauge mediated model have been given by one of us in Ref . [145] . 
We shall not repeat these expressions here. For our purpose, it suffices to 
stress some important details. 

As wc mentioned before, after renormalization and mixing effects are 
added, the right-handed stop is the lightest squark in the spectrum. The large 
value of the left-handed stop mass increases the negative Yukawa dependent 
radiative corrections to the right-handed stop mass. Therefore, for A2 larger 
than a few times A3 and > 1, the lightest stop is lighter than the gluino 
and the wino. Even for A^ = 1 this tends to be true, for not too small values 
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of the messenger scale (S'a^lO^ GeV). Notice that, for A2 larger by an order 
of magnitude than A3 and/or large values of A^, the right-stop (or, in extreme 
cases the right-sbottom) becomes the lightest standard model superpartner. 
We shall thus concentrate on cases where A2, though larger than A3, is still 
of the same order. 

The sleptons do not play a relevant role in stop decays as long as the 
charginos are heavier than the lightest stop. One has only to ensure that the 
lightest slepton does not become lighter than the bino. The hierarchy of the 
right-handed slepton mass and the bino mass is approximately given by 

where ri is the appropriate renormalization group factor relating the masses 
at the messenger scale to the masses at the weak scale. This factor is about 
1.5-1.7 [144] in the region of interest for this work and therefore we are 
led to conclude that, as long as < 3, the bino is the next-to-lightest 
supersymmetric particle. Observe that this constraint on A^ is somewhat 
weaker than in the minimal model, due to the appearence of the factor 3/5 
in front of the dominant A2 contribution. 

In all of the above, we have not discussed the process of radiative elec- 
troweak symmetry breaking and the determination of the ^ parameter. This 
is due to the fact that in the minimal gauge mediated models of the kind 
we described above, the generation of a proper value of the parameter Bfj, 
requires the presence of new physics that necessarily modifies the value of 
the Higgs mass parameters [146]. It is, therefore, justified to treat /i as an 
independent parameter in such models. 

To summarize, we have seen that a mild modification of the simplest 
gauge mediated models leads to a model where the stop is lighter than all the 
other squarks and also lighter than the weak and strong gauginos, while the 
bino remains the next-to-lightest supersymmetric particle. The unification 
relations are preserved as long as the number of messenger families is less than 
or equal to 3. The only condition for this to happen is that the effective scale 
A2 is a few times larger than A3, but still of the same order of magnitude. 
The relation between the stop mass and the bino mass will be governed by 
the hierarchy between A2 and A3. This modification of the minimal gauge 
mediated models is well justified in order to get consistency of a light stop 
with electroweak precision observables. Some of the sleptons might be lighter 
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than the stop, but, as long as the charginos remain heavy, they have no 
impact on the stop phenomenology. 

Let us stress that the above should be considered only as a simple exam- 
ple in which a light stop, mainly right-handed, can appear in the spectrum 
in low energy supersymmetry breaking models. As the nature of supersym- 
metry breaking is unknown, so is the exact spectrum. We only make the 
simplifying assumptions that the left- and right-handed sfermions receive an 
approximately common mass at the messenger scale and that the gaugino 
masses are of the same order of magnitude as the squark masses. The truly 
defining feature of our assumption is that the wino and the gluino are heavier 
than the squarks and that the bino is the ncxt-to-lightest supersymmetric 
particle. Although the left-handed squarks tend to be heavier than the right- 
handed ones, the exact spectrum obtained in the simple model detailed above 
does not differ markedly from the simplified spectrum that we choose to work 
with. In this kind of models, apart from the somewhat lighter stops, there 
will be ten squarks approximately degenerate in mass, which will predom- 
inantly decay into a quark and a bino, which will subsequently decay into 
photons (or Z bosons) and missing energy. 

5.4 Top squark production at Tevatron Run II 

Top squarks are produced at hadron colliders overwhelmingly via the strong 
interaction, so that the tree-level cross sections are model independent and 
depend only on the stop mass. The production modes of the lightest stop, 
ti, at the Tevatron are qq — > t\i\ and gg t\t\. The cross sections for these 
processes are well known at leading order (LO) [147, 148], and the next-to- 
leading order (NLO) QCD and SUSY-QCD corrections have been computed 
[149] and significantly reduce the renormalization scale dependence. The 
NLO cross section is implemented numerically in PROSPINO [150, 149]. 

We generate stop events using the LO cross section evaluated at the 
scale /i = m^, improved by the NLO K-factor ^ obtained from PROSPINO 
[150, 149] (see Fig. 5.1). The K-factor varies between 1 and 1.5 for mi 
decreasing from 450 to 100 GeV. We use the CTEQ5 parton distribution 
functions [151]. We assume that the gluino and the other squarks are heavy 

^Although gluon radiation at NLO leads to a small shift in the stop distribution 
to lower pt values [150, 149], we do not expect this shift to affect our analysis in any 
significant way. 
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enough that they do not affect the NLO cross section. This is aheady the 
case for gluino and squark masses above about 200 GeV [149]. 




100 150 200 250 300 350 400 450 
Stop mass (GeV) 

Figure 5.1: LO (dashed) and NLO (sohd) cross sections for stop pair produc- 
tion in pp colhsions at Tevatron Run I (1.8 TeV) and Run II (2.0 TeV), from 
PROSPINO [150, 149]. Cross sections are evaluated at the scale /i = m^. 

Top squarks can also be produced via cascade decays of heavier super- 
symmetric particles, with a highly model dependent rate. To be conservative, 
we assume that the masses of the heavier supersymmetric particles are large 
enough that their production rate at Tevatron energies can be neglected. 

5.5 Previous studies of stops at Run II 

A number of previous studies have considered the prospects for stop discovery 
at Run II of the Tevatron, which we summarize here. In general, the most 
detailed SUSY studies have been done in the context of supergravity; in this 
case SUSY is broken at the Planck scale so that the gravitino plays no role 
in the collider phenomenology. Then the lightest neutralino is the LSP and 
ends all superparticle decay chains. The reach of the Tevatron for a number 
of stop decay modes has been analyzed in Refs. [152, 153]. The signal 
depends on the decay chain, which in turn depends on the relative masses 



108 



of various SUSY particles. For sufficiently heavy stops, the decay t — > 
will dominate. This channel is of limited use at Run 11 because the stop pair 
production cross section falls rapidly with increasing stop mass, and this 
channel requires > rrit + m^o , which is quite heavy for the Tevatron in 
the case of minimal supergravity.^ For lighter stops, if a chargino is lighter 
than the stop then t — > bxt tends to dominate (followed by the decay of the 
chargino). The details of the signal depend on the Higgsino content of Xi- If 
Xi has a mass larger than rrif—mi,, the previous decay does not occur and the 
three-body decay i — > bW~^Xi dominates; this decay proceeds through the 
exchange of a virtual top quark, chargino, or bottom squark. If the stop is 
too light to decay into an on-shell W boson and Xi, then the flavor- changing 
decay t cXi tends to dominate. Finally, if a sneutrino or slepton is light, 
then t M^Ui or t bi^ug, respectively, will occur (followed by the decays 
of the slepton or the sneutrino, if it is not the LSP). At Run 11 with 2 (20) 
fb~^ of total integrated luminosity, in the context of minimal supergravity 
one can probe stop masses up to 160 (200) GeV in the case of the flavor 
changing decay, while stop masses as high as 185 (260) GeV can be probed 
if the stop decays into a bottom quark and a chargino [153]. A similar reach 
holds in the case of a light sneutrino [153] and in the case of large tan /? when 
the stop can decay into bruxi [154]. 

One can also search for stops in the decay products of other SUSY parti- 
cles [153]. In top decays, the process t ^ ig is already excluded because of the 
existing lower bound on the gluino mass ^. Other possibilities are x~ ~^ bt* 
and g ^ tt*. Finally, the decays b — >■ tW~ and tH~ have to compete with the 
preferred decay, b bxi, and so may have small branching ratios (depending 
on the masses of f, 6, Xi ^iid The signals for these processes at the 

Tevatron Run II have been considered in minimal supergravity in Ref . [153] . 

If R-parity violation is allowed, then single stop production can occur via 
the fusion of two quarks. Single stop production is kinematically favored 
compared to stop pair production and offers the opportunity to measure 
R-parity violating couplings. This has been considered for the Tevatron in 
Ref. [157], which showed that the stop could be discovered at Run II for 
masses below about 400 GeV provided that the R-parity violating coupling 

^In low-energy SUSY breaking scenarios, however, the mass range rrif > mt + m^o is 
interesting at Tevatron energies because the distinctive signal allows backgrounds to be 
reduced to a very low level, as we will show. 

■''The bound on the gluino mass is, however, model dependent. Under certain conditions, 
an allowed window exists for gluino masses below the gauge boson masses [155, 156]. 



109 



A" > 0.02 — 0.1 and that the stop decays via ii — > bxt (followed by xt ~^ 

Relatively few studies have been done in the context of low-energy SUSY 
breaking with a gravitino LSP. A study of GMSB signals performed as part 
of the Tevatron Run II workshop [131] considered the decays of various SUSY 
particles as the NLSP. As discussed before, the NLSP in such models will 
decay directly to the gravitino and Standard Model particles. If the stop 
is the NLSP in such a model, then it will decay via i t^*^G — > bW~^G 
(for rrif > nib + rnw)- Note that because G is typically very light in such 
models, mi > mw + iTib is sufficient for this decay to proceed with an on-shell 
W boson. Rcf. [153] found sensitivity at Run II to this decay mode for stop 
masses up to 180 GeV with 4 fb~^. This stop decay looks very much like a top 
quark decay; nevertheless, even for stop masses near rrit, such stop decays 
can be separated from top quark decays at the Tevatron using kinematic 
correlations among the decay products [158]. 

Finally, Ref. [133] considered general GMSB signals at the Tevatron of 
the form 77 + X. The authors of Ref. [133] provide an analysis of the 
possible bounds on the stop mass coming from the Run 1 Tevatron data. 
They analyze the stop decay mode into a charm quark and a neutralino, and 
also possible three body decays, by scanning over a sample of models. They 
conchidc that stop masses smaller than 140 GeV can be excluded already by 
the Run 1 Tevatron data within low energy supersymmetry breaking models 
independent of the stop decay mode, assuming that m^o > 70 GeV. 

5.6 Top squark decay branching ratios 

The decay properties of the lighter stop depend on the supersymmetric parti- 
cle spectrum. Of particular relevance are the mass splittings between the stop 
and the lightest chargino, neutralino and bottom squark. In our analysis, we 
assume that the charginos and bottom squarks arc heavier than the lighter 
stop, so that the on-shell decays t bW and t x^b are kincmatically 
forbidden. Then the details of the stop decay depend on the mass splitting 
between the stop and the lightest neutralino. If mf < mw + mb + two 
decay modes are kincmatically accessible: 

1. the flavor-changing (FC) two-body decay t ex?- This two-body 
decay proceeds through a flavor- changing loop involving , or 
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X+ exchange or through a tree-level diagram with a,t — c mixing mass 
insertion; 

2. the four-body decay via a virtual W boson, t — > W'^*bxi — > jjbxi or 
Mx\ [159]. 



The branching ratio of the stop decay into charm and neutralino strongly 
depends on the details of the supersymmetry breaking mechanism. In models 
with no flavor violation at the messenger scale, the whole effect is induced 
by loop effects and receives a logarithmic enhancement which becomes more 
relevant for larger values of the messenger mass. In the minimal supcrgravity 
case, the two-body FC decay branching ratio tends to be the dominant one. 
In the case of low energy supersymmetry breaking, this is not necessarily 
the case. Since the analysis of the four-body decay process is very similar 
to the three-body decay described below for larger mass splittings between 
the stop and the lighter neutralino, here we shall analyze only the case in 
which the two-body FC decay t — > cx5 is the dominant one whenever mi < 
mw -\-mb-\- m^o . 

For larger mass splittings, so that mw + m^ + m^o < m^ < m^ + m^o, 

the three-body decay i — > W'^bxi becomes accessible, with Xi iG. This 

stop decay proceeds through a virtual top quark, virtual charginos, or virtual 
sbottoms. Quite generally, this 3-body decay will dominate over the 2-body 
FC decay in this region of phase space. 

For still heavier stops, m^ > mt+m^a, the two-body tree-level decay mode 
t — > txi becomes kinematically accessible and will dominate. Although the 
3-body and 2-body FC decays are still present, their branching ratios are 
strongly suppressed. 

Let us emphasize that, since the bino is an admixture of the zino and the 
photino, a pure bino neutralino can decay into either 76* or ZG (see Eq. 5.4). 
If the lightest neutralino is a mixture of bino and wino components, then the 
relative zino and photino components can be varied arbitrarily, leading to 
a change in the relative branching ratios to 7G and ZG. If the hghtest 
neutralino contains a Higgsino component, then the decay to h'^G is also 
allowed. We show in Fig. 5.2 the branching ratio of the lightest neutralino 
into jG as a function of its mass and Higgsino content. 
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Figure 5.2: Branching ratio for the decay of the hghtest neutrahno into a 
photon and a gravitino as a function of the neutrahno mass. Shown are 
the branching ratio if the neutrahno is a pure bino (sohd hne) and for 20% 
and 50% Higgsino admixtures (long and short dashed hues, respectively) 
assuming that m/jo = 120 GeV and that the other MSSM Higgs bosons are 
very heavy. For the Higgsino admixture in the lightest neutralino, we choose 
the H1-H2 mixing so that the field content is aligned with that of and the 
longitudinal component of the Z boson in order to minimize the branching 
ratio to photons. 



5.7 Top squark signals in low-energy SUSY 
breaking 



As discussed in the previous section, the decay properties of the lighter stop 
depend primarily on the mass splitting between the stop and the lightest 
neutralino. In this section, we proceed with the phenomenological analysis 
of the signatures of top squark production associated with the different decay 
modes. In Sec. 5.7.1 we shall analyze the signatures associated with the two- 
body FC decay, which after the neutralino decay leads to t — > cjG. In 
Sec. 5.7.2 we shall analyze the signatures associated with the three-body 
decay, which after neutralino decay leads to t — > hW^'~^G. The two-body 
decay t txi, which typically dominates for > + m^o, leads to the 
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same final state as the three-body decay and therefore the discussion of this 
case will be included in Sec. 5.7.2. Finally, in Sec. 5.7.3 we consider the 
possibility that stops are produced in the decays of top quarks. 

5.7.1 Two-body FC decay: i ^ cjG 

With the stop undergoing the aforementioned decay, the final state consists 
of a pair each of charm jets, photons and (invisible) gravitinos. As we will 
show below, the backgrounds to this process are small enough that we need 
not require charm tagging. Thus the signal consists of 

2(jets) + 27+ ^T- 

The selection criteria we adopt are: 

1. each event must contain two jets and two photons, each of which should 
have a minimum transverse momentum (px > 20 GeV) and be con- 
tained in the pseudorapidity range —2.5 < rj < 2.5; 

2. the jets and the photons should be well separated from each other; 
namely, 

5Rjj > 0.7, 6R^^ > 0.3, 5Rj-y > 0.5 

where SR^ = 5rf + 50^, with 57] {6(f)) denoting the difference in pseu- 
dorapidity (azimuthal angle) of the two entities under consideration; 

3. the invariant mass of the two jets should be sufficiently far away from 
the W- and Z-masses: 

TJijj ^ (75 GeV, 95 GeV); 

4. each event should be associated with a minimum missing transverse 
momentum (j/t > 30 GeV). 

The photons and jets in signal events tend to be very central; in partic- 
ular, reducing the pseudorapidity cut for the two photons to —2.0 < 1] < 2.0 
would reduce the signal by less than about 3%. (This change would re- 
duce the background by a somewhat larger fraction.) Apart from ensuring 
observability, these selection criteria also serve to eliminate most of the back- 
grounds, which are listed in Table 5.1. 

A primary source of background is the SM production of jj^'^ViVi where 
the jets could have arisen from either quarks or gluons in the final state 
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Background 


Cross section after cuts 


after 7 ID 


33TlZ, Z ^vv 


~ 0.2 fb 


~ 0.13 fb 


ii^vv + 7 radiation 


~ 0.002 fb 


~ 0.001 fb 


6677, CC77 


< 0.1 fb 


~ 0.06 fb 


ihi 


~ 0.2 fb 


~ 0.13 fb 


Backgrounds with fake photons: 






jj(ee 77) 


~ 5 X 10-4 fb 


~ 5 X 10-4 fb 


33l{3 7) 


~ 0.8 fb 


~ 0.8 fb 


33{33 11) 


~ 0.8 fb 


~ 0.8 fb 


Total 


~ 2 fb 


~ 2 fb 



Table 5.1: Backgrounds to tt* jj'yy ^t- The photon identification effi- 
ciency is taken to be = 0.8 for each real photon. See text for details. 

of partonic subprocesses. A full diagrammatic calculation would be very 
computer-intensive and is beyond the scope of this work. Instead, we consider 
the subprocesses that are expected to contribute the bulk of this particular 
background, namely pp ^ 2j + 2'y + Z + X with the Z subsequently decaying 
into neutrinos. These processes are quite tractable and were calculated with 
the aid of the helicity amplitude program Madgraph [160]. On imposition 
of the abovementioned set of cuts, this background at the Run II falls to 
below ~ 0.2 fb. 

An independent estimate of the jj^^i'iPi background may be obtained 
through the consideration of the single-photon variant, namely jj'ji'ii'i pro- 
duction, a process that Madgraph can handle. After imposing the same 
kinematic cuts (other than requiring only one photon) as above, this process 
leads to a cross section of roughly 0.2 fb. Since the emission of a second hard 
photon should cost us a further power of aem, the electromagnetic coupling 
constant, this background falls to innocuous levels. We include both this 
estimate and the one based on Z production from the previous paragraph in 
Table 5.1. Though a naive addition of both runs the danger of overcounting, 
this is hardly of any importance given the overwhelming dominance of one. 

A second source of background is 6677 (0077), with missing transverse 
energy coming from the semileptonic decay of one or both of the b (c) mesons. 
In this case, though, the neutrinos tend to be soft due to the smallness of the 
b and c masses. Consequently, the cut on ]/t serves to eliminate most of this 
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background, leaving behind less than 0.1 fb. This could be further reduced (to 
< 0.001 fb) by vetoing events with leptons in association with jets. However, 
such a lepton veto would significantly impact the selection efficiency of our 
signal, which contains cc, thereby reducing our overall sensitivity in this 
channel. 

A third source of background is j.777 production, in which the jet (light 
quark or gluon) and/or photon energies are mismeasured leading to a fake 
^T- To simulate the effect of experimental resolution, we use a (very pes- 
simistic) Gaussian smearing: 6Ej/Ej = 0.1 + 0.6/ ^/Ej{GeV) for the jets and 
SE^/E^ = 0.05 + 0.3/^/K^{GeV) for the photons.^ While the production 
cross section is much higher than that of any of the other backgrounds con- 
sidered, the ensuing missing momentum tends to be small; in particular, our 
cut on reduces this background by almost 99%.^ On imposition of our 
cuts, this fake background is reduced to ~ 0.2 fb.^ 

Finally, we consider the instrumental backgrounds from electrons or jets 
misidentified as photons. Based on Run I analyses [162] of electron pair 
production and taking the probability for an electron to fake a photon to be 
about 0.4%, we estimate the background due to electrons faking photons to 
be of order 5 x 10^"^ fb. More important is the background in which a jet 
fakes a photon. Based on a Run I analysis [161] and taking the probability 
for a jet to fake a photon to be about 0.1%, we estimate the background due 
to jjjj in which one of the jets fakes a second photon to be about 0.8 fb. 
Similarly, we estimate that the background due to jjjj in which two of the 
jets fake photons is somewhat smaller; to be conservative we take it to be of 
the same order, i.e., 0.8 fb. 

Having established that the backgrounds are small, let us now turn to 
the signal cross section that survives the cuts. In Fig. 5.3, we present these 
as contours in the mf-m^o plane. We assume that the branching ratio of 
i — > cxI dominates in the region of parameter space that we consider here. 

^We have also performed similar smearing for the other background channels as well 
as for the signal. However, there it hardly is of any importance as far as the estimation of 
the total cross section is concerned. 

'^This reduction factor is in rough agreement with that found for 7 + j events in CDF 
Run I data in Ref. [161]. 

^We note that in the squark searches in supergravity scenarios the signal is jj ^t; 
in this case a similar background due to dijct production with fake is present. This 
background is larger by two powers of a than the ^,777 background, yet it can still be 
reduced to an acceptable level by a relatively hard cut onpr (see, e.g., Ref. [153]). 
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Figure 5.3: Cross sections in fb for stop pair production in Run II with 
i — > C'jG, after cuts. No efficiencies have yet been applied. The black area is 
excluded by non-observation of jj'j'j events in Run I. 

The branching ratio of Xi ^ iG is taken from Fig. 5.2 assuming that x? is 
a pure bino.^ All our plots are made before detector efficiencies are applied. 
With a real detector, each photon is identified with about = 80% efficiency; 
thus the cross sections shown in Fig. 5.3 must be multiplied by = 0.64 
in order to obtain numbers of events. While the production cross sections 
are independent of the neutralino mass, the decay kinematics have a strong 
dependence on m^o. For a given m^, a small mass splitting between mi 
and m^o would imply a soft charm jet, which would often fail to satisfy our 



^If Xi is not a pure bino, its branching ratio to will typically be reduced somewhat 
when its mass is large (see Fig. 5.2). This will lead to a reduction of the signal cross 
section at large m^o^ by typically a few tens of percent. 

^°The diphoton trigger efficiency is close to 100%, so we neglect it here. 
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selection criteria. This causes the gap between the cross section contours 
and the upper edge of the parameter space band that we are exploring in 
Fig. 5.3. This is further compounded at large m^o by the fact that a large 

neutralino mass typically implies a smaller Xi ~^ iG branching fraction (see 
Fig.5.2). On the other hand, a small m^o results in reduced momenta for the 
gravitino and the photon, once again resulting in a loss of signal; however, 
this is important only for m^o < 70 GeV, which is not relevant in this search 
channel. 

The signal cross sections are fairly substantial. In particular, the dark 
area in Fig. 5.3 corresponds to a signal cross section of 50 fb or larger at Run 
I of the Tevatron. Taking into account the identification efficiency of 64% for 
the two photons, such a cross section would have yielded 3 signal events in 
the 100 pb^^ collected in Run 1 over a background of much less than 1 event. 
Run 1 can thus exclude this region at 95% confidence level. In particular, 
we conclude that Run I data excludes stop masses up to about 200 GeV, for 
large enough mass splitting between the stop and the neutralino. When the 
stop-neutralino mass splitting is small {i.e., less than about 10 — 20 GeV), 
the charm quark jets become too soft and the signal efficiency decreases 
dramatically. For comparison, a D0 search [163] for inclusive pp — X2 + ^ 
with X2 ~^ 7X1 ill the context of minimal supergravity yields a limit on the 
production cross section of about 1 pb for parent squark masses of order 150- 
200 GeV. Interpreting this in terms of stop pair production with — * 7%° 
reidentified as Xi ~^ iG increases the signal efficiency by a factor of ~ 2.7 
because every SUSY event now contains two photons [163]; the D0 analysis 
then yields a stop mass bound of about 180 GeV, in rough agreement with 
our result. In addition, Ref. [133] projected a Run I exclusion of stops 
in this decay channel for masses below about 160-170 GeV, again in rough 
agreement with our result. 

To claim a discovery at the 5(T level, one must observe a large enough 
number of events that the probability for the background to fluctuate up to 
that level is less than 5.7x 10~^. Because the number of expected background 
events in this analysis is small, we use Poisson statistics to find the number 
of signal events required for a 5(T discovery. Taking the total background 
cross section to be 2 fb from Table 5.1, we show in Table 5.2 the expected 



i^Note, however, that the non-iicgligible mass of the LSP of about 35 GeV assumed in 
Rcf. [163] leads to kinematics tliat differ significantly from those in our analysis, in which 
the gravitino is essentially massless. 
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maximum stop discovery mass reach at Tevatron Run II for various amounts 
of integrated luminosity. In particular, with 4 fb~^ a stop discovery can 
be expected in this channel if < 285 GeV, with S/ B of more than 2/\}^ 
Including the effects of mixing in the composition of the lightest neutrahno, 
a 50% reduction in the — > branching ratio compared to the pure bino 
case would reduce the stop mass reach by only about 20 GeV. For such a 
reduction to occur in the relevant neutralino mass range of about 200 — 250 
GeV, the lightest neutralino would have to be less than half bino. 





B 


S for a 5(7 discovery 




Maximum stop mass reach 


2 fb-i 


4 


14 


7.0 fb 


265 GeV 


4 fb-i 


8 


18 


4.5 fb 


285 GeV 


15 fb-i 


30 


31 


2.1 fb 


310 GeV 


30 fb-i 


60 


42 


1.4 fb 


325 GeV 



Table 5.2: Number of signal events (S) required for a 5(7 stop discovery at 
Tevatron Run II in the j j77 -^t channel and the corresponding signal cross 
section after cuts and efficiencies and maximum stop mass reach. We assume 
a photon identification efficiency of e-y = 0.80. The number of background 
events (B) is based on a background cross section of 2 fb from Table 5.1. 



5.7.2 Three-body decay: i hW+^G 

For a large enough splitting between the stop and neutralino masses, the 
signature of stop pair production would be^^ jjWW^^ ^t- In this analysis, 
we consider only the dominant hadronic decay mode of both W bosons. 

"'^^The cross sections and numbers of events required for discovery are quoted in terms 
of integrated luminosities at a single detector. If data from the CDF and D0 detectors 
are combined, the integrated luminosity of the machine is effectively doubled. 

^^For comparison, in the case of minimal supergravity a reach of < 180 GeV can 
be expected in the i ex? channel with 4 fb~^ at Tevatron Run II; the same reach is 
obtained in low-energy SUSY breaking scenarios in which the stop is the NLSP rather 
than the neutralino [153]. In both of these cases, the signal consists of jj ^t, with no 
photons in the final state. 

^^The backgrounds are again small enough after cuts in this channel that we do not need 
to tag the b quarks. In the case of a discovciy, one could imagine tagging the b quarks 
and reconstructing the W bosons in order to help identify the discovered particle. 
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This decay mode of the stop proceeds via three Feynman diagrams, in- 
volving an intermediate off-shell top quark, chargino or sbottom. We use the 
full decay matrix elements as given in Ref. [164]. Although this introduces 
several additional parameters into the analysis, a few simplifying assump- 
tions may be made without becoming too model dependent. For example, 
assuming that the lightest stop is predominantly the superpartner of the 
right-handed top-quark (t/?), eliminates the sbottom exchange diagram alto- 
gether. Even if the stop contains a mixture of and t^, under our assump- 
tion that the hghter stop is the next-to-hghtest Standard Model superpartner, 
the sbottom exchange diagram will be suppressed by the necessarily larger 
sbottom mass. As for the chargino exchange, the wino component does not 
contribute for a decay. Thus the chargino contribution is dominated by 
its Higgsino component. Furthermore, if we concentrate on scenarios with 
large values of the supersymmetric mass parameter // and of the wino mass 
parameter (in which case the charginos are heavy and the neutralino is al- 
most a pure bino) , the chargino exchange contribution is also suppressed and 
the dominant decay mode is via the diagram involving an off-shell top quark. 
To simplify our numerical calculations, we have then assumed that only this 
diagram contributes to the stop decay matrix element. We have checked for 
a few representative points though, that the inclusion of the sbottom and 
chargino diagrams does not significantly change the signal efficiency after 
cuts as long as we require that > mi. 

As the W bosons themselves decay, it might be argued that their polar- 
ization information needs to be retained. However, since we do not consider 
angular correlations between the decay products, this is not a crucial issue; 
the loss of such information at intermediate steps in the decay does not lead 
to a significant change in the signal efficiency after cuts. This is particularly 
true for the hadronic decay modes of the W , for which the profusion of jets 
frequently leads to jet overlap, thereby obscuring detailed angular correla- 
tions. It is thus safe to make the approximation of neglecting the polarization 
of the W bosons in their decay distributions, and we do so in our analysis. 

Before we discuss the signal profile and the backgrounds, let us elaborate 
on the aforementioned jet overlapping. With six quarks in the final state, 
some of the resultant jets will very often be too close to each other to be 
recognizable as coming from different partons. We simulate this as follows. 
We count a final-state parton (quark or gluon) as a jet only if it has a mini- 
mum energy of 5 GeV and lies within the pseudorapidity range —3 < ?7 < 3. 
We then merge any two jets that fall within Sk 5R separation of 0.5; the mo- 



119 



mentum of the resultant jet is the sum of the two momenta. We repeat this 
process iteratively, starting with the hardest jet. Our selection cuts are then 
applied to the (merged) jets that survive this algorithm. 
The signal thus consists of: 

n jets + 27+ {n < 6) 

Hence, all of the SM processes discussed in the previous section yield back- 
grounds to this signal when up to four additional jets are radiated. Now, the 
radiation of each hard and well separated jet suppresses the cross section by a 
factor of order ag — 0.118. Then, since the jj77 backgrounds are already 
quite small after the cuts applied in the previous section (sec Table 5.1), the 
backgrounds with additional jets are expected to be still smaller. 

There exists a potential exception to the last assertion, namely the back- 
ground due to ^77 production. To get an order of magnitude estimate, the 
cross section for tt production at Tevatron Run II is 8 pb [165]. If both of 
the photons arc required to be energetic and isolated, we would expect a 
suppression by a factor of order a^,^, leading to a cross section of the order of 
0.5 fb. This cross section is large enough that we need to consider it carefully. 

If both W bosons were to decay hadronically, then a sufficiently large 
missing transverse energy can only come from a mismcasurement of the jet or 
photon energies. As we have seen in the previous section, this missing energy 
is normally too small to pass our cuts, thereby suppressing the background. 
If one of the W bosons decays leptonically, however, it will yield a sizable 
amount of ^t- This background can be largely eliminated by requiring that 
no lepton (e or /x) is seen in the detector. This effectively eliminates the W 
decays into e or /i or decays to r followed by leptonic r decays; the remaining 
background with hadronic r decays is naturally quite small without requiring 
additional cuts. Considerations such as these lead us to an appropriate choice 
of criteria for an event to be selected: 

1. At least four jets, each with a minimum transverse momentum pTj > 20 
GeV and contained in the pseudorapidity interval of —3 < rjj < 3. Any 
two jets must be separated by 6Rjj > 0.7. As most of the signal 
events do end up with 4 or more energetic jets (the hardest jets coming 

^^For the backgrounds in which one or both of the photons are faked by misidentified 
jets, wc have taken into account the larger combinatoric factor that arises when more jets 
are present to be misidentified. 
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typically from the W boson decays), this does not cost us in terms of the 
signal, while reducing the QCD background significantly. In addition, 
the events with both W^s decaying leptonically are reduced to a 
level of order 10~^ fb by this requirement alone. 

2. Two photons, each one with pT-y > 20 GcV and pseudorapidity —2.5 < 
T]^ < 2.5. The two photons must be separated by at least 6Ry.y > 0.3. 

3. Any photon-jet pair must have a minimum separation of SRj^ > 0.5. 

4. A minimum missing transverse energy > 30 GeV. 

5. The event should not contain any isolated lepton with > 10 GeV 
and lying within the pseudorapidity range — 3 < 77 < 3. 

As in the previous section, the cut on serves to eliminate most of the 
background events with only a fake missing transverse momentum (arising 
out of mismeasurement of jet energies). In association with the lepton veto, 
it also eliminates the bulk of events in which one of the W bosons decays 
leptonically (including the r channel). A perusal of Table 5.3, which summa- 
rizes the major backgrounds after cuts, convinces us that the backgrounds to 
this channel are very small, in fact much smaller than those for the previous 
channel. 

The signal cross section after cuts (but before photon identification effi- 
ciencies) is shown in Fig. 5.4 as a function of the stop and Xi masses. We 
assume that the branching ratio of f — > bWxi dominates in the region of 
parameter space under consideration. The branching ratio of x° — > is 
again taken from Fig. 5.2 assuming that Xi is a pure bino. The signal effi- 
ciency after cuts is about 45%. For small neutralino masses {Tn^o < 50 GeV), 
though, both the photons and the gravitinos {^t) tend to be soft, leading to 
a decrease in the signal efficiency. For small stop masses (as well as for large 
stop masses when the stop-neutralino mass difference is small), on the other 
hand, the jets are soft leading to a suppression of the signal cross section 
after cuts. As one would expect, both of these effects are particularly pro- 
nounced in the contours corresponding to large values of the cross section. 
The additional distortion of the contours for large neutralino masses can, 
once again, be traced to the suppression of the — > 76* branching fraction 
(see Fig. 5.2). 

The non-observation of jjWW^'-f events at Run 1 of the Tevatron 
already excludes the region of parameter space shown in black in Fig. 5.4. 
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Background 


Cross section after cuts 


after 7 ID 




~ 0.003 fb 


~ 0.002 fb 


33^^11 + 2i 


~ 3 X 10-^ fb 


~ 2 X 10-5 fb 


{blrf^, CC77) + 2j 


~ 0.001 fb 


~ 0.0006 fb 


3311 + 2i 


< 0.003 fb 


< 0.002 fb 


tt77, WW jjjj 


< 10-4 fb 


< 10-4 fb 


tt'j'y, WW jjiu, i = e, fi, or T ^ £ 


~ 0.001 fb 


~ 0.0006 fb 


ti'yj, WW jjTu, T ^ j 


< 0.01 fb 


< 0.006 fb 


Backgrounds with fake photons: 






jj{ee -f-f) +2j 


~ 7 X 10-6 fb 


~ 7 X 10-6 fb 


33lU ^ 7) + 2 J 


~ 0.02 fb 


~ 0.02 fb 


33(33 11) + "^3 


~ 0.03 fb 


~ 0.03 fb 


Total 


< 0.07 fb 


- 0.06 fb 



Table 5.3: Backgrounds to tt* jjWW'j'j with both W bosons decaying 
hadronically. The photon identification efficiency is taken to be e-y = 0.8 for 
each real photon. See text for details. 

As in the previous section, in this excluded region at least 3 signal events 
would have been produced after cuts and detector efficiencies in the 100 pb-^ 
of Run I data, with negligible background. In particular. Run I data excludes 
stop masses below about 200 GeV, for neutralino masses larger than about 
50 GeV. 

We show in Table 5.4 the expected maximum stop discovery mass reach 
at Tevatron Run II for various amounts of integrated luminosity. partic- 
ular, with 4 fb-^, a stop discovery can be expected in this channel if rrif < 320 
GeV.^^ If the lightest neutralino is not a pure bino, the reach at large neu- 
tralino masses will be reduced. However, the maximum stop mass reach 
quoted here will not be affected, because it occurs for m^o ~ 50 — 100 GeV; 
in this mass range the lightest neutralino decays virtually 100% of the time 
to 76* due to the kinematic suppression of all other possible decay modes. 

Again, if data from the CDF and D0 detectors are combined, the integrated luminosity 
of the machine is effectively doubled. 

^^For comparison, in the case of minimal supergravity a reach of < 190 GeV can be 
expected in the i bWxf channel with 4 fb~^ at Tevatron Run II [153]. 
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Figure 5.4: Cross section in fb for stop pair production with t — > bW'jG, 
after cuts. Both W bosons are assumed to decay hadronically. The black 
area is excluded by non-observation of jjWW'j'j events in Run I. 

unless its photino component is fine-tuned to be tiny. 

Including a separate analysis of stop production and decay with one or 
more of the W bosons decaying leptonically would yield an increase in the 
overall signal statistics; however, we do not expect this increase to dramati- 
cally alter the stop discovery reach. 

5.7.3 Stop production in top quark decays 

If m^ + m^o < rrit, then stops can be produced in the decays of top quarks. As 
we will explain here, most of the parameter space in this region is excluded 
by the non-observation of stop events via direct production or in top quark 
decays at Run I of the Tevatron. However, some interesting parameter space 
for this decay remains allowed after Run I, especially if the lighter stop is 
predominantly t^. 
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B 


S for a 5cT discovery 




Maximum stop mass reach 


2 fb-i 


0.1 


5 


2.5 fb 


300 GeV 


4 fb-i 


0.2 


6 


1.5 fb 


320 GeV 


15 fb-i 


0.9 


8 


0.53 fb 


355 GeV 


30 fb-i 


1.8 


10 


0.33 fb 


375 GeV 



Table 5.4: Number of signal events (S) required for a 5a stop discovery 
at Tevatron Run II in the jjWWj'y channel and the corresponding 
signal cross section after cuts and efficiencies and maximum stop mass reach. 
We take e-y = 0.80. The number of background events (B) is based on a 
background cross section of 0.06 fb from Table 5.3. 

For mi < nify + rriiy + m^o, so that the stop decays via i —>■ cXi, the 
region in which t ixi is possible is almost entirely excluded by the limit 
on stop pair production at Run I, as shown in Fig. 5.3. A sliver of parameter 
space in which the stop-neutralino mass splitting is smaller than about 10 
GeV remains unexcluded. For rrif > + mv7 + m^o, so that the stop decays 
via t — > bWxi, the signal efficiency in the search for direct stop production 
is degraded for light ncutralinos with masses below about 50 GeV and for 
stops lighter than about 150 GcV. This prevents Run I from being sensitive 
to stop pair production in the region of parameter space in which top quark 
decays to stops are possible with the stops decaying to bWxl, as shown in 
Fig. 5.4. In what follows, we focus on this latter region of parameter space. 

As discussed before, if the lightest ncutralino is mostly bino, the con- 
straints on its mass are model-dependent. The constraints from Tevatron 
Run I are based on inclusive chargino and neutralino production [133, 166] 
under the assumption of gaugino mass unification; the cross section is dom- 
inated by production of xfxf ^^'^ XiX2- ^ the assumption of gaugino mass 
unification is relaxed, then Run I puts no constraint on the mass of Xi- At 
LEP, while the pair production of a pure bino leads to an easily detectable 
diphoton signal, it proceeds only via t-channcl sclcctron exchange. The mass 
bound on a bino Xi from LEP thus depends on the selectron mass [137]. In 
particular, for selectrons heavier than about 600 GeV, bino masses down to 
20 GeV are still allowed by the LEP data. If xi contains a Higgsino admix- 
ture, it couples to the Z and can be pair-produced at LEP via Z exchange. 
For an NLSP with mass between 20 and 45 GeV, as will be relevant in our 
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top quark decay analysis, the LEP search results limit the H2 component to 
be less than 1%. Such a small H2 admixture has no appreciable effect on 
the top quark partial width to ixi- We thus compute the partial width for 
t — > ixi assuming that the neutralino is a pure bino. Taking the lighter stop 
to be ii — i^cos 61 + Ir sin O^, we find, 



4 1 

- sin^ 9f-\ cos^ 9f 

9 * 36 * 



COS^ 9w TTlt 



(5.10) 



where Ej^ = {ml + m| — m?)/2mt. The numerical factors in the square 
brackets in Eq. 5.10 come from the hypercharge quantum numbers of the 
two stop electroweak eigenstates. Clearly, the partial width is maximized if 
the lighter stop is a pure state; it drops by a factor of 16 if the lighter 
stop is a pure II state. In any case, the branching ratio for t ^ tB does not 
exceed 6% for > 100 GeV and m^o > 20 GeV.^« 

The signal from top quark pair production followed by one top quark 
decaying as in the Standard Model and the other decaying to txi, followed 
by the stop 3-body decay and the neutralino decays to ^G, is bbWlV'-ff ^t- 
This signal is the same (up to kinematics) as that from stop pair production 
in the 3-body decay region. As in the case of stop pair production followed 
by the 3-body decay, we expect that the background to this process can be 
reduced to a negligible level. Run I can then place 95% confidence level 
exclusion limits on the regions of parameter space in which 3 or more signal 
events are expected after cuts and efficiencies are taken into account. Using 
the Run I top quark pair production cross section of 6 pb [165] and a total 
luminosity of 100 pb~^, we compute the number of signal events as a function 
of the stop and neutralino masses and the stop composition, assuming various 
values of the signal efficiency after cuts and detector efficiencies. If ti = t^, 
then Run I excludes most of the parameter space below the kinematic limit 
for this decay even for fairly low signal efficiency ~ 20%, as shown in Fig. 5.5. 
If ii — II, on the other hand, the signal cross section is much smaller and 
Run I gives no exclusion unless the signal efficiency is larger than 75%, which 
would already be unfeasible including only the identification efficiencies for 
the two photons; even 100% signal efficiency would only yield an exclusion 
up to rrii ~ 118 GeV. 

At Run II, the top quark pair production cross section is 8 pb [165] and 
the expected total luminosity is considerably higher. This allows top quark 

iSpQj. neutralino masses below m.z., as are relevant here, the branching ratio into jG is 
virtually 100% (see Fig. 5.2), almost independent on the neutralino composition. 
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Stop mass (GeV) 

Figure 5.5: 95% confidence level exclusion limits for top quark decays to stops 
from Run I for various signal efficiencies, in the case that ti = t^, which gives 
the largest event rates. The area to the left of the curves is excluded. The 
case ti = ti gives no exclusion for the signal efficiencies considered and is 
not shown here. The solid line from upper left to lower right is the kinematic 
limit for = 175 GeV. The solid line from upper right to lower left separates 
the regions in which the 2-body FC decay and the 3-body decay of the stop 
dominate. 



decays to txi to be detected for stop and neutralino masses above the Run 

1 bound. For ii = Ir, top quark decays to stops will be probed virtually 
up to the kinematic limit, even with low signal efficiency ~10% and only 

2 fb~^ of integrated luminosity. If ti = t/,, so that the signal event rate 
is minimized, top quark decays to stops would be discovered up to within 
10 GeV of the kinematic limit for signal efficiencies > 20% and 4 fb~^ of 
integrated luminosity (see Fig. 5.6). In this region of parameter space, stops 
would also be discovered in Run II with less than 2 fb~^ via direct stop pair 
production (see Fig. 5.4). 
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Figure 5.6: 5a discovery contours for top quark decays to stops at Run II 
with 4 fb~^ for various signal efficiencies, in the case that ti = ti, which gives 
the smallest event rates. The case ti = would be discovered virtually up 
to the kinematic limit even with 10% signal efficiency, and is not shown here. 
The solid lines are as in Fig. 5.5. 



5.8 Ten degenerate squarks 



Having concentrated until now on the production and decay of the top squark, 
let us now consider the other squarks. In the most general MSSM, the 
spectrum is of course quite arbitrary. However, low energy constraints [167] 
from flavor changing neutral current processes demand that such squarks be 
nearly mass degenerate, at least those of the same chirality. Interestingly, in 
many theoretical scenarios, such as the minimal gauge mediated models, this 
mass degeneracy between squarks of the same chirality happens naturally; in 
addition the mass splitting between the left-handed and right-handed squarks 
associated with the five light quarks turns out to be small. For simplicity, 
then, we will work under the approximation that all of these 10 squarks 
(namely, ul^r, dL,R, cl^r, sl^r and are exactly degenerate. 
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5.8.1 Production at Tevatron Run II 



While the cross sections for the individual pair-production of the cl,r, sl,r 
and hL,R are essentially the same as that for a top-squark of the same mass, 
the situation is more comphcated for squarks of the first generation. The lat- 
ter depend sensitively on the gluino mass because of the presence of t-channcl 
diagrams. Moreover, processes such as uu ulu*j^ or dd dL,RdL,R become 
possible and relevant. Of course, in the limit of very large gluino mass, the 
squark production processes are driven essentially by QCD and dominated 
by the production of pairs of mass eigenstates, analogous to the top squark 
production considered already. In particular, at the leading order, the total 
production cross section for the ten degenerate squarks of a given mass is 
simply ten times the corresponding top squark production cross section. 

Since a relatively light gluino only serves to increase the total cross section 
(see Fig. 5.7), it can be argued that the heavy gluino limit is a conservative 
one. To avoid considering an additional free parameter, we shall perform 
our analysis in this limit. To a first approximation, the signal cross sections 
presented below will scale^^ with the gluino mass approximately as shown in 
Fig. 5.7. 

Like top squarks, the 10 degenerate squarks can also be produced via 
cascade decays of heavier supersymmetric particles. To be conservative, we 
again neglect this source of squark production by assuming that the masses of 
the heavier supersymmetric particles are large enough that their production 
rate at Tevatron energies can be neglected. 

The NLO cross sections for production of ten degenerate squarks includ- 
ing QCD and SUSY-QCD corrections have been implemented numerically 
in PROSPINO [150]. We generate squark production events using the LO 
cross section evaluated at the scale ^ = rriq, improved by the NLO K-factor 
obtained from PROSPINO [150] (see Fig. 5.8), in the limit that the gluino 
is very heavy. The K-factor varies between 1 and 1.25 for niq decreasing 
from 550 to 200 GeV. As in the case of the top squark analysis, we use the 
CTEQ5 parton distribution functions [151] and neglect the shift in the pr 
distribution of the squarks due to gluon radiation at NLO. 



^^That the gluino exchange diagram has a different topology as compared to the (dom- 
inant) quark-initiated QCD diagram indicates that corresponding angular distributions 
would be somewhat different. Thus, the efficiency after cuts is not expected to be strictly 
independent of the gluino mass. For the most part, though, this is only a subleading effect. 
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Figure 5.7: Gluino mass dependence of the NLO cross section for production 
of ten degenerate squarks in 2 TeV pp collisions, from PROSPINO [150]. 
Shown are the cross sections for qq* production normalized to the value at 
large rrig, for common squark masses of 200, 300, 400 and 500 GeV. Produc- 
tion of qq is small at a pp collider and is neglected here; it yields an additional 
3-15% increase in the total cross section at low rrig for this range of squark 
masses. Cross sections are evaluated at the scale /i = rriq. 



5.8.2 Signals in low-energy SUSY breaking 

The decays of the ten degenerate squarks are very simple. As long as they are 
heavier than the lightest neutralino,^^ they decay via q — qxi QlG- The 
signal and backgrounds are then identical to those of the two-body FC stop 
decay discussed in Sec. 5.7.1, and consequently we use the same selection 
cuts. In fact, in view of the tenfold increase in the signal strength, we could 
afford more stringent cuts so as to eliminate virtually all backgrounds, but 
this is not quite necessary. 

The signal cross section after cuts (but before efficiencies) for production 
of ten degenerate squarks in the heavy gluino limit is shown in Fig. 5.9 

^"As before, we do not allow the possibility of cascade decays through other neutrali- 
nos/charginos. Were we to allow these, the channel we are considering would be somewhat 
suppressed, but additional, more spectacular, channels would open up. 
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Figure 5.8: LO (dashed) and NLO (solid) cross sections for the production of 
ten degenerate squarks in the heavy gluino hmit in pp colhsions at Tevatron 
Run I (1.8 TeV) and Run II (2.0 TeV), from PROSPINO [150]. Cross sections 
are evaluated at the scale /i = rriq. 

as contours in the niq-m^o plane. We assume that the squark q decays 
predominantly into gxi- The branching ratio of x? ~^ iG is taken from 
Fig. 5.2 assuming that x\ is a pure bino. Clearly, the effect of the kinematic 
cuts on the signal is very similar to that in the case of the 2-body FC decay 
of the stop. The mass reach, of course, is much larger due to the tenfold 
increase in the total cross section; also, unlike in the case of the stop, the 
2-body decay is dominant throughout the entire parameter space. 

The non-observation of j j77 events at Run I of the Tevatron excludes 
the region of parameter space shown in black in Fig. 5.9. As in our stop 
analysis, in this excluded region, at least 3 signal events would have been 
detected in the 100 pb~^ of Run I data, with negligible background. In 
particular, we estimate that Run I data excludes the ten degenerate squarks 
up to a common mass of about 280 GeV. As in the case of 2-body FC stop 
decays, when the mass sphtting between the squarks and the neutralino is 
too small {i.e., less than about 40 GeV), the jets become soft and the signal 
efficiency decreases dramatically, leaving an unexcluded region of parameter 
space. The D0 search for inclusive — X2 + ^ with X2 ~^ 7Xi [163] yields 
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Figure 5.9: Cross section in ih for production of 10 degenerate squarks in 
Run II in the heavy gluino limit with q —>■ qjC, after cuts. The black area 
is excluded by non-observation of jj'j'j events in Run I. 
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a limit on the production cross section of about 0.5 pb for parent squark 
masses above about 250 GeV. Interpreted in terms of pair production of 
10 degenerate squarks in the large rrig limit with X2 ~^ IXi reidentified as 
Xi — > 7G again increases the signal efficiency by a factor of ~ 2.7 because 
every event contains two photons [163]; the D0 analysis then yields a bound 
on the common squark mass of about 275 GeV, again in rough agreement 
with our result. 

Taking the total background cross section to be 2 fb (see Table 5.1), 
we show in Table 5.5 the expected maximum discovery mass reach for ten 
degenerate squarks at Tevatron Run II for various amounts of integrated 
luminosity. In particular, with 4 fb~^ a squark discovery can be expected 
in this channel if rriq < 360 GeV. Again, if the lightest neutralino is not a 
pure bino, the reach at large neutralino masses will be reduced. However, 
this will have very little effect on the maximum squark mass reach quoted 
here, because this maximum reach occurs for m^o ~ 100 — 150 GeV, where all 

neutralino decay modes other than 7G suffer a large kinematic suppression. 
From Fig. 5.2, it is evident that the neutralino branching fraction into 7G 
will be reduced by no more than 10% in this mass range as long as the 
neutralino is at least 50% bino; such a reduction in the neutralino branching 
fraction will lead to a reduction of only a few GeV in the maximum squark 
mass reach. 



^^Ref. [163] quotes a squark mass bound of 320 GeV for the case rrig <C rng in the 
context of low-energy SUSY breaking; we expect that this is due to the contribution of 
chargino and neutraUno production to the total SUSY cross section in their analysis. 

Again, if data from the CDF and D0 detectors are combined, the integrated luminosity 
of the machine is effectively doubled. 
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B 


S for a 5(7 discovery 




Maximum squark mass reach 


2 fb-^ 


4 


14 


7.0 fb 


345 GcV 


4 fb-^ 


8 


18 


4.5 fb 


360 GcV 


15 fb-i 


30 


31 


2.1 fb 


390 GeV 


30 fb-i 


60 


42 


1.4 fb 


405 GeV 



Table 5.5: Number of signal events (S) required for a 5cr squark discovery at 
Tevatron Run II, assuming production of 10 degenerate squarks in the limit 
that the gluino is very heavy, and the corresponding signal cross section after 
cuts and efficiencies and maximum squark mass reach. Wc take = 0.80. 
The number of background events (B) is based on a background cross section 
of 2 fb from Table 5.1. 
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Chapter 6 

Concluding remarks 



We have examined the possibihty of the existence of a non-minimal Higgs 
sector. In particular, we consider the next to minimal Higgs sector consist- 
ing of two Higgs doublets. There are several motivations to explore those 
scenarios, among them we have: 1) The fact that the Higgs sector is totally 
unknown so far, 2) the unnatural hierarchy predicted by SM between the 
Yukawa couplings of the third family of quarks, 3) the possibility of gener- 
ating either spontaneous or explicit CP violation, 4) the generation of flavor 
changing neutral currents (inspired in the increasing evidence on neutrino 
oscillations) still compatible with the strong experimental limits on them, 
and 5) the possibility that some models with larger symmetries end up at 
low energies in a non minimal Higgs sector as in the case of SUSY models. In 
particular if the SUSY particles are heavy enough, the MSSM Higgs sector 
acquires the form of a constrained 2HDM type II. 

The two Higgs doublet model predicts the existence of five Higgs particles: 
One CP odd Higgs boson (A^), two CP even ones {h^, H^), and two charged 
Higgs particles {H"^). If a charged Higgs were discovered first, it would be a 
clear signature of the presence of a non-minimal Higgs sector. The detection 
of a CP-odd Higgs boson would also yield evidence on physics beyond the 
SM because of the significant deviations of its couplings to SM particles 
respect to the ones of a CP even Higgs. In particular, it worths to mention 
that couplings of A'^VV should be absent, or at least highly suppressed. Of 
course, another possible signal would be the discovery of more than one Higgs. 
Finally, if only one CP even Higgs is discovered and the same experiment is 
unable to detect other Higgs modes, deviations on its couplings from the ones 
predicted by the SM can contain the information on certain underlying Higgs 
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spectrum, in the case of the 2HDM even in the framework of the decouphng 
hmit in which the couphngs at tree level of the lighter Higgs boson are SM- 
like, precision measurements of the couplings could detect such deviations 
via radiative corrections. If a SM-like Higgs is detected, one of the most 
promising couplings to look for possible deviations from the SM behavior is 
the trilinear self coupling h^h^h^. 

The most important sources for production of Higgs bosons in Hadron 
Colliders are the gluon-gluon fusion and qq' production, while in e~^e~ colliders 
the most important sources consists of WW, ZZ fusion and Higgstraghlung. 
Additionally, detection depends on the subsequent decay of the Higgs bosons, 
the dominant decay channels are highly model dependent and the estimated 
rate production should be calculated for specific regions of the free param- 
eters. Finally, other interesting perspectives for Higgs production might be 
provided by 77 fusion in 77 colliders dominated by pW~ or tt loop, and also 
by fj,/! colliders if they run at the Higgs resonance ^/s — mn- 

On the other hand, the presence of FCNC could be an indirect signature 
of the existence of an extended Higgs sector, perhaps the strongest motivation 
to look for these rare processes lies on the evidence of neutrino oscillations. 
Despite processes with FCNC are strongly suppressed by some underlying 
principle still unknown, there is an important region of parameters in which 
those processes are still compatible with experimental constraints. As for the 
quark sector, the FCNC involving the first family are strongly constrained, 
while the bounds on the other mixings are considerably softened. In this work 
we concentrate on FCNC in the lepton sector. If LFV is expected in the neu- 
tral sector of leptons because of the neutrino oscillations, it is reasonable to 
consider that those processes are also present in the charged lepton sector as 
well. Useful constraints involving the vertices C,fiT,^eT,^TT,C,fifi can be derived 
from the {g — 2)^ factor, and the decays fi 67, r fi'-f, and r /jfifi. 
In turn, based on the information gotten for these vertices some upper lim- 
its on the decay widths r(T — > 67), and r(T — > eee) are obtained. Those 
estimations were calculated by assuming ~ llbGeV and uiao ^ rriho. 
Since the most recent estimations of Aa^, still provides an important window 
for new Physics, we obtain from it very interesting lower and upper bounds 
on the mixing vertex at 95% C.L. Specifically, an allowed interval of 
7.62 X 10~^ < < 4.44 x 10~^ was found in a quite wide region of pa- 
rameters. Of course, we should realize that both SM test and experimental 
measurements of still being scrutinized and current results are not 

definitive at all. However, if not severe changes occur in forthcoming experi- 
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ments and/or SM estimations, these constraints could continue being valid at 
least at a lower confidence level. Future improvements on both estimations 
should elucidate this point. 

Based on these constraints, and on the leptonic decays // — >^ 67, r — > 
//7, and r — > /i/iii we got the following conservative and quite general bounds 
of LFV vertices 

< 2.77xl0-i\ 
-1.8 X 10-2 < ^rr < 2.2 X 10-' ; < 0.13 for m^o » m;,o 
\^rr\ < 1-0 X 10-2 ; l^^^l < 0.15 for ttiao ^ niho 

It is remarkable the very strong hierarchy among the allowed values for 
and Moreover, these bounds on ^t-t and ^^i^ are considerably re- 
laxed for certain specific values of m^o. In that case, the room available 
for them is such that they permit either a strong enhancement or a strong 
suppression of the couplings Htt and/or HfxJI for any neutral Higgs boson. 

Furthermore, we estimate upper limits on the decay widths r(r 67), 
and r(r — > eee), finding that they are basically hopeless to look for, in near 
future experiments, at least in the framework of the 2HDM type III with 
heavy Higgs bosons. 

In addition, we see that processes with fiavor changing charged currents 
are intimately related to the processes with fiavor changing neutral currents. 
The same mixing elements that produce the FCNC correct the FCCC in 
the quark sector (respect to the CKM elements) while producing the FCCC 
in the leptonic sector. It worths to point out that the study of the FCCC 
at tree level in the 2HDM (III) depend on a less number of free param- 
eters since only one scalar particle is exchanged, and the Higgs boson in- 
volved does not couple through a mixing angle. For instance, based on 
the decay width F (/x^ — » Uf,e~V^) we can get the upper limit \ien/i^H+\ < 
3.8 X 10-^Gey-^ without any further assumption on the free parameters of 
the model. 

It worths to say that the phenomenological constraints described so far, 
were calculated in the fundamental parametrization of the 2HDM(III) with 
only one VEV. Nevertheless, it is possible to use other paramctrizations for 
the 2HDM type HI in which both VEV's are taken different from zero. In 
that case, the part of the Lagrangian involving the Higgs sector is written 
in terms of the spurious parameter tan/5, in such a way that the Yukawa 
Lagrangian type HI is expressed as the one in the model type I plus some 
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FC interactions, or as the model type II plus some FC interactions. These 
particular parametrizations are equivalent to the fundamental one in which 
one of the VEV is zero and tan {3 vanishes. However, these non trivial bases 
could be useful when we want to compare the model type III with the models 
type I or II. Additionally, these parametrizations facilitates the study of the 
decoupling limit of the 2HDM. By setting a fixed value of tan /9 wc arc Tising 
a particular basis to write the model, in that sense, the parameter tan /3 acts 
as a "gauge fixing" term. 

Now, based on non trivial parametrizations of model III in which models 
type I and II become apparent, we get bounds on the Higgs boson masses, by 
making reasonable assumptions on the FC vertices. Since the contribution 
of the pseudoscalar Higgs boson to Aa^ is negative while the estimated value 
is positive definite at 95% CL, the estimated value of Aa^ impose lower 
bounds on the Pseudoscalar Higgs mass. Specifically, we have taken for 

7/ ( M the geometric average of the Yukawa couplings, and we also utilized 



values one order of magnitude larger and one order of magnitude smaller. 
Taking these three values for the FC vertex we find that the smaller value 

for ( M the more stringent lower bounds for m^o . Additionally, assuming 



ruHO — rriho , we show that in the limit of small (large) tan j3 the lower bound 
of m^o becomes merely m^o ~ m^o for parametrization of type I (II). In 
the case of different scalar masses, there is still a lower asymptotic limit for 
m^o. Notwithstanding, these lower constraints on m^o should be considered 

carefully, since for 77 ( M we can only make reasonable estimations but they 



are unknown so far. 

On the other hand, as a matter of addendum, the production and decays 
of top-squarks at the Tevatron collider is examined in the framework of gen- 
eral low energy supersymmetry breaking models. In models of low-energy 
SUSY breaking, signatures of SUSY particle production generically contain 
two hard photons plus missing energy due to the decays of the two neutralino 
NLSPs produced in the decay chains. Standard Model backgrounds to such 
signals arc naturally small at Run 11 of the Tevatron. We studied the pro- 
duction and decay of top squarks at the Tevatron in such models in the case 
where the lightest Standard Model superpartner is a hght neutrahno that 
predominantly decays into a photon and a light gravitino. We considered 2- 
body fiavor- changing and 3-body decays of the top squarks. The reach of the 
Tevatron in such models is larger than in the standard supergravity models 
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and than in models with low-energy SUSY breaking in which the stop is the 
NLSP, rather than the neutralino. We estimate that top squarks with masses 
below about 200 GeV can be excluded based on Run I data, assuming that 
50 GeV < m^o < mi — 10 GeV. For a modest final Run II luminosity of 4 
fb~^, stop masses up to 285 GeV are accessible in the 2-body decay mode, 
and up to 320 GeV in the 3-body decay mode. 

Top squarks can also be produced in top quark decays. We found that, 
within the context of low-energy SUSY breaking with the stop as the next- 
to-next-to-lightest SUSY particle, the region of parameter space in which 
stop production in top quark decays is possible is almost entirely excluded 
by Run I data if the lighter stop is predominantly right-handed; however, 
an interesting region is still allowed if the lighter stop is predominantly left- 
handed, due to the smaller branching ratio of t ^ ^lXi- Run II will cover 
the entire parameter space in which top decays to stop are possible. 

We also studied the production and decay of the ten squarks associated 
with the five light quarks, assumed to be degenerate. In models of low- 
energy SUSY breaking, the decays of the ten degenerate squarks lead to 
signals identical to those for the 2-body flavor- changing stop decays. The 
cross section for production of ten degenerate squarks at the Tevatron is 
significantly larger than that of the top squark. We estimate that the 10 
degenerate squarks with masses below about 280 GeV can be excluded based 
on Run I data, assuming that m^o < rriq — 40 GeV. For a final Run II 
luminosity of 4 fb~^, squark masses as large as 360 GeV are easily accessible 
in the limit that the gluino is very heavy. The production cross section, and 
hence the discovery reach, increases further with decreasing gluino mass. 
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Appendix A 

Determination of the minima of 
the potential and the mass 
eigenstates 

A.l Minimum condition for the potential 



The potential can be written in terms of the following hermitian gauge in- 
variant operators 



The most general renormalizable and SU (2) x U {1) invariant potential is 
given by 



+X^AB + AelC + X^AD + X^BC + X^BD + XwCD (A.l) 



For the sake of simplicity, we shall use in this appendix the following parametriza- 
tion for the doublets 





(A.2) 
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we will assume from now on, that both VEV can be taken real (i.e. there 
is not spontaneous CP— violation) such that (^3) = vi/\/2 ; (^7) = V2/V2. 
Now, we examine the minimum conditions (tadpoles). 



= 

<l>3=Vl/V2,cf>r=V2/V^ 

where 0j, with i = 1, .., 8 denotes each gauge eigenstate of the scalars defined 
in (A. 2). We obtain the following non trivial equations 




= T3 = ^Ar^s + Xivf + -\6vlv2 - nlvi + -Xsvivi + -XqVivI - -i4v2 

= T4 = (-2/i4 + Xgvl + As^i + Aiot;2'yi) V2 

3 1111 
= T7 = -XyvIvi - iilv2 + X2VI - -i4vi + + -A3V1V2 + -X^V2vl 

= ^8= (2At4 + Xgvl + X^vl + Aiot;2^^i) vi (A.3) 



Now, the mass matrix elements are gotten from the terms quadratic in 
the scalars 



1 9V 



2 d(t)id(t)j 



4>3=v\ I yj2,4>7=V2 1 \/2 



(A.4) 



from (A.4, A.l), and (A. 2) the 8x8 mass matrix becomes 



(Ml, 











Ml 


Ml 








\ 





Ml 








Ml 


Ml 
















Ml 


Ml 








Ml 


Ml 










Ml 


Ml 








Ml 


Ml 




Ml 


Ml 








Ml 













Ml 


Ml 











Ml 
















Ml 


Ml 








Ml 


Ml 




V 





Ml 


Ml 








Ml 


Ml 


I 



(A.5) 



where the matrix elements are given by 
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1 1 
-nl + Xivl + -A5V2 + -X6V1V2 

1 2 Ix Ix 2 Ix 2 

-3/^3 + ^^3ViV2 + -Xev^ + -Xjv^ 

1 2 Ix Ix 2 Ix 2 

-^fJ'i + ^AioVi-ys + -XsV^ + -XgV2 

-fll + Xivf + ^X^vl + ^XqViV2 

Ix 2 Ix 1 2 Ix 2 

-^A9f^2 - -^MoViV2 + -i^i - -Ag-yi 

1 2 Ix 2 Ix 2 Ix 

~2^^ + 1^^^2 + |A6Wi + -Ast'lWs 

-/^i + ?,Xivl + ^A3t;2 + ^A5V2 + ^A6^;i^;2 

-]^XsViV2 - ^XiqvI 



Ill 

-/^i + Ai^i + -A4t'2 + 2^5^2 + 2-^6t'lt'2 

-|^f^2 - ^Aio'yii;2 + ^/^4 - i^gv^ 

1 2 Ix 2 Ix 2 Ix Ix 
"2/^3 + ^'^7'y2 + ^'^e^^l + -X2,ViV2 - -X4V1V2 

-nl + X2v\ + + \^X^vxV2 

-nl + A2'y2 + ^AsVi + ^A7t;it;2 

113 
-nl + 3A2W2 + 2^3t'i + 2-^5^^! + 2'^7^^1^^2 

-XgV2Vx + |Aiofi 

-Ai2 + A2V2 + ^A4Vi + ^As^i + ]^X^VlV2 
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Observe that V2 = 0, is a valid solution for Eqs. (A. 3). Now, if we solve 
using f2 = 0, and replace the other minimum conditions into the mass terms 
(A. 6) we get 



Ml, 



Ml, 



1 



; M^^s = ; Ml, = -X^vl ; MI^ = ; M 



2 

25 



0; M|3 = 2\^vl ■ Mi, = ; M^, = ^X,vl ; M^, = Xsv^ 



Ml 



Ml, 



M^s 



; = -^Ag^;? ; M^g = ; M^^ = -//^ + iA5^;^ 

-/^2 + ^-^5-^1 ; = -III + ^As^i + ^As^i 



We can see that the potential A.l violates CP explicitly by means of the 
terms Mi,, M25, M34, M38, M/^, Afyg in the mass matrix (A. 5), since they 
mix real parts with imaginary parts of the complex neutral fields (p^, see 
Eqs. (A. 2, 2.1). As we shall see below, these terms vanish in the potentials 
(2.4, 2.5, 2.8) avoiding explicit CP violation. Notwithstanding, potentials 
(2.4, 2.6) can exhibit spontaneous CP violation while potentials (2.5, 2.8) do 
not allow spontaneous CP violation either. 

Before continuing with some special cases of the potential, we shall intro- 
duce some useful formulae. If we have a 2 x 2 real symmetric matrix 



a c 
c b 



its eigenvalues and orthonormal eigenvectors are given by 



cos 5 
sin S 



sin 2(5 



a + b±y/{a-bf + Ac^ 
— sin S 



u I -ir^ ki ; 
2c 



^(a-6)' + 4c2 



cos 5 
COS 25 



k2 

{a-b) 
a - bf + 4c2 



(A.7) 



144 



from which we get 



JJ _ f 111 \ — ( ^^^^ ^^^^ 
~ V '^2 J \ —sin5 cos 5 

further, if we start from a doublet of gauge eigenstates 

n = 



and settle the following definitions 



we can check that, owing to the unitarity of the rotation matrix U the fol- 
lowing identity is held 

n^An = H^MH 

and since M is diagonal, we say that the doublet H contains the mass eigen- 
states Hi, H2. They are obtained by rotating the original gauge eigenstates 
through the unitary matrix U; and ki, k2 correspond to the squared masses 
of the mass eigenstates Hi, H2 respectively. These results will be applied 
widely in foregoing issues. 

From now on, we shall examine some special cases of the potential (A.l). 



A. 1.1 The potential with C— invariance 

According to section (2.2.1), if wc demand from the potential to be C— invariant, 
wc ought to settle yu| = Ag = Ag = Aio = 0. Additionally, wc can always make 
a rotation on the two doublets in such a way that only one of them get a 
VEV (see appendix C for details). So we can take (^s) = V1/V2 ; (^y) = 
V2/V^ = 0, without any loss of generality. 

With the assumptions above, the minimum conditions are reduced to only 
two equations 

2x2 2 

III = XiVi ; = ^- 
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replacing them into the mass matrix elements we get 



Ml 



1 2 2 
~2^3 + 4^6^^! 



Mi, 



2X,vl ■ Ml, = ]^XqvI 



-lA + ^^3^? + ]^Kvl ; Mis = -//2 + \KvI + ]^X^vl 



the other terms vanish. Now, we can relabel the matrix in the following way 
1, 2, 3, 4, 5, 6, 7, 8 ^ 1, 2, 5, 6, 3, 7, 4, 8 the new positions of the above elements 
in this new matrix are^ 



/o 

























































Ml 




























Ml, 




























Ml 


Ml 






















Ml 


Ml 




































\o 




















Ml 


) 



(A.8) 



which can be decomposed in the following way 



/ 





\ 






Ml 







Ml 



J 



Ml 



Ml 



Ml Ml 




M|« 



(A.9) 



Additionally, since Ml — Ml we get for the first matrix 



/O 














Ml 








Ml 





Ml 



'2x2 



(A.IO) 



With this decomposition it is easier to diagonalize the matrix (A.8) since 
it is equivalent to diagonalize each of the submatrices given above. Let 



^It is equivalent to relabel the scalar fields ^ 
content of the mass matrix. 



which clearly does not affect the physical 
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us start by diagonalizing the matrix (A. 10) corresponding to the indices 
1, 2, 5, 6. Its orthonormahzed eigenvectors and eigenvalues are given by 



|( 1 ( 1 
|( 1 ( 







)1 

n 



m 



G+ 



= 



thus each eigenvalue is doubly degenerate. Therefore, they correspond to 
charged eigenstates. The eigenvalue rnl+ = correspond to a would be 
Goldstone boson (masslcss scalar) while = —fJ^"^ + ^^5^1 correspond to 

the mass of a charged Higgs boson. According to the eigenvectors obtained 
in (A. 11) the mass eigenstates are given by 



/ ^1 ^ 




( 1 








\ 




/ 01 


\ 


G2 







1 










02 




Hi 










1 







05 








^0 








1^ 




\ 06 


/ 



But taking into account that to define a scalar charged Higgs boson we 
require two degrees of freedom, we realize that the first two scalar fields 
define the charged would be Goldstone boson G^, which is written as a linear 
combination of the scalars 0i, 02. The assigment of both degrees of freedom 
can be done as (0i, 02)^^ — > 0i+i02 = 0^; in the same way (Gi, G'2)'^ Gi + 
iG2 = G"*", something similar happens to (05,06)'^ = 02 ^-nd {Hi,H2)^ = 
H^. Therefore, the relation among the gauge and mass eigenstates of the 
charged scalar fields can be simplified as^ 

G+ 
H+ 



0r 



02+ 



The second matrix in (A. 9) correspond to the indices 3,7. 

IXevf -III + + IX^vj 

applying Eqs. (A. 7) the eigenvalues and orthonormahzed eigenvectors are: 



(Mi, 


Ml ^ 


H 


\Ml 


Ml J 





V 




Ai - -A_ 
2 



+ 



^Noteworthy, this is the opposite procedure respect to the one developed in Sec. (2.2.2), 
in which wc extend the dimension of the representation, see for example Eq. (2.17). In 
this case we are "shrinking" the dimension of the representation. 
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COS a \ 2 / — sm a 

sma J \ cos a ' 

therefore, the diagonahzation process is performed by the following rotation 
of the gauge eigenstates 



7203 -vi\ f 



cos a sm a 

- sin q; cos 0( J ^ V2(f)r — V2 J \ h*^ 
but according to the parametrizations (A. 2) and (2.2) of the Higgs doublets 

7203,7 - ^1,2) /il,2 

obtaining 

cos« sin a \ / hi \ _ f 
- sin a cos a J \ h2 J \ hP 

Finally, we diagonalize the submatrix corresponding to the elements 4,8 
in Eq. (A.9) 

MI4 ML \ f 



ML Mi, J \^ -//i + i (A4 + As) vl 
whose eigenvectors and eigenvalues are 



the first eigenvalue correspond to a neutral, massless would Goldstone boson, 
while the second one is associated to another neutral Higgs boson. 

mco = ; niAo = -fj^l + ^ (A4 + A5) vf 

The mass eigenstates are given by 

1 0\f<P,\_fG' 
1 [ 6, -[ 
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A. 1.2 The potential with a invariance 



As explained in section (2.2.1) a way to avoid spontaneous CP violation 
consists of demanding invariance under the Z-i symmetry $i ^ $i, $2 — ^ 
— $2 (as well as the C— invariance) and correspond to setting /^s = = 
Ag = A7 = As = Ag = Aio = in (A.l). The minimum conditions (A. 3) are 
reduced to two equations: 



vx + \xv\ + A+v^] = 
'^1 [—lA + ^2vl + A+Vi] = 



where A_| 
a) 



I (A3 + A5) . So we have two sets of independent solutions 



^2 ^ A2/X1 - A+/x| 
A1A2 — A^ 



A1A2 — A?_ 



or b) 



First set of solutions 

The decomposed mass matrix, after using the minimum conditions, rela- 
behng and taking into account that Mf^ — M|g, Mf^ — M22, and Mfg = 
M|g becomes 



tot 



Ml 



Ml, 



Ml, Ml, 



'2x2 



MI3 Ml, 
Ml, Ml, 



Ml 



Ml 



Ml Ml 



Ml, = 
Ml = 

Ml = 



1 

2'' 



—^^2^3 ; 



Ml = ^A3'Ui^;2 ; Ml = -^1^2 -^3 ; Ml = ^A37;i^;2 



2Xxvl ; Ml = 2X+V1V2 ; Ml = X_vl ; M|g = -A_'t;i't;2 



Ml 



-\vlX. 



Ml 



2X2VI ; Ml = X_vl 
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where A_ = | (A4 — A3). We first diagonalize the submatrix corresponding 
to 1,5. The eigenvalues and orthonormal eigenvectors are 



^ , 




1 (vf + v^) (//| - A3^i-t;2) 

2 V1V2 



which can be written as 



cos/3 A ^ f -smP \ 1 2> 1 2\ 

smP I ' V cos/3 y 2 ^ 2 ^ 



COS /3 = = ; smp — 



So the mass Higgs bosons, and mass eigenstates are 



G+ \ _ / cos/3 sin/3 
H+ J ~ \-smp cos/3 

ma+ = , = -^A3 {v^ + vj) 



2 

As for the indices 3,7; the eigenvectors and eigenvalues are evaluated from 
(A.7) 



'mHO,h° = ^ivj + A2V2 =•= ^ i^ivf - A2t'|)^ + 4:X\vlvj 

. o 2V1V2X+ 

tan 2a — — — ^ ; — ^ 

(Aiwf - A2WI) 

\ _ / cos a sin a \ / hi 
J Y — sina cos a ) \ h2 

Finally, for the indices 4,8; it is obtained 

cos/3 sin/3 W 04 \ _ / 
-sin/3 cos/3 J \ 08 / V ^° 

mgo = 0, m\o = ^ (A4 - A3) [vl + vl) 
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Second set of solutions 

The mass matrix combined with the minimal conditions is given by: 



Ml 
Ml 






Ml, 
Mi, 





'2x2 



2x2 



-i4 + Ihvf 





Mk 
Mh 



Mir 
Mh 
2Xivl 




Ml Mi, 
Mis Mi, 



-A«2 + + k^bVf 



So in this case the matrix is already diagonal and consequently the mass 
eigenstates coincide with the gauge eigenstates, the elements in the diagonal 
are the eigenvalues i.e. the squared masses of the Higgs bosons. Since Mfg = 
M|g, they correspond to two particles degenerate in mass (i.e. two charged 
particles), in addition, other two degrees of freedom should be associated to 
their corresponding Goldstone bosons. So we get 



^3 
1 



H 



. 



) v'7 
2 . ^2 



H 



. 



^0 



A 



. 



+ o (^3 + A5) vl 



m^o = ; m\o = -^l + ^ (A4 + A5) vl 



A. 1.3 The potential with the global f/ (1) symmetry 

As explained in section (2.2.1) another way to avoid spontaneous CP violation 
consists of imposing invariance under the global C/ (1) symmetry $2 — > 
(as well as the C— invariance) and correspond to setting Ae = Ay = 
As = Ag = Aio = //| = and A3 — A4 in (A.l). In this case, we should assume 
that in general Vx^v^ 7^ 0. The minimum conditions (A. 3) are reduced to two 
equations: 



Ai-^i - //iVi + -^X-iv\vx + -^X,vxv\ - -f4v2 = 
-IJ,lv2 + X2V2 - ^IJ'lvi + ^X3vlv2 + ^AsVa^i = 
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and after reorganizing the matrix elements as 1, 2, 3, 4, 5, 6, 7, 8 — > 1, 2, 5, 6, 3, 7, 4, 8; 
using the minimum conditions, and taking into account that Mfr, — M|g, Mf^ — 
M22, = M|g the decomposed matrix becomes 



ML 



Ml 



Ml, 



Ml, Ml 



'2x2 



MI3 Ml, 
Mir ^77 



Mi, Mis 
Mis Mis 



with 



Ml, 



Mi, 



Ml,- 



-2^2 1. 2 
TT/^a- - 2A3.2 

^fil^^+2X,vl ■ Mi, 



1 

2> 

1 

2' 



1 2 Ix 
"2/^3 + 2'^3'yi'y2 



1 2 

-0/^3 + A5t'lt'2 + X3V1V2 



Ml = 2-"'^ ' = -3/^3 ; Ml, = -/x^-i - -Ast;? 
Mh = ^/4^ + 2A2^;^ M|3 = i^^^ 

Z ^2 Z ^2 

the eigenvalues and the matrix of rotation corresponding to the submatrix 
of indices 1, 5 are 

( G+ \ _ f cos/3 sin/3 \ f <l>t \ 
[h+ J - [ -sin/3 cos/3 J [ cj>t J 

mG+ = , mH+ = 

2 V1V2 

for the indices 3,7, the eigenvalues and the rotation angle read 



tan2Q; 



Xivl + \2v\ + (tan /3 + cot /3) ± i?A 

2wit;2A+ - |/i3 
Ai-yf — A2t'i + iA*i (tan/3 — cot/3) 



Aiw^ — \2v\ + ^/^3 (tan /3 — cot /3) 



+ ^2vi^;2A+ - 
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finally for the matrix 4,8, the rotation angle is P once again 

cos/? sinp \ f \/2(/)4 -vi \ _ f G° 
-sin/3 cos/3 J \ x/2(/)8 - V2 J ~ \ A'^ 

and the eigenvalues are 

1 + 

rriGO = 0, m^o = -/^s 

2 ■i;ii;2 

Finally, it worths to say that the parametrization of the Higgs Hunter's 
Guide (in the CP conserving case i.e. ^ = 0, and A5 = Xq in Eq. (4.8) of Ref. 
[17]), can be gotten from the following associations 



- 


[2X^vl + 2Xs{vl + vl)] 


i4 - 


[2A2^;2 + 2A3 + ^;2)] 


f4 - 


2A5t'it'2 ; Ai ^ Ai + A3 


A2 - 


-> A2 + A3 ; A3 — > A5 — A4 


A4 - 


-> Ae — A4 = A5 — A4 


A5 - 


A4 + 2A3 



where the parameters on left are the ones in our parametrization Eqs. (2.8, 
2.7), and the parameters on right correspond to the ones in the Higgs Hunter's 
Guide. 
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Appendix B 

Rotation of the Yukawa 
Lagrangian in the 2HDM (III) 

We start defining two Higgs doublets with VEV 

^'1,2 = ^'^y j and = v,,, (B.l) 

and writing the Yukawa Lagrangian as 

-£y - V^'^L^\U^n + V^''^L^[D% + ^^^i^^^^ 

+lepton sector + /i.e. (B.2) 

we can make a rotation between the doublets 

V *2 y V -sin^ cos^ / V *2 / 

We shall deal with the quark sector only henceforth, the results that we 
are going to obtain are also valid for the lepton sector with the appropiate 
replacements. In terms of $1, $2, the Yukawa Lagrangian could be rewritten 

as 

-£y = ^'"^"^(cos^Si-sin^S^)^/?^ 

+S'°^L(sm^$i+cos^$2)t/?i^ 

(sin ^$1 + cos ^$2) + h.c. 
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-£y = Q-i(cose^'Vsin^g'°)$iC/°^ 

+01*^ (- sin + cos e^f) 

+g°^ (- sin e?jf.' + COS $2l^°i, + /i.e. 

defining 

V j " V - cos^ J 1^ J ^^-^^ 

the Lagrangian could be written as 

-£y - ^L^f^^U + ^L^'^^D^R 

+^^Li■f^2U% + -^,,i^''^2D% + /..C. (B.5) 

with the same form as the original Lagrangian if we forget the prime nota- 
tion. Consequently, the combined rotations (B.3,B.4) do not have physical 
consequences since it is basically a change of basis. In particular we can 
choose 9 = (3 such that 



= cos + sin ^ ^/l^ ^ ^vf + vl^v 

($2) = - sin /?($;)+ cos /3($'2) = (B.6) 

In whose case we managed to get ($2) = 0. Since this lagrangian contains 
exactly the same physical information as the first one, we conclude that in 
model type III the parameter tan (3 is totally spurious and we can assume 
without any loss of generahty that one of the VEV is zero ^. We shall settle 
9 — P henceforth. 

On the other hand, it is possible to reverse the steps above and start from 
the representation in which ($2) = (the "fundamental representation") and 
make a rotation of the Higgs doublets from which the tan f3 parameter arises. 



^Indeed the rotation should be made in aU the Lagrangians involving Higgs doublets, 
consequently, we have to ensure that the physical content of the potential is not altered by 
the rotation, in appendix (C) we see that the rotation can be performed in the potential 
V defined in Eq. (2.4). 
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If we assume the fundamental representation as our "starting point", we see 
that by making a rotation of the doublets, the mixing matrices are rotated 
such that 

( \ _ ( cos/? -sin^ \ / r]lf'^^'^ \ 

from this rotation we can see that Tj[j^'^^''^ and C,i^'^^'^ depend on the tan/? 

parameter (since rjl^'^^''^ and ^[f'^^''^ come from the fundamental representa- 
tion, they do not depend on tan/3). In the non trivial parametrization the 
VEV's are vi and V2 where tan/3 — V2/V1, sin/3 = V2I \/v\ + f|, and the 
Yukawa Lagrangian in terms of the mass eigenstates can be written in the 
following ways (see section (2.2.3)) 

-£l^ = -i^y (Type I, tan /3) (r]^,^^, tan /3) 

{rf ,e (B.8) 
= -i:Y(Type II, tan/3) + ^^(77^,^^,tan/3) 

+P(r;^,e'',tan/3) (B.9) 

the fundamental parametrization is clearly independent on tan /?, and since 
all these parametrizations are physically equivalent, the parametrizations 
defined by Eqs. (B.8, B.9) cannot depend on tan/3 either. I shall restrict the 
discussion to —£lj but the same ideas are applicable to — £y. In — , it is 
very clear that — £y (Type II, tan/3) depends on tan/3 exphcitly. However, 
as the whole Lagrangian must be independent of it; the mixing matrices 
should depend on tan/3 in such a way that they cancel the dependence on 
this parameter from —£y (Type I, tan/3) . This emphasize the fact that the 
mixing matrices are basis dependent (changing the basis means changing 
tan /3) . 

Let us see explicitly the relation among the parameters in the fundamental 
representation {E}^'^ , rf^'^ , $1, $2, cx) and the ones in the non trivial represen- 
tation (,^^''^,?7'^''^,$'-L,$2,«',/3). First of all, Eqs. (B.3) provides us with the 
relation between ($1,^2) and {^'^,^'2). On the other hand, Eq. (B.4) pro- 
vides the relation among (^^(u,D)fi ^^{u,D)fi^ ^^(u,D),o^~(u,D),oy Notwith-^ 

standing, the most useful relations would be the ones among {^^'^, rf'^) and {^^'^ , rj^'^) 
i.e. the mixing matrices when the Lagrangian is written in terms of mass 
eigenstates. In the non-trivial basis the relation between gauge and mass 
Higgs eigenstates can be taken from (2.9) with the appropiate change of 
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notation. 




performing the rotation (B.3) 



/ cos /? -smp\ f G+ \ 

\ sin/3 cos^ J \H+ J 

f cos a' —sin a' \ f \ 

\ sin a' cos a' J \ J 

( cos(5 - sin /3 \ / (7° \ 

V sin/3 cos/3 ) \ ) 

and using (B.l) we get 



$1 = cos /3$i + sin /3$2 

^ '°'^( (/li + i+lo/v^ ) ^'^^{{h', + v!l\g',)/V2 
<l>t 

{hi + v + igi) /V2 
$2 = - sin + cos /3$2 

^ " '''' ^ ( + i + li) / V2 ) + ^ ( (/.'2 + i + /V2 

<t>t 

(/i2 + ^y2) /\/2 

so the conversion from the gauge eigenstates in the non trivial basis to the 
gauge eigenstates in the fundamental representation is given by 



cos/3 sin/3 \ / (0+)' 

-sin/3 cos/3 ) \ (0+)' 

hi \ _ f cos/3 sin/3 \ / h'l 

h2 J \ — sin /3 cos /3 J \ h'2 

9i \ _ ( cos/3 sin/3 \ / g[ 



92 J \ —smp cos/3 J \ 92 J 

And the relation between the gauge eigenstates and the mass eigenstates in 
the fundamental parametrization, becomes particularly simple. 
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91 

92 

hi 
h2 



h2 



cos P sin P 

— sin P cos P 

cos P sin P 

— sin P cos /9 

cos /? sin P 

— sin /5 cos P 

cos a — sin a 
sin a cos a 



where we have defined 



cos P — sin /3 
sin P cos /3 



h' 



cos a 
sin a' 



hP 



— sm a 
cos a' 



a = a — P 

On the other hand, from (B.4) and (2.44) we get 



H+ 



H+ 



(B.12) 
(B.13) 



?7^'° = cos/3^^'° + sin/5e^'° = cos/3^^'° + sin/5 | ^tIMuTr - ^77^-° 
and after applying the unitarity transformations defined in (2.42) we obtain 



nrf^'Tl, = cos/5 (nv'^^'Tl) + sin/5 - ^ (t^^^'^T^) 



we should note that the unitary matrices Tl,r and Vl,r defined in Eq. (2.42) 
are independent on the basis since they indicate the rotations of the spinors 
from gauge to mass eigenstates, and rotations (B.3,B.4) do not modify such 
spinors. Consequently the transformation T^X^^Tj^ defines the conversion 
of the mixing matrices from the gauge eigenstates to the mass eigenstates, 
therefore 



cos/5^ + sin/5 ( —Mu - -7)^ 
V2 V2 



vi^^ vi^^ ^2 

— Tj H Mu Tj = — 



M, 



u 
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so we obtain finally 



V 



i.e. the matrix t]^ is associated only to the mass of the fermion as expected. 
Similarly from Eqs. (B.4, 2.47, 2.42) 



= -sin/3^'^ + cos/3e'''^ = -sin/3 



7u,o 



Vi Vi 



+ cos 



U,0 



/2 ^ ^ 

Tl^^^tI = -tan/3— Mc/ + sin/3tan/5e^ + cos/5e^ 



R 
'-U 



V2tan/3 



u 



furthermore, Eq. (2.47) is also valid for rj^^'^')'^ by setting V2 — 0, Vi 
V. Using Eq. (2.47) applied to rj^'^ and Eq. (B.4) we get 



= sin/3 



tIMuTr^ 

V 



similarly 



V2sin/3 



Ml 



u 



+ cos /3C^'° 
+ cos/3C^ 



^ = M[/— cos /3 - sin /3^^ 



for the down sector is exactly the same. Summarizing, we get the fol- 
lowing links among the mixing matrices (?7^'^,^^'^) in the fundamental 

parametrization with the ones in the non trivial basis {rj^'^ 

= M^,^^cos/3-sin/3e^'^ ; ^^.^ = ^^^^sin^M^,^ + cos/J^^'^ 



U,D 



V2 



M, 



U,D 



y j^ _ v^tan/3 



U,D 



(B.14) 



and the relations among the gauge Higgs eigenstates and mass Higgs eigen- 
states in the fundamental parametrization, are given by 
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hi \ _ f COS a — sin q; \ / 
h2 J \ sin a; cos a J y hP 

a = a' — (3 

while the analogous relations for the non-trivial bases, are given by Eqs. 
(2.9). 

As a proof of consistency, we shall use the equations B.14 and the Yukawa 
couplings in the non trivial parametrization Eqs. (2.45, 2.46, 2.48, 2.49) in 
order to obtain the Yukawa couplings in the fundamental parametrization 
Eq. (2.39). Let us check first the couplings DDA^ in the parametrization 
in which the model type II becomes apparent, where D refers to the three 
down fermions. From Lagrangian (2.49) and expressions (B.14) 

v2cosp 



igtan/?— ^ i 

^ '^-DMn^.DA^ + -D 

V 2 cos p 



2Mi 



DDA 







w 



Dei.DA' 



V2sin/? 



Md + cos I3i 



D 



1,DA^ 



a/2 



now let us examine DDH^ 
DDH^ 



^ DMdD cos a' + 



2Mw cos f3 



^ D^^ D sin {a' - f3)H^ 



\/2cos/3 



^ DMdD cos a' 



2Mw cos /3 



D 



V^sinp 



Md + cos p^'- 



D sin {a' -(3) 



V2 cos f3 

cos {a' - (5) DH^ + D sin {a' - (3) 

-4^DMdDH^ cos a + -^Di^DH^ sin a 
2Mw V2 
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now DDh^ 
DDh'^ : 



9 



= D 



2Mw cos p 
9 

'2Mwcos(3 

V2COS/9 

qMd 



DMdD sin a'h" + 



V^cos/3 



D^^D cos (a' - /5) 



\/2sin/3 



Md + cos/3C^ 



Dcos(q;'-/3) /i° 



21/^^ cos (3 



(sin a' — sin (3 cos (a' — 



+^L'C^L'cos (a'-/?) /i° 
v2 



9 -n, 



1 ^. 



-DMflD sin (a' - /3) /i° + -^D^^D cos (a' - /5) 

= ^DMoDh^ sin a + -^Di^Dh^ cos a 

2Mh/ V2 

finally for the charged Higgs, in the Lagrangian where model type II becomes 
apparent i.e. (2.45) plus (2.49) the Yukawa coupling is given by 



UDH- 



^^^^ UMuKPlDH^ + -^UrfKPLDH+ 

V2Mw cos/? 
9C0tp-^^ ^^^^ gtanP-^^^, 



\/2M, 



-UMuKPlDH+ - 



w 



\/2M, 



UKM^PrDH^ 



w 



1 



V2g 



^sin/3^ \2Mw 
1 



cos P 



UK 



Mu cos p - sin pf ] KPlDH+ 

PrDH+ 



2Mw 

= -U^^KPlDH+ + UK^^PrDH+ 
= U {K^''PR-fKPL)DH+ 

all these couplings coincide with the ones describe in the Lagrangian for the 
fundamental parametrization, Eq. (2.39). Therefore, with this procedure 
we can check that all Lagrangians generated from Eqs. (2.45, 2.46, 2.48, 
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i.L,. ■^y([/) <^y(D)' °^y{c/) '^Y(Z3)' ■^y(i/) '^y(D)' °^y(!7) 
'^y(D)' coincides with the Lagrangian (2.39) if we take into account the ex- 
pressions (B.14). 

The results expressed by Eqs. (B.14), show us that the value of the flavor 
changing vertices is basis dependent, though the value of the couplings are 
basis independent as it must be. The mixing angles between Higgs gauge 
eigenstates and Higgs mass eigenstates are also basis dependent as expected. 
The transformations (B.3, B.4) reveals an 5*0(2) global symmetry of the 
model type HI. This is like a "global gauge invariance of the 2HDM type IIP 
in which tan (3 fixes the gauge. 

In that sense, we can realize that the models type I and II have a re- 
markable difference respect to the model type 111. since it is well known that 
the former two ones are highly dependent on the tan/5 parameter while the 
latter is not. In writing the parametrization <£(/ we can see easily the rea- 
son: — £y (Type II) clearly depend on tan/3 but when we add the mixing 
parameters, they acquire the precise values to cancel such dependence. In 
other words, the model type II does not have mixing parameters at the tree 
level to absorb their tan/? dependence. We can see the difference from the 
point of view of symmetries, the 2HDM is constructed in such a way that we 
make an exact "duphcate" of the SM Higgs doublet. 

^-^-^ 

These doublets have the same quantum numbers and are consequently in- 
distinguishable (at least at this step). Owing to this indistinguibility we can 
perform the rotation 

/ \ _ / cos^ sin^ W ^ 

V *2 y V -si^^ cos^ y V *2 y 

over an arbitrary angle 6 without any physical consequences (it is in fact a 
change of basis). It means that the model is invariant under a global SO (2) 

transformation of the "bidoublct" ( $1 $2 ) • However, it is very common 
to impose a discrete symmetry on the Higgs doublets ($1 — > $1, $2 —$2) 
or a global U {1) symmetry ($1 — > $1, $2 — ^ to prevent dangerous 

FCNC. In that case, we are introducing a distinguibility between the dou- 
blets, because they acquire very different couplings to the fermions (models 
type I and II). Of course, we could have defined the symmetry in the oppo- 
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site way ($1 — > — $1, $2 ^2) but once we have chosen one of them, we 
cannot interchange $1 ^ $2 anymore, without changing the physical con- 
tent. Such fact breaks exphcity the SO (2) symmetry of the "bidoublet". On 
the other hand, it is precisely this symmetry what allows us to absorb the 
tan P parameter, and since models type I and II do not have that symmetry, 
we are not able to absorb it properly. 
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Appendix C 



Rotation in the Higgs potential 



C.l Transformation of the parameters in the 
potential 

Let us start from an arbitrary parametrization in which both VEV are in 
general different from zero. The most general renormalizable and gauge 
invariant potential read 

Vg = - JJlB' - JgC' - JllD' + + X2B'^ + XsC'^ + X^D'^ 

+X5A'B' + XeA'd' + XsA'D' + XjB'C' + XgB'D' + XwC'D' (C.l) 

where we have defined four independent gauge invariant hermitian operators 

A' ^ , B' ^ d'^^ + $;t$;) ^ Re ($;t$',^ , 

the doublets and the VEV are denoted as 

^'1,2 = ( Sti-^l' ] = ( f.St!::L.y ] ^nd = v[^, (C.2) 



V2 



In appendix B, wc have seen that in the Yukawa Lagrangian type III, we are 
able to make a rotation of the doublets as 

$1 \ / cos^ sin^ \f^[\ .p.. 

*2 y V -sin^ cos^ y V *2 / 
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without changing the physical content of the Lagrangian. However, we must 
demonstrate that such rotation can be carried out in the potential without 
changing its physical ^ontent either. In order to show it, we shall calculate 
the way in which /Ij, parameters transform under this rotation. First, we 
calculate the way in which the operators A', B' , C", D' transform. Taking 
into account Eq. (C.3) we get 

A = = (^$1 cos e-^l sin 9^ ($1 cos 9-^2 sin 9) 

= cos^ 9 - $|$2 cos 6* sin 6* - $^$1 cos 6* sin ^ + $|$2 sin^ 9 

= cos^ 9-2cos9 sin 9 |^£i^i±^2^ j + $t ^^^2 q 

= Acos^9 + BsiTi^9 -sm29d 

similarly, we obtain the transformation for the other operators, the results 
read 

A' = lcos^^ + Ssin2^-Csin2^ 
B' = Asrn^ 9 + B cos^ 9 + C s\n29 

C' = -Asi's\29--Bsin29^Ccos29 
2 2 

D' = D 



^ A? cos^ 9 + P sin" 9 + C'' siu^ 29 + ]-AB siu^ 29 

-2 AC sin 29 cos^ 9 - 2BC sin^ 9 sin 29 

= ^sin^^ + S2cos^^ + C^sin2 2^ + -lSsin2 2e 

2 

\2AC sin 29 sin^ 9 + 25(5 sin 2^ cos^ 9 
= i(l2 + ^2^sin2 2^ + a2cos2 2^-^lBsin2 2^ 

+-ACsin4e - -SCsin4e 
2 2 
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A'B' = (^^A^ + -C^^ sm"" 20 + AB {cos'' e + sm^'e) 
+ (^AC - BC^ sin 29 cos 29 

A'd' = -A:^sm29cos'^9- -Psiri^ 9 sm29-d^sm29 COS29-- AB sin i9 
2 2 4 

+AC (4 cos^ ^ - 3) cos^ 9 + BC{A cos^ 9 - l) sin^ 9 
A'D' = ADcos^ 9 + BDsin^ 9- CD sin29 



B'C' = -A?sin29sin^9 --B'^ sin29cos^ 9 + -C"^ sin A9+-AB sin A9 
2 2 2 4 

+AC (cos 29 + 2 cos^ 9) sin^ 9 + BC (cos 29-2 sin^ 9) cos^ 9 
B'D' = ADsin^ 9 + BDcos^ 9 + CD sin29 

C'D' = ^ (^15 - sin 29 + c5 cos 29 (C.4) 
Now, we can build up a new parametrization of tiie potential such that 

Vg = -^xlA-^xlB-^llc-^llD + x^A? + X2p + XiC'^ + XAD'^ 

+X^AB + XqAC + AglS + A7-BC + X^BD + AioC5 (C.5) 

in order to find the values of //j, in terms of /Ij, Aj, we use the Eqs. (C.l), 
and (C.4) to write e.g. the coefficient proportional to the operator A, and 
these terms are compared with the term proportional to the operator A in 
Eq. (C.5) obtaining 

—jji\A = (— /^i cos^ ^ ~ A*2 si^^ ^ ~ A*3 sin 9 cos 9) A 

therefore, the coefficient /if is related to the parameters Aj in the following 
way 

by the same token, the other sets of /Xj, Aj parameters are related to the 
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Aj parameters in the following way 



= ^//^ cos 6* + //2sin^6' + -//3sin26' 

A*2 — (^1 sin^ + Ji^ cos^ 9 — -Jll sin 29 

A*3 — ("/^i sin 2^ + Jll sin 2^ + Jll cos 2^) 

2 ~2 



Ai = ^AiCos^^ + A2sin^e+^ p3 + A5)sin2 2^ 

+^ (^Ae cos^ 6' + A7 sin^ 9^ sin 261 
A2 = ^Aisin^e + A2C0S^^ + ^ (^A3 + A5) sin22e 

(^Ae sin^ 9 + X7 cos^ ^) sin 2^^ 
A3 = ^ (Xi + A2 - A5) sin^ 29 + A3 cos=^ 2^ + ^ (^Ay - Ag) sin 4^^ 



A4 — A4 



A5 = Q (^Ai + A2 - A3) sin^ 29 + A5 (cos^ ^ + sin^ 9) 

+^ (Ar-Ae) sin 4^^ 

As = 2 (^2 sin^ ^ - Ai cos^ 9^ sin 2^ + ^ (^Aa + A5) sin 49 

+Xe (4 cos^ ^ - 3) cos^ 9 + X7 (cos 2^ + 2 cos^ 9) sin^ 9 
As = ^As cos^ 61 + 3^9 sin^ 9 + ^ Aio sin 26*^ 
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At = 2 (^Aa cos^ ^ - Ai sin^ 6^ sin 2^ - ^ (^Aa + A5) sin 49 

+Xe (4 cos^ ^ - 1) sin^ 9 + X7 (cos 2^-2 sin^ 9) cos^ ^ 
Ag = ^Ag sin^ ^ + Ag cos^ ^ — -Aio sin 2^ 

Aio = (|(^Ag - Ag) sin2e + Aiocos2^) (C.6) 
C.1.1 Tadpoles 

From now on, we shall consider the potential with invariance under charge 
conjugation, i.e. 114 = = Xg = Aio = (see section 2.2.1). In that case the 
tadpoles are given by 

/I 113 

^ = ( -A«i^i - 2^*3^2 + Xivf + -XsVivj + -X^vivl + -Xevfvi 

1 \ f 1 1 

+ ^^7V2 ) ^1 + ( - + X2V2 + -X3VIV2 

1 1 3 \ 

+ -X5VIV2 + -^K^l + -A7V2V1 j /l2 



/I 1 

= ( - ^fJ'lv2 + Xivf + -X3V1V2 

13 1 
+^X5Vivl + -XQvfv2 + ^Ay^a 

f 1 1 

^7=1 -IJ'lv2 - -/^3^^1 + X2VI + -X3VIV2 

1 1 3 \ 

+^X5vlv2 + -Xevf + -Xjvlvi j (C.7) 

these tadpoles coincide with the minimum conditions as we see from Eq. 
(A. 3) and applying fi4 = Xg = Xq — Aio = 0. Now, wc find the relation 
among the tadpoles in both parametrizations by using Eq.(C.7), and the 
third of Eqs. (B.12). 
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T^hi + T7/12 = T3 {h'l cos 9 + sin 9) + T7 (— /i-i sin 9 + cos ^^2) 
= (T3 cos ^ - T7 sin 9) h[ + (T3 sin 9 + T7 cos ^) /i^ 
= T;,h[+T^h', 

from which we see that the tadpoles in both parametrizations are related 
through the rotation 

\ - f cos9 -sm9 \ f T3 \ , . 

n ) = [ sm9 cos9 )\T, ) ^^"^^ 

As a proof of consistency, we can check that from (C.7), and from (C.6) the 
following relation is gotten after a bit of cumbersome algebra 

T3 cos 6* - T7 sin 6* = (^-Jilv'i - + ^v'l + ^^^svWi 

+ 2'^5ViV2 + -XqV^ V2 + -X7V2 

T3 sin 9 + Ttcos9 = ( -/x^t;^ - -fxlv[ + Aa^^a^ + -Ast;^ t;^ 

+l'^5v'M + ^A6< + Ai) (C.9) 

where 

v\ \ _ f cos 6* — sin ^ \ f 

sin 9 cos 9 J \ V2 

relates the VEV between both parametrizations^. The terms on right of Eqs. 
(C.9) are precisely the tadpoles in the "prime parametrization" , then the 
relations given by Eqs. (C.8) are held as expected. Therefore, tadpoles are 
preserved by the rotation. 

C.1.2 Higgs boson masses 

Another important proof of consistency is to verify that both parametriza- 
tions predict the same masses for the Higgs bosons. We shall use once again, 

^Observe that, such rotation keeps invariant the quantity Vi+V2= v'l + v'2 = 
as it should be. 
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the potential with C— invariancc Eq. (2.4). In a general paramctrization, the 
minimal conditions can be taken from A. 3 by using /i4 = Ag = Ag = Aio = 0, 
as it corresponds to the Lagrangian in Eq. (2.4). Those conditions are re- 
duced to 

/I 1 1 3 1 

II2V2 = (^-^fJ'SVi + X2V2 + ^X3vlv2 + ^X5vlv2 + ^Xev^ + ^XyVivI 

the mass matrix is obtained from (A. 5), and (A. 6) by using once again /i4 = 
As = Ag = Aio = 0. Let us start with the matrix elements corresponding to 
ttT'HO , frihP ■ If we assume that both VEV are different from zero and utilize 
the minimum conditions, we obtain the following mass matrix 



Ml Ml 
Ml Ml 



(C.IO) 



with 



Ml = ^ (2//2v2 + 8A1V? + 2>XqvIv2 - X7vl) 

Ml = ~l4, + ^'^7^^2 + + X2,ViV2 + A5V1V2 

Ml = -^[2iJilvi^%X2vl-XQvl^?>X7vlvi) (C.ll) 

4^2 

For the sake of simplicity, we just show that the determinant of this matrix 
(i.e. the product of the squared masses), coincides for two parametrizations 
connected by a transformation like (C.3). The mass matrix in any other 
parametrization with both VEV different from zero, have the same form as 
(C.IO, C.ll) but replacing //f — > /If, Aj — > A^. It is a fact of cumbersome 
algebra to demostrate that 



MlMl - [Mlf = MlMl - [mI) 



2 



this demostration is carried out by taking into account the relations (C.6) 
among the parameters in both bases. In a similar fashion, we can show that 
the eigenvalues coincide in both bases. Therefore, the mass Higgs bosons 
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arc equal in both parametrizations as it must be. Finally, if the angle of 
rotation is chosen such that one of the VEV is zero, (e.g. V2 = 0) in one of 
the bases, then the minimum conditions and mass matrix elements become 
much simpler (see Sec. A. 1.1), and the equality is easier to demonstrate. 

By the same token, we can check that for the other Higgs mass matri- 
ces the determinants and eigenvalues are invariant under the transformation 
(C.3). Showing that the observables are not altered by this change of basis. 
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